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ABSTRACT 


A  key  element  of  the  U.S.  Navy’s  transition  to  an  electric  naval  force  is  an 
Integrated  Power  System  (IPS)  that  provides  continuity  of  service  to  vital  systems  despite 
combat  damage.  In  order  to  meet  subsequent  survivability  standards  under  a  reduced 
manning  constraint,  the  IPS  system  must  include  a  fault  tolerant  control  scheme,  capable 
of  achieving  automated  graceful  degradation  despite  major  disruptions  involving 
cascading  failures.  Toward  this  objective,  online  model-based  residual  generation 
techniques  are  proposed,  which  identify  explicitly  defined  faults  within  a  stochastic  DC 
Zonal  Electrical  Distribution  System  (DC  ZEDS).  Two  novel  polynomial  approaches  to 
the  design  of  unknown  input  observers  (UIO)  are  developed  to  estimate  the  partial  state 
and,  under  certain  conditions,  the  unknown  input.  These  methods  are  shown  to  apply  to  a 
larger  class  of  systems  compared  to  standard  projection  based  approaches  where  the  UIO 
rank  condition  is  not  satisfied.  It  is  shown  that  the  partial-state  estimate  is  sufficient  to 
the  computation  of  residuals  for  fault  diagnosis,  even  in  such  cases  where  full-state 
estimation  is  not  possible.  In  order  to  reduce  the  complexity  of  the  system,  a  modular 
approach  to  Fault  Detection  and  Isolation  (FDI)  is  presented.  Here,  the  innovations 
generated  from  a  bank  of  Kalman  filters  (some  of  them  UIOs)  act  as  a  structured  residual 
set  for  the  stochastic  DC  ZEDS  subsystem  modules  and  are  shown  to  detect  and  isolate 
various  classes  of  faults.  Certain  mathematical  models  are  also  shown  to  effectively 
identify  input/output  consistency  of  systems  in  explicitly  defined  fault  conditions. 
Numerical  simulation  results  are  based  on  the  well-documented  Office  of  Naval  Research 
Control  Challenge  benchmark  system,  which  represents  a  prototypical  U.S.  Navy 
shipboard  IPS  power  distribution  system. 
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EXECUTIVE  SUMMARY 


The  incorporation  of  an  Integrated  Power  System  (IPS)  in  the  DDG  1000 
Destroyer  represents  a  first  step  in  a  transfonnation  to  an  Electric  Naval  Force  [72]. 
Moreover,  the  U.S.  Navy  is  presently  assessing  alternative  power  and  control  system 
architectures  for  the  Next  Generation  Integrated  Power  System  (NGIPS)  [395],  These 
systems  must  be  survivable  and  increasingly  capable  of  brokering  propulsion  and  pulsed 
power  demands  while  providing  ships  service  power  continuity  to  vital  mission  systems 
[91].  In  an  era  of  asymmetric  warfare  and  reduced  manning  constraints,  the  NGIPS  must 
provide  automated  continuity  of  service  control  to  vital  systems  despite  major  disruptions 
involving  cascading  failures.  Toward  this  objective,  online  physics-based  fault  detection 
and  isolation  (FDI)  techniques  are  proposed  and  are  expected  to  operate  within  a  fault 
tolerant  control  (FTC)  framework.  The  methods  seek  to  identify  explicitly  defined  faults 
within  a  DC  Zonal  Electrical  Distribution  System  (DC  ZEDS)  portion  of  an  IPS.  Here, 
mathematical  models,  together  with  sensors  providing  real-time  data,  are  used  to  validate 
the  condition  of  the  system.  The  results  of  this  research  demonstrate  various  methods  of 
generating  fault-indicating  signatures  for  classes  of  sensor,  actuator,  and  plant  faults.  In 
the  case  of  system  actuator  FDI,  several  estimation  techniques  based  on  a  polynomial 
realization  are  developed.  Numerical  simulation  results  are  based  on  the  well- 
documented  Office  of  Naval  Research  Control  Challenge  benchmark  system,  which 
represents  a  prototypical  U.S.  Navy  shipboard  IPS  power  distribution  system. 
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I.  INTRODUCTION 


In  keeping  with  the  vision  set  forth  by  the  Chief  of  Naval  Operations,  Naval 
warships  in  the  Twenty-First  Century  must  have  the  ability  to  “change,  adapt  and 
transform”  to  meet  new  threats.  Indeed,  Sea  Power  21  defines  a  Navy  capable  of 
constantly  improving  its  ability  to  “project  offensive  power,  defensive  assurance,  and 
operational  independence  around  the  globe”  [111].  To  achieve  the  vision,  the  CNO 
outlines  three  fundamental  pillars;  among  them  is  a  continuing  transfonnation  of 
improved  technology. 

Achieving  these  goals  require  warships  with  a  higher  degree  of  survivability  to 
enable  decisiveness,  sustainability,  responsiveness  and  agility,  which  are  the  four 
fundamental  qualities  in  Sea  Power  2 1 .  From  an  operational  point  of  view,  survivability 
is  the  ability  to  avoid  and  withstand  a  hostile  environment.  From  a  system-level  point  of 
view,  it  can  be  thought  of  as  a  “fight  hurt”  or  “fight-through”  capability  -  that  is,  the 
warship’s  ability  to  maintain  continuity  of  service  despite  major  disruptions  and 
cascading  failures  in  the  presence  of  catastrophic  events. 

Within  this  framework,  an  automated  fault-tolerant  control  scheme  is  essential  to 
a  warship’s  fight-through  capability.  Of  particular  importance  in  achieving  a  graceful 
degradation  of  warship  systems  under  these  conditions  is  a  design  strategy  that  integrates 
combat  systems  with  advanced  hull,  mechanical  and  electrical  capabilities.  According  to 
the  U.S.  Navy  Survivability  Design  Handbook  for  Surface  Ships,  developing  affordable 
survivability  measures  under  a  distributed,  dynamically  interdependent  system  of  systems 
framework  will  “constitute  the  hallmark  of  the  Navy's  new  approach  to  reliably  automate 
damage  control  functions  while  operating  with  fewer  crewmembers”  [343]. 

Take  for  example  an  Integrated  Fight  Through  Power  (IFTP)  system,  which  seeks 
to  maintain  power  continuity  resulting  from  a  calamitous  exogenous  event.  According  to 
numerous  Naval  Research  Advisory  Committee  (NRAC)  reports  and  naval  power 
systems  and  automation  experts,  there  is  a  growing  need  to  develop  cost-effective  ways 
to  solve  power  automation  issues  by  taking  into  account  the  following  aspects: 
[390]  [391]  [91]  [264]  [72]  [76] 
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1.  Automation  will  be  essential  to  ensure  survivable  and  effective  surface 
combatants  in  the  warfare  environment  [271]. 

2.  Personnel  costs  comprise  over  50%  of  operating  and  supporting  costs, 
therefore  it  is  necessary  to  reduce  the  number  of  shipboard  war  fighters. 
Smart  Ship  demonstrated  that  technology  and  procedural  changes  can 
improve  ship  operations  [272]. 

3.  Well-trained  crew  members  are  decision-makers  that  require  real-time  data  to 
make  informed  decisions.  No  longer  able  to  rely  on  a  well-staffed  crew  for 
data  and  unable  to  react  to  rapid  and  complex  nonlinear  system 
interdependencies,  an  automation  strategy  must  instead  respond  to  cascading 
system  faults  and  failures  [392]. 

4.  Life  Cycle  Costs  (LCC),  which  include  research  and  development,  can  be 
reduced  by  implementing  simulation-based  design  environments  -  especially 
in  cases  where  revolutionary  technologies  are  being  developed  -  where  no 
historical  cost  data  exists  [273]. 

5.  Advanced  distributed  simulations  and  simulation-based 
design/manufacturing  to  assess  new  technologies  and  ship  design  options  is 
encouraged  [273] [392], 

Given  this  trend,  it  is  most  desirable  and  advantageous  to  employ  an  on-line 
analytically  redundant,  robust  Fault  Tolerant  Control  (FTC)  scheme  designed  to  detect 
unpermitted  deviations  from  a  system’s  expected,  acceptable  or  standard  condition  and 
perform  actions  as  necessary  to  maintain  continuity  of  service  -  in  this  case,  power 
distribution.  Herein,  the  first  of  a  two-step  FTC  process  is  presented  using  robust  model- 
based  Fault  Detection  and  Isolation  (FDI)  methods.  Characteristics  of  faults  such  as 
time,  location,  size,  and  relative  influence  are  detennined  through  residual  generation  and 
analysis  methods  developed  herein.  The  methods  are  applied  to  a  DC  Zonal  Electrical 
Distribution  System  (DC  ZEDS)  testbed  designed  to  resemble  dc  power  distribution 
systems  currently  being  researched  for  installation  on  future  all-electric  warships  for 
reasons  to  be  discussed. 
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A.  ROADMAP  TO  AN  ELECTRIC  NAVAL  FORCE 

The  NRAC  report  Roadmap  to  an  Electric  Naval  Force  identified  several  major 
motivations  for  the  navy  to  shift  away  from  a  separated  power  and  propulsion  platform 
and  towards  an  all-electric  warship  [271]  [390] [391]  [392]: 


•  Electric  weapons  and  advanced  sensors  for  superior  firepower  range  and 
resolution  require  far  greater  pulse  power  than  currently  available. 

•  Electric  propulsion  and  auxiliaries  possess  superior  mobility,  stealth  and 
endurance. 

•  Integrated  Power  System  (IPS)  allows  provision  for  common  electric 
power  system  for  real-time  power  allocation,  reconfigurability,  and 
superior  survivability. 

•  Electric  ship  platform  provides  support  for  off  board  weapons  and  sensors 
for  superior  reach  and  warfighter  sustainment. 


These  reasons  among  many  others  make  transitioning  to  an  all-electric  drive 
warship  as  significant  a  change  as  shifting  from  steam  to  gas  turbine  and  possibly  even 
from  sail  to  steam.  With  this  change  comes  the  challenge  of  managing  far  greater  power 
requirements  -  on  the  order  of  tens  of  megawatts,  which  is  only  the  beginning 
[11] [72] [90].  Already  committed  to  fielding  an  electric  warship,  the  U.S.  Navy  is 
currently  building  DDG  1000  as  the  baseline  for  future  DC  zonal  electric  power  systems 
[271].  Again,  survivability  requirements  and  current  trends  in  power  electronics 
technology  obviate  the  need  for  a  power  automation  scheme  capable  of  providing 
continuity  of  service  under  temporal  and  spatial  bursts  of  natural  or  hostile  disruption 
[390].  A  Fault  Tolerant  Control  scheme,  then,  provides  a  framework  upon  which  to  build 
a  desirable,  realizable,  and  tractable  solution. 
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B.  FAULT  TOLERANT  CONTROL 

A  fault  can  be  described  as  any  deviation  from  nominal  system  behavior  that  can 
negatively  impact  overall  desired  system  performance.  A  mathematical  model-based 
approach  to  the  detection  of  a  fault,  therefore,  is  sought.  From  a  system  automation  point 
of  view,  a  representative  Fault  Tolerant  Control  architecture  can  be  depicted  as  in  Figure 
1.1  [44]. 


/V 


/ 


Figure  1.1:  Fault  Management  Architecture 


Notice  the  fault  management  architecture  contains  an  execution  level,  which 
resembles  common  model-based  state  feedback  control  system  architecture  and  a 
supervisory  level,  which  depicts  two  major  processes  -  diagnosis  and  controller  re¬ 
design.  The  plant  includes  modeling  uncertainties  defined  as  a  vector  of  unknown  inputs 
d  .  A  vector  of  faults  / ,  cause  the  plant  to  deviate  from  a  nominal  condition. 
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Assuming  linearity,  a  plant  model  (depicted  in  Figure  1.1)  can  be  expressed  in 
continuous  time  as 


x  =  Ax  +  Bu  +  Ed  +  Gf 
y  =  Cx  +  Fd  +  Hf 


where  state  vector  x  e  M",  input  u  e  R'",  and  output  y  e  W  and  matrices  A,  B,  C,  E, 
F ,  G  and  H  are  of  appropriate  dimensions.  The  terms  Ed  and  Gf  model  uncertain 
(unknown)  inputs  and  faults,  respectively,  in  the  actuator(s)  and  plant  model.  In  the 
output  equation,  Fd  and  F[f  represent  uncertain  (unknown)  inputs  and  faults  in  the 
sensors.  Note  that  matrices  A,  B  ,  and  C  represent  the  nominal  system  where  all  faults 
and  uncertainties  are  assumed  to  be  contained  within  /  and  d  .  It  is  further  assumed  that 
matrices  E ,  F ,  G  and  H  are  constant  and  known,  while  /  and  d  are  typically 
unknown. 

Looking  at  the  reconfiguration  level  in  Figure  1.1,  the  main  objective  of  the  on¬ 
line  Fault  Diagnosis  process  block  is  to  detennine  if  the  system  is  subjected  to  a  fault  and 
its  time  of  occurrence.  Then,  all  subsequent  actions  should  focus  on  characterizing  the 
fault(s)  to  the  greatest  extent  possible.  These  actions  should  include  estimating  the  fault’s 
location,  identification,  and  severity.  If  a  history  of  fault  activity  is  needed  for 
conditioned-based  maintenance  or  scenario  reconstruction,  then  developing  fault  statistics 
would  be  advantageous.  This  step  is  represented  as  the  output  to  the  Fault  Diagnosis 
process  block  and  is  annotated  by  / . 

The  Controller/Plant  adjustment  process  block  then  receives  the  fault  information 
and  performs  appropriate  controlling  action  as  depicted  in  the  figure.  This  action  might 
be  to  adjust  the  parameters  of  the  controller  or  reconfigure  the  controller  altogether.  If  a 
fault  becomes  a  failure  in  the  nominal  plant,  then  it  might  be  necessary  to  operate  in  a 
reduced  capability,  damage  control  reconfiguration  mode  in  order  to  maintain  continuity 
of  service.  If  this  is  the  case,  then  it  might  be  advantageous  to  have  the  Controller/Plant 
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adjustment  process  block  initiate  a  reconfiguring  of  the  nominal  plant  and  associated 
mathematical  models  in  the  Fault  Diagnosis  block. 

Some  work  has  been  done  at  the  execution  level  within  a  Fault  Tolerant  Control 
framework.  To  a  limited  degree,  fault  tolerance  can  be  achieved  under  a  robust  or 
adaptive  control  design.  In  Robust  Control,  passive  fault  tolerance  is  achieved  by  a  fixed 
controller  that  is  designed  to  satisfy  pre-defined  performance  specifications.  Adaptive 
control  methods,  on  the  other  hand,  are  considered  active  fault  tolerance  measures  since 
controller  parameters  are  adjusted  to  changing  plant  parameters  (possibly  caused  by 
faults)  as  a  function  of  time.  In  both  cases,  however,  the  classes  of  faults  are  restricted  in 
the  sense  that  they  must  be  “well-behaved”  according  to  their  fault  definitions  and 
mathematical  fonnulations.  Adaptive  control,  for  example,  performs  well  only  for  linear 
systems  with  slowly  changing  parameters  [44].  Neither  method  attempts  to  diagnose  a 
broader  range  of  faults  nor  do  they  attempt  to  change  the  controller  or  plant  structure  to 
mask,  compensate,  correct,  re-configure,  re-design  or  otherwise  “tolerate”  faults.  This 
work  looks  beyond  the  restricted  class  of  faults  to  generate  residuals  at  the  supervisory 
level. 


C.  OVERVIEW  AND  CONTRIBUTIONS 

As  described  in  Figure  1.1,  the  flow  in  the  reconfiguration  level  of  a  representative 
FTC  architecture  begins  with  Fault  Diagnosis.  It  is  the  intent  of  this  work  to  provide  the 
output  of  the  Fault  Diagnosis  process  block  -  that  is,  the  necessary  fault  information  for 
further  processing  to  obtain  the  Navy’s  goal:  achieve  automated  fight-through  capability 
despite  major  disruptions  and  cascading  failures.  Common  error  masking  and  error 
recovery  methods  are  poorly  suited  for  a  system  description  under  this  framework 
[390]  [392] . 

The  goal  of  this  work  is  to  develop  Fault  Diagnosis  strategies  applied  to  the  DC 
Zonal  Electrical  Distribution  System  (DC  ZEDS)  by  generating  model-based  residuals 
from  estimates  of  the  partial  state.  The  following  chapters  will  establish  a  methodology 
for  establishing  the  DC  ZEDS  models,  estimation  methods,  and  Fault  Detection  and 
Isolation  (FDI)  methods. 
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Chapter  II  details  the  modules  of  the  DC  ZEDS  system  to  include  how  some  of 
the  nonlinear  dynamics  are  averaged.  Since  a  deductive  approach  will  be  taken  to 
identify  the  effects  of  a  system  fault,  a  close  observation  of  fault  effects  will  be 
monitored.  It  will  be  shown  that  output  sensor  measurements  are  selected  to  guarantee 
observability  and  are  consistent  with  algebraic  relationships  and  conservation  of  energy 
laws.  Chapter  III  then  explains  the  methodology  for  choosing  estimation  methods, 
including  several  novel  polynomial  unknown  input  observer  methods.  An  explanation 
for  using  linear  methods  on  nonlinear  models  will  also  be  addressed.  Chapter  IV 
identifies  common  mode  failures  and  identifies  classes  of  faults  through  simulation,  then 
provides  a  methodology  for  choosing  the  Fault  Diagnosis  schemes  to  developing  fault 
indications  in  the  residuals.  Finally,  Chapter  V  provides  a  summary  of  results. 

Overall,  the  contributions  of  this  work  can  be  categorized  into  2  main  areas.  The 
first  is  in  a  novel  polynomial  approach  to  the  design  of  unknown  input  observers  (UIO) 
with  particular  attention  paid  to  partial-state  and  input  estimation.  The  second  area  is  in 
the  FDI  techniques  applied  to  the  DC  ZEDS  application. 

A  structural  residual  set  -  predicated  upon  physics-based  mathematical  model 
estimation  -  is  developed  to  decouple  the  effects  from  a  defined  set  of  faults  in  order  to 
determine  fault  indications.  The  unknown  input  observer  (UIO)  developed  in  this 
research  will  be  part  of  this  task. 

Even  though  UIOs  are  used  for  Fault  Detection  and  Isolation  (FDI),  little 
application  research  has  been  done  to  take  advantage  of  input  estimation,  which  is  shown 
to  enhance  analytical  redundancy.  In  addition,  it  will  be  shown  that  the  proposed 
polynomial  based  approach  can  be  applied  to  a  more  general  class  of  systems,  compared 
to  the  standard  projection  based  approach  of  [175]. 

A  polynomial  approach  to  a  general  UIO  for  state  and  input  estimation  for  both 
discrete  and  continuous  time  systems  is  presented  for  the  purpose  of  fault  detection  and 
isolation.  Unlike  UIO  methods  that  are  based  on  a  subspace  projection  of  the  whole  state 
onto  a  subspace  of  known  state  components  (such  as  [73],  [80]  and  [175]),  this  robust 
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estimation  method  is  based  on  an  estimate  of  the  partial-state  representation  of  the 
system.  This  is  sufficient  to  the  computation  of  a  residual,  even  in  cases  when  full-state 
estimation  is  not  possible. 


8 


II.  DC  ZONAL  ELECTRICAL  DISTRIBUTION  SYSTEM 


Among  all  shipboard  distributed  systems,  it  can  be  argued  that  the  power 
distribution  system  is  the  most  vital.  If  not  for  its  tightly  coupled  interdependence  with 
every  shipboard  distributed  system  (including  propulsion),  it  would  be  in  terms  of 
potential  vulnerability.  Indeed,  an  all-electric  drive  warship  depends  absolutely  upon  its 
power  system  to  provide  quality  and  continuity  of  service  while  operating  in  a  hostile 
environment.  And  with  future  warships  potentially  requiring  eight  to  ten  times  the  power 
generated  on  today’s  warships,  it  will  only  increase  its  stock. 

The  prototypical  power  electronics  based  dc  zonal  electrical  distribution  system  is 
uniquely  designed  as  an  efficient  high  power  density,  fast  transient  response,  tightly 
coupled  finite  inertia  system.  The  sophisticated  and  robust  power  converter  control 
elements  provide  nearly-ideal  load  regulation  and  transient  perfonnance  (tens  of 
milliseconds)  with  a  high  degree  of  fault  tolerance.  The  system  is  also  characterized  as 
having  a  high  degree  of  automation  and  nearly  instantaneous  reconfiguration  capabilities. 
For  example,  a  dc/dc  buck  converter  -  which  is  the  Ship’s  Service  Converter  Module 
(SSCM)  in  the  representative  DC  ZEDS  system  herein  -  is  capable  of  maintaining  a 
constant  output  voltage  regardless  of  input  disturbances.  From  the  output  side,  this 
attribute  is  highly  desirable  for  power  quality-sensitive  loads. 


A.  ONR  INTEGRATED  POWER  SYSTEM  TESTBED 

The  nonproprietary,  well-characterized,  reduced-scale  ONR  IPS  Testbed  DC 
ZEDS  model  serves  as  an  ideal  resource  for  initial  validation  of  electric  warship  power 
distribution  analysis  methods  and  control  schemes.  Through  sponsorships  from  the 
National  Science  Foundation  (NSF)  and  the  Office  of  Naval  Research  (ONR),  the 
MATLAB  Simulink  based  model  was  developed  by  a  collaboration  of  academic 
communities  from  Purdue  University  and  the  U.S.  Naval  Academy.  The  distribution 
system  topology  is  also  representative  of  the  Naval  Combat  Survivability  Testbed,  a 
reduced-scale  hardware  testbed  designed  to  achieve  the  same  purpose. 
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For  the  purpose  of  this  research,  the  DC  ZEDS  model  -  in  both  average  value  and 
detailed  model  forms  -  is  used  to  apply  novel  estimation  and  residual  generation 
techniques  for  robust  model-based  fault  diagnosis.  When  considering  the  topology,  the 
DC  ZEDS  system  is  spatially  designed  to  improve  survivability  and  redundancy  by 
splitting  the  DC  power  distribution  into  a  two-bus  system  -  one  on  the  port  side,  the  other 
on  the  starboard  -  each  being  fed  by  a  dedicated  power  supply.  The  power  distribution 
system  is  subdivided  into  3  zones  (numbered  forward  to  aft)  that  are  separated  by  water¬ 
tight  compartments  as  shown  in  Figure  II.l. 


Figure  II.l:  DC  Zonal  Electrical  Distribution  System  Topology 


Within  each  zone,  two  Ship’s  Service  Converter  Modules  (SSCM)  are  fed  from 
the  port  bus  and  starboard  bus,  respectively.  Both  of  their  regulated  voltage  outputs  are 
compared  and  all  the  zonal  load  current  demand  is  supplied  from  the  bus  with  the  higher 
voltage.  The  converter  modules  have  the  capability  to  share  power  requirements  in  droop 
mode,  but  is  not  considered  in  this  study.  Opposing  diodes  on  the  lines  leading  to  the 
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load  prevent  one  bus  from  providing  power  to  the  other.  The  load  in  each  zone  is 
identical,  consisting  of  a  Ship’s  Service  Inverter  Module  (SSIM)  feeding  an  ac  load  bank. 

If  either  power  supply  fails,  the  remaining  bus  will  supply  the  entire  load  without 
interruption.  Current  limits  on  the  SSCMs  alleviate  faults  at  the  oring  (diode)  function 
while  some  faults  within  the  modules  are  mitigated  through  robust  controls.  It  should  be 
noted  here  that  the  control  scheme  improves  the  robustness  of  the  system,  but  this  can 
have  a  masking  effect  on  detecting  and  isolating  faults.  Further  discussion  on  this  will  be 
addressed  in  a  later  chapter. 

As  stated  before,  there  are  two  time-domain  simulation  models  considered  in  this 
work:  the  Average  Value  Model  (AVM)  and  the  more  accurate  Detailed  Model.  Power 
electronics  modeling  experts  have  noted  that  a  so-called  ‘detailed’  model  is  unfortunate 
because  the  term  is  rather  arbitrary,  but  have  also  documented  its  sufficient  accuracy  as  a 
truth  model  for  initial  system  analysis  and  control  architecture  simulations.  For  the 
purpose  of  this  dissertation,  which  is  consistent  with  the  literature,  ‘detailed’  refers  to  a 
simulation  in  which  semiconductor  switching  action  is  included,  even  if  only  on  an  on  or 
off  basis.  By  comparison,  the  nonlinear  average  value  model  refers  to  simulations  where 
the  switching  is  represented  on  an  average  value  basis  resulting  in  state  variables  constant 
in  the  steady  state.  Although  the  average  value  model  removes  the  switching  dynamics 
present  in  the  detailed  model,  both  models  are  comparable  in  every  other  respect,  which 
includes  parameter  values,  control  architecture,  dc  steady  state,  and  transient  response. 

By  comparison,  the  detailed  model  is  the  more  representative  simulation  of  the 
DC  ZEDS  system,  but  the  average  value  model  provides  the  basis  for  the  state  space 
models  used  for  Fault  Diagnosis.  The  remaining  part  of  the  chapter  describes  these 
modules  in  more  detail  starting  with  the  detailed  model,  followed  by  the  nonlinear  state 
space  models  of  the  average  value  model.  It  should  be  noted  here  that  the  red  markings 
on  each  circuit  diagram  represent  sensor  measurements.  By  convention,  a  red  sensor 
marking  with  a  negative  slope  indicates  a  current  sensor,  while  a  positive  slope  indicates 
a  voltage  sensor  (always  measured  with  respect  to  ground).  Finally,  a  frequency  analysis 
of  the  two  models  is  presented  and  a  method  is  proposed  to  generate  residuals  based  on 
the  AVM  state  space  models  and  measurement  data  from  the  Detailed  Model. 
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B.  MODULES 

1.  Power  Supply 

a.  Detailed  Model 

Each  15  kW  power  supply  is  fed  from  a  dedicated  3-phase  ac  source  that 
may  vary  from  480  -  560  V  line-to-line  rms  at  60  Hz  and  provides  500  V  dc  output  to  the 
bus.  The  power  supply  consists  of  an  isolating  transformer  (modeled  as  part  of  the 
leakage  inductance),  a  controlled  3-phase  bridge  rectifier,  and  dc  link  capacitance, 
inductance  and  resistance  as  shown  in  Figure  II. 2. 


For  the  simulations  herein,  the  uncontrolled  rectifier  mode  is  considered 
with  a  zero  firing  angle  (a  =  0) .  The  output  voltage  vout  is  regulated  by  a  buck 

converter/controller  scheme  connected  to  the  output  of  the  rectifier.  Parameters  for  the 
system  are  shown  in  Table  II.  1 . 
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Parameter 

Description 

Value 

Units 

h 

Source  Commutating  Inductance 

1.0 

m  H 

Rdc 

Rectifier  Line  resistance  (dc  link  resistance) 

0.03 

Q 

Ldc 

Rectifier  line  inductance  (dc  link  inductance) 

0 

H 

cr 

Rectifier  output  (dc  capacitor) 

500.0 

af 

°>g 

Input  ac  voltage  frequency 

377 

rad  /  sec 

RL 

M out 

Buck  controller  line  resistance 

0.03 

Q 

Rout 

Buck  controller  line  inductance 

3.0 

fflH 

Cout 

Buck  controller  output  capacitance 

500.0 

af 

Table  II.  1:  Power  Supply  Parameters 


b.  Average  Value  Model 

The  average  value  model  for  the  power  supply  is  developed  in  [220]  and 
makes  the  assumptions  that  the  rms  amplitude  of  the  ac  source  voltages  E  is  constant, 
only  one  commutation  occurs  at  a  time,  and  the  dc  load  current  idc  is  constant  (or  at  least 
near-constant).  The  average  value  model  makes  several  reference  frame  transformations 
of  the  input  voltage  to  obtain  the  time-dependent  peak  voltage  of  the  input,  y]2E .  A 
more  detailed  derivation  of  the  transformation  is  shown  in  Appendix  B.  Since  the  output 
waveform  of  the  bridge  rectifier  circuit  is  periodic  over  a  ;r/3  interval  of  the  angular 
position  of  the  source  voltage  6,r ,  the  average  dc  voltage  can  be  expressed  as 


(0-2) 
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With  further  substitution  and  simplification,  the  average  dc  voltage  from 
the  rectifier  can  be  expressed  as 


v 


d 


3^3 


V2 E 


cos  a 


(0.3) 


and  the  average  dc  current  id  is  a  state  in  the  average  value  model  as  shown  in  Figure 
II.3. 


71 


I 


out 


Figure  II.3:  Power  Supply  Average  Value  Model 

From  the  model  in  Figure  II. 3,  the  power  supply  rectifier  and  buck 
converter  are  dynamically  coupled  through  a  current  and  voltage  source  dependent  on  the 
duty  cycle  d  ,  which  is  the  control  input  used  to  regulate  the  output  voltage  vout . 
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c.  Power  Supply  Control 

The  power  supply  controller  design  is  shown  in  Figure  II. 4  with  its 
associated  parameters  listed  in  Table  II. 2. 
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Figure  II.4:  Power  Supply  Control 

The  nonlinear  controller  allows  for  voltage  and  current  regulation  as  it 
calculates  the  control  signal  d  that  detennines  the  on/off  state  of  the  switch.  For  the 
average  value  model,  d  is  the  duty  cycle  of  the  input  state.  The  measured  output  voltage 
vout  is  compared  with  the  reference  voltage  Vref  and  the  error  signal  is  fed  through  a 

proportional  plus  integral  (PI)  controller.  To  provide  a  fast  response  to  changes  in  load, 
the  load  current  iout  is  added  to  the  output  of  the  PI  controller  and  is  multiplied  by  the 

output  of  nonlinear  stabilizing  control  (NSC),  which  has  a  scaling  effect  on  the  influence 
of  the  output  voltage  and  current.  The  resulting  reference  current  signal  is  bounded  and 
compared  with  the  measured  current  i,  in  the  current  regulator.  The  output  of  the 

‘-'out 

current  regulator  is  used  for  the  hysteresis  modulator  so  that  the  actual  current  closely 
tracks  the  measured  current. 
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Parameter 

Description 

Value 

Units 

Kef 

Reference  voltage  (command  output  voltage) 

500.0 

V 

i 

max 

Maximum  command  current 

40.0 

A 

kp 

PI  controller,  proportional  gain 

1.0 

Kt 

PI  controller,  integral  gain 

100.0 

Kih 

Current  regulator,  integral  gain 

100.0 

n 

NSC  exponent 

1.0 

Krr 

Hysteresis  error  bandwidth 

1.0 

A 

V 

i,max 

PI  controller,  integral  gain  output  max 

100.0 

V 

v/,min 

PI  controller,  integral  gain  output  min 

-100.0 

tf 

J  max 

NSC  Transfer  Function  output  max 

1000.0 

tf . 

J  mm 

NSC  Transfer  Function  output  min 

400.0 

Table  II.2:  Power  Supply  Control  Parameters 


d.  Nonlinear  State  Space  Model 

Looking  at  the  power  supply  average  value  model  circuit  in  Figure  II. 3, 
the  inputs  to  the  system  can  be  defined  as 


Ml 

4ie 

u2 

di, 

^ out 

n3 

dvr 

U4\ 

_  hM  _ 
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If  the  controller  is  considered  external  to  the  model,  the  control  signal  d 
can  be  viewed  as  an  external  input  dynamically  coupled  to  states  of  the  system.  In  this 
manner,  nonlinear  terms  are  defined  as  external  inputs  u2  and  w3 .  The  model  is  then 

expressed  as  a  linear  state  space  model  driven  by  inputs  w,  and  u4  with  a  decoupled 
nonlinear  term  as  shown 


where 


u. 

U, 

x  =  Apsx  +  BXps 

t 

+  B2  ,ps 

z 

_^4  J 

w3 

(0-4) 
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1 

C~ 
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0  0 

0  0 


0 

0 

-R, 


out 

1 


C„ 
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3y/3 

0 

'  0 
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-1 

0 

0 

0 

,  B ;  = 

0 

0 

’  2  ,ps 

0 

1 

0 

-1 

L out 

0 

0 

c 

^  out  _ 

and  the  output  equations  are  chosen  as 


y  = 


1 

0 

0 

o' 

0 

0" 
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0" 
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0 

0 

0 

0 

0 
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0 

0 
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0 
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-1 

0 

1 

0 

0 

0 
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1 

+ 

2 

0 

0 

1 

0 

0 

0 
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0 
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U3 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

1 

0 

0 

-R, 

0 
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0 

-1 

0 

0 

0 

0 

1 

0 

— 

(0-5) 
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Although  the  dynamic  model  has  only  four  states,  we  included  a  larger 
number  of  observations  for  the  purpose  of  redundancy  and  fault  detection.  Figure  II. 5 
shows  a  simulation  of  the  power  supply,  plohing  the  input  and  output  voltage  for  both  the 
detailed  and  average  value  models.  Therein,  the  input  voltage  is  shown  to  have  the  same 
dc  value  between  the  two  models,  but  the  detailed  model  shows  a  small  amount  of  ac 
content  and  switching  dynamics.  The  output,  however,  is  nearly  identical. 
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Figure  II.5:  Comparison  of  Detailed  and  Average  Value  Models  of 

the  Power  Supply 


2.  Bus 

The  bus  model  connects  the  power  supply  to  the  SSCM  in  each  zone.  It  is  a 
medium  length  Nominal-T  circuit  transmission  model  with  switches  separating  the  three 
zones  as  shown  in  Figure  II. 6.  The  series  impedance  and  shunt  admittance  parameters 
are  equivalent  in  each  zone  and  are  shown  in  Table  II. 3. 
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Zone  1 


Zone  2 


Zone  3 


Figure  II.6:  Bus  Model 


Parameter 

Description 

Value 

Units 

R 

Bus  line  resistance 

34.0 

mfi 

L 

Bus  line  inductance 

2.26 

//H 

C 

Bus  shunt  line  capacitance 

0.5 

//F 

G 

Bus  shunt  line  admittance 

10.0 

Q  1 

Rs 

Bus  zone  short  line  resistance 

1.0 

Q 

Table  II.3:  Bus  Model  Parameters 


Included  in  the  bus  model  is  a  shunt  resistance  Rs  used  to  model  bus  shorts. 
During  no-fault  operations,  the  shunt  resistance  is  infinite  to  model  an  open  line.  While 
not  part  of  the  power  electronics  modules,  the  bus  model  has  proven  to  be  quite  an 
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important  link  between  modules.  Generalized  fonns  of  the  bus  model  are  developed  such 
that  n  zones  can  be  considered  on  larger  scale  systems. 

a.  Six-state  Bus  Model 

Consider  three  operating  modes  of  the  bus  system  where  Mode  1 
represents  the  power  supply  providing  power  to  Zone  1  only,  Mode  2  is  where  Zone  1 
and  Zone  2  are  energized,  and  Mode  3  is  where  all  3  zones  are  energized.  The  mode 
switch  configurations  for  each  zone  are  described  in  Table  II. 4. 


Switch  .s'. 

Switch  s2 

Mode  1 

Open 

Closed 

Mode2 

Closed 

Open 

Mode  3 

Closed 

Closed 

Table  II.4:  Bus  Zone  Switch  Configuration 


It  is  desired  to  express  the  bus  model  in  matrix  equation  form.  From 
Kirchhoff  s  conservation  laws,  the  following  equations  are  derived  from  the  circuit  in 
Figure  II. 6: 


R  G  1  c 
x  =  -—x  +  -Qmvz+-Su 


i  =  —Q,x - v.  -  Mu 

Z  ZZ'Hl Z  - 


V.  =  —l,  +  Z 

G  ~ 


z  =  —K, 

~  C 


(0.6) 

(0.7) 

(0.8) 

(0.9) 
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where 


_ 

^bz{ 

u\ 

U2 

_ 

l bz2 

w3 

^ bz3 

U4\ 

ps 

x  =  [x(  x2  xj  =  bus  current  through  zones  , 
z  =  [zj  z2  z3  ]?  =  voltage  across  zone  shunt  capacitance , 
v.  =  vZj  vz_  J  =  zone  bus  voltage ,  and 
i  =  i  i  i  =  zone  shunt  admittance  line  current . 

Z  z2  Z3  J 


From  equations  (0.6)  through  (0.9),  it  can  be  shown  that  the  six-state 
model  can  be  expressed  as 


R 

L 


(r  1  rA 

G  X\rn  " 
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~  M 

CG 

RSCG 
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(0.10) 


with  output  equations 
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(0.11) 


x 

z 


- 

06x4 

X 

1 

+ 

--M 

[zj 

G 

G 

— -M 

G 

where 


G  =  G  +  — , 
R 


(0.12) 


M  = 


10  0  0 
0  10  0 
0  0  10 


(0.13) 


and 


5  = 


0  0  0  0 
0  0  0  0 
0  0  0  1 


(0.14) 


The  effective  conductance  G  shown  in  (0.12)  includes  the  shunt  resistance  term  used  in 
modeling  a  bus  short.  Under  no-fault  conditions,  the  shunt  resistance  line  in  Figure  II. 6  is 
modeled  as  an  open  circuit  ( Rs  =  oo) .  The  model  also  includes  the  three  mode 

configurations  with  the  matrix  Qm  defined  as 
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T 0 

0 

-f 

i 

-1 

0 

,  m  =  1 

(mode 

1) 

[-1 

0 

0 

r  o 

1 

-fi 

0 

-1 

0 

,  m  =  2 

(mode 

2) 

L-l 

0 

0 

r  o 

1 

-1] 

i 

-1 

0 

,  m  =  3 

(mode 

3) 

[-! 

0 

0 

Note :  Q2  =  Q[  and  Q3  =  Ql. 


State  ordering  assignments  for  the  port  and  starboard  busses  are  also 
considered  in  the  model  using  the  transfonnations 


"1 

0 

o' 

"0 

0 

f 

0 

1 

0 

and  Tstbd  = 

0 

1 

0 

0 

0 

1 

1 

0 

0 

to  account  for  opposite-side  state  numbering.  To  demonstrate  the  bus  model,  a  baseline 
scenario  is  identified  for  all  fault  diagnosis  experiments.  In  further  chapters,  any  induced 
system  faults  will  reference  the  baseline  system  dynamics  for  analysis.  Table  II. 5  defines 
the  nominal  plant  configuration. 


Zone  1 

Zone  2 

Zone  3 

Load  on  (sec) 

0.4 

0.3 

0.4 

Stbd  SSCM  vref  (V) 

400 

400 

390 

Port  SSCM  vref  (V) 

380 

410 

410 

Table  II.5:  DC  ZEDS  Baseline  Scenario  Parameters 
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Herein,  both  busses  are  in  Mode  3  with  all  six  SSCMs  on  at  the  start  of  the 
simulation.  Each  zone  supplies  an  identical  constant  power  load  of  5.65  kW.  From  the 
Oring  Function  dynamics,  current  is  supplied  from  the  starboard  bus  in  zone  1  and  from 
the  port  bus  in  zones  2  and  3.  Figures  II. 7  -  II. 9  show  an  example  of  the  bus  current 
distribution  of  the  DC  ZEDS  baseline  scenario  for  both  the  detailed  and  average  value 
models. 
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Figure  II.7:  DC  ZEDS  Baseline  Scenario  Bus  Current  Distribution,  Zone  1 

From  the  results  in  Figure  II. 7,  it  can  be  seen  that  the  load  in  each  zone 
requires  about  12  Amps.  For  Zone  1,  the  port  bus  (bottom)  is  directly  connected  to  the 
port  power  supply,  which  shows  current  flowing  through  to  Zone  2  at  t=0.3  seconds.  An 
additional  12  Amps  is  added  to  the  port  bus,  which  is  being  sent  to  Zone  3  at  t=0.4  sec. 
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The  starboard  bus  (top)  is  at  the  end  of  the  line  since  it  is  directly  connected  to  Zone  3  at 
the  opposite  end.  Here,  Figures  II. 7  -  II. 9  show  that  the  starboard  bus  only  feeds  Zone  1. 


Starboard  Bus  Current,  Zone  2 


Figure  II.8:  DC  ZEDS  Baseline  Scenario  Bus  Current  Distribution,  Zone  2 
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Starboard  Bus  Current,  Zone  3 


Time  (sec) 

Figure  II.9:  DC  ZEDS  Baseline  Scenario  Bus  Current  Distribution,  Zone  3 


b.  Three-state  Bus  Model 

The  three-state  bus  model  uses  the  same  parameter  values  as  the  six-state 
model  as  shown  in  Table  II. 3  with  the  exception  of  completely  passive  line  impedance  in 
each  zone.  This  is  equivalent  to  setting  the  line  inductance  L  to  zero.  From  a  modeling 
and  simulation  point  of  view,  the  states  to  describe  the  dynamics  of  the  bus  current  can 
no  longer  be  found  by  integration,  but  by  algebraic  relationship.  The  six-state  model  is 
essentially  a  double-cascaded,  medium  length  Normal-T  circuit  transmission  line  model 
described  in  [323].  The  energy  storage  elements  are  very  small  and  equate  to  transient 
dynamics  three  orders  of  magnitude  faster  than  that  of  the  power  electronics.  It  also 
implies  that  most  of  the  power  is  transferred  from  the  power  supply  to  the  SSCMs  with 
very  little  energy  storage  in  the  bus,  which  is  desired.  Reducing  the  number  of  states  may 
improve  the  realism  of  the  model,  but  might  incur  unresolved  algebraic  loops  in  the 
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model.  It  will  be  shown  that  a  three-state  model  can  be  derived  and  is  demonstrated  to  be 
useful  for  input  estimation  on  the  unknown  input  observers. 

To  begin  the  model  derivation,  the  same  matrix  definitions  are  used  except 
this  time  the  three  modes  are  expressed  as  subspace  projections  of  the  full  mode. 
Consider  the  bus  model  in  Figure  II. 6  and  the  switching  configurations  in  Table  II. 4.  If 
the  bus  is  operating  in  the  full  mode  (Mode  3),  it  can  be  represented  as  an  orthogonal 
basis  in  three  dimensions  as 


h  ’ 


with  each  column  (x,y,z  from  left  to  right)  representing  the  bus  current  through  Zone  3, 
Zone  2,  and  Zone  1  (from  left  to  right)  respectively.  The  numbering  convention  in 
reverse  order  is  an  unfortunate  result  of  having  2  power  supplies  entering  the  DC  ZEDS 
system  from  opposite  zones.  Each  lesser  mode  can  be  expressed  as  a  subspace  basis  of 
q3 .  For  example,  the  Mode  2  basis  spans  the  y-z  plane  as 


<h  = 


0  0 
1  0 
0  1 


while  the  Mode  1  basis  spans  the  x-z  plane  as 


<h  = 


1  0 
0  0 
0  1 


27 


In  a  physical  sense,  the  Mode  2  basis  describes  opening  the  switch 
between  Zone  2  and  Zone  3,  which  sets  the  current  in  Zone  3  to  zero  and  the  dynamics 
follow.  Likewise  the  Mode  1  basis  is  like  opening  the  switch  between  Zone  1  and  Zone 
2,  which  sets  the  current  in  Zone  2  to  zero.  Now  define  the  mode  projection  matrix  as 


P„  = 


q 


T 

m 


where  m  =  1,  2,  3  denotes  the  mode.  Removing  the  inductance  in  the  bus  line  changes 
the  dynamics  of  the  bus  and  shunt  line  current  Kirchhoff  law  equations  to 


x  =  —  Qv=+  —  Su  (0.15) 

R  R 


and 


4  =  —Qx  -  Mu  . 


(0.16) 


With  some  substitution  and  algebra  on  equations  (0.8),  (0.9),  (0.15)  and  (0.16),  the  three- 
state  model  is  expressed  as 


z  =  — 


CG 


T,QC,  +  j-I 


z  -- 


CG 


(: TxQDl+M)u 


(0.17) 


Where  G  ,  G  ,  Tx ,  C ,  R  ,  Rs,  and  M  are  the  same  as  in  the  six-state  model  and 
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-1 
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The  corresponding  output  equations  (defined  similarly  to  the  six-state  model)  are 
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where 
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The  properties  of  the  3 -state  model  show  that  the  D  matrix  is  full  rank, 
which  turns  out  to  be  important  for  full  input  estimates  from  unknown  input  observers  as 
will  be  explained  in  the  following  chapter. 

3.  Ship’s  Service  Converter  Module 

The  purpose  of  the  Ship’s  Service  Converter  Module  (SSCM)  is  to  step  the  bus 
voltage  down  to  an  acceptable  power  quality  and  rating  for  zonal  loads.  As  shown  in 
Figure  II.  1,  the  DC  ZEDS  system  has  2  SSCMs  in  each  zone  -  one  receiving  input  power 
from  the  port  bus,  and  the  other  from  the  starboard.  All  six  SSCMs  are  identical  in 
topology  and  parameters.  Rated  at  8  kW,  the  SSCM  accepts  a  500  V  input  (that  may 
vary)  and  provides  a  clean,  400  V  output  and  a  maximum  of  20  Amps. 

a.  Detailed  Model 

The  SSCM  topology  is  a  standard  Buck  Converter  as  shown  in  Figure 
II.  10  with  parameters  listed  in  Table  II. 6.  The  signal  d  for  the  detailed  model  controls 
the  semiconductor  switch  position  to  achieve  the  desired  command  voltage. 
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Figure  11.10:  Ship's  Service  Converter  Module  Circuit 


Parameter 

Description 

Value 

Units 

h 

Input  filter  inductor 

357.0 

/u\\ 

ra. 

Input  inductor  resistance 

0.2 

Q 

CA 

Input  filter  capacitor  1 

500.0 

//F 

Rjc 

Input  filter  series  resistor 

1.0 

Q 

CL 

Input  filter  capacitor  2 

45.0 

//F 

cout 

Buck  Converter  output  capacitor 

500.0 

M  F 

^ out 

Buck  Converter  dc  inductor 

3.0 

mR 

Rl 

^ out 

Resistance  of  the  dc  inductor 

0.5 

Q 

P 

- 1  out 

Rated  power  output 

8.0 

kW 

Table  II.6:  SSCM  Parameters 
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b.  Average  Value  Model 

Similar  to  the  power  supply  model,  the  SSCM  average  value  model 
simulates  a  fast-average  of  the  semiconductor  switching  dynamics  of  the  detailed  model 
such  that  state  variables  are  constant  in  steady  state.  Again,  the  average  value  model 
maintains  transient  response  accuracy  and  implements  the  same  model  and  control 
parameters.  The  model  topology  depicted  in  Figure  II.  11  shows  the  averaged  command 
signal  (duty  cycle)  d  being  multiplied  by  states  of  the  system,  which  acts  as  a  dependent 
voltage  and  current  source. 


w 


1 


Load 


Figure  11.11:  SSCM  Average  Value  Model  Topology 

c.  SSCM  Control 

Because  of  the  fast  semiconductor  switching  action  of  the  SSCM  and 
near-ideal  nature  of  the  control  scheme,  the  desired  output  voltage  can  be  met  even  in  the 
presence  of  input  disturbances.  For  this  DC  ZEDS  model,  the  controller  design  is 
depicted  in  Figure  11.12  with  associated  parameters  listed  in  Table  II. 7. 
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Figure  11.12:  SSCM  Control  Architecture 


Parameter 

Description 

Value 

Units 

vref 

Output  reference  voltage 

400.0 

V 

KP 

PI  controller,  proportional  gain 

0.1 

K, 

PI  controller,  integral  gain 

23.0 

h, 

Current  feedforward  proportional  gain 

0.015 

f 

switching  frequency 

20.0 

kHz 

Table  II.7:  SSCM  Control  Parameters 


Here,  the  output  voltage  is  regulated  through  a  PI  controller  design  (with 
gains  Kp  and  K. ).  For  simultaneous  current  regulation  the  error  signal,  based  on  the 

difference  between  iout  and  iL  and  scaled  by  hi ,  is  placed  on  the  feedforward  path  and 
added  to  the  PI  control  output. 
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d.  Nonlinear  State  Space  Model 

Similar  to  the  power  supply  model,  the  control  signal  d  is  viewed  as 
external  to  the  model  and  dynamically  coupled  to  states  of  the  system  such  that  the  inputs 
can  be  defined  as 


u\ 
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l _ 

u2 

di, 

^ out 

u3 

dvcfi 

_UA_ 

1 

. 

_ 1 

In  the  same  manner  as  the  power  supply  model,  nonlinear  tenns  are 
expressed  as  unknown  inputs  while  the  remaining  linear  dynamics  can  be  viewed  as  2 
subsystems  (pre-filter  and  buck  portions)  coupled  together  as  shown: 
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with  plant  dynamics 
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and  known  and  unknown  input  matrices  associated  with  linear  and  nonlinear  inputs 
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The  states  of  the  SSCM  depicted  in  Figures  II.  10  and  II.  11  are  the  energy  storage 
elements  from  left  to  right  in  the  model  as 


It  should  be  noted  here  that  there  is  a  limit  to  the  direction  of  the  current 
through  the  inductor  path  such  that  iL  is  always  in  the  positive  direction  defined  in 

Figure  II.  10  and  Figure  11.11.  The  model  handles  this  state  restriction  by  setting  a  lower 
bound  on  the  integrator  to  zero.  The  output  equations  are  defined  as 
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Notice  that  the  feed-through  term  in  (0.20),  defined  as  the  D  matrix,  is 
full  rank,  which  implies  that  full  input  estimation  can  be  achieved  from  an  unknown  input 
observer.  The  SSCM  model  will  be  used  as  an  example  in  the  following  chapter  to 
demonstrate  the  unknown  input  observer  methods  with  input  reconstruction.  Simulation 
results  from  the  Zone  1,  Starboard  SSCM  are  shown  in  Figure  11.13  where  the  DC  ZEDS 
baseline  scenario  depicted  in  Table  II. 5  is  used. 
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SSCM  Stbd  Zone  1  Input:  V, 


Figure  11.13:  Zone  1  Starboard  SSCM  Simulation  Regulates  400  V  Output  (bottom) 

from  500  V  Input  (top) 


4.  Oring  Function 

The  Oring  function  is  the  model  in  Figure  II.  1  that  determines  which  bus  supplies 
the  power  to  the  load  in  each  zone.  The  port  and  starboard  SSCMs  have  slightly  different 
voltage  regulation  set  points  so  that  the  SSCM  output  voltage  with  the  higher  set  point 
picks  up  all  of  the  required  power  to  the  load  in  the  zone.  The  set  points  are  close  so  that 
if  the  primary  SSCM  is  lost,  then  the  alternate  SSCM  picks  up  the  load  requirement  at 
nearly  the  correct  voltage.  As  an  example,  for  the  nominal  simulation  configuration 
shown  in  Table  II. 5,  Zone  1  load  current  is  supplied  from  the  starboard  bus  as  shown  in 
Figure  11.14. 
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Or  Zone  1  Output:  i. 


Load.s 


Time  (sec) 

Figure  11.14:  Oring  Function  Simulation  Results 


The  function  is  modeled  as  two  voltage  sources  on  opposite  ends  of  a  circuit 
providing  current  to  a  common  constant  power  load  through  opposing  diodes  as  shown  in 
Figure  11.15  and  parameters  detailed  in  Table  II. 8.  If  the  inputs  and  states  of  the  state 
space  model  are  defined  as  (respectively) 
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then  the  state  space  model  is  defined  as 
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with  output  equations 
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The  opposing  diodes  in  the  model  are  simulated  as  lower  bounds  on  the  state  integrators. 
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Figure  11.15:  Oring  Function  Circuit  Model 
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Parameter 

Description 

Value 

Units 

K 

Port/stbd  line  resistance 

1.0 

mQ. 

Kr 

Port/stbd  line  inductance 

100.0 

//H 

Table  II.8:  Oring  Function  Parameters 


5.  Inverter  Module  and  Constant  Power  Load 

According  to  the  DC  ZEDS  model,  the  load  in  each  zone  is  a  Ship’s  Service 
Inverter  Module  -  converting  dc  voltage  back  to  ac  -  followed  by  a  constant  power  load 
bank.  Although  the  simulation  model  is  more  detailed,  the  relevant  dynamics  may  be 
represented  by  a  capacitor  and  effective  series  resistance  in  parallel  with  an  ideal  constant 
power  load  of  5.65  kW.  Herein,  the  FDI  analysis  is  conducted  on  the  DC  side  of  the 
power  distribution  system,  which  encompasses  the  rectifier  dc  input  voltage  of  the  Power 
Supply  to  the  dc  voltage  and  current  leading  to  the  SSIM  in  each  zone. 


C.  METHOD  FOR  ANALYZING  DC  ZEDS  SYSTEM 

With  respect  to  the  complete  set  of  state  space  models,  including  power  supply 
and  SSCM  control,  the  DC  ZEDS  model  is  comprised  of  21  distinct  state  space  models 
with  a  total  of  83  inputs,  159  outputs  and  68  states.  Among  the  inputs,  2  are  external  plus 
8  reference  voltage  set  points.  The  remaining  88%  of  inputs  are  either  output 
measurements  from  connecting  modules  or  nonlinear  components  that  are  a  function  of 
control  signals  and  internal  states.  This  fact  demonstrates  the  high  degree  of  dynamic 
interdependence  between  the  modules  and  also  the  model’s  propensity  to  propagate  faults 
from  one  module  to  another.  Among  all  the  output  measurements,  only  about  30%  of 
them  are  directly  used  as  inputs  to  controllers  or  other  DC  ZEDS  modules.  The 
positioning  of  the  remaining  70%  of  output  measurements  is  based  solely  on  enabling 
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observer  sufficiency  conditions  for  full-state  and  input  estimation.  In  other  words,  almost 
70%  of  the  sensors  throughout  the  DC  ZEDS  system  are  used  exclusively  for  detecting 
faults  and  analytical  redundancy. 

It  is  shown  that  the  average  value  model  accurately  represents  the  signature  of  the 
detailed  model  -  even  in  the  transient  response.  In  steady  state,  a  frequency  analysis 
between  the  two  models  reveals  additional  frequency  content  in  the  detailed  model  due  to 
the  switching  action  and  it  can  be  made  to  closely  resemble  the  AVM  dynamics  with  a 
simple  lowpass  fdter  (cutoff  frequency  of  100  Hz). 
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III.  A  POLYNOMIAL  APPROACH  TO  DYNAMIC  MODELING 

AND  FAULT  DETECTION 


The  fundamental  framework  to  detecting  fault  characteristics  in  this  work  is  based 
on  analytical  redundancy.  Indeed,  the  large-scale  dynamically  interdependent  DC  ZEDS 
system  must  rely  on  mathematical  models  to  provide  fault  indications  for  two  main 
reasons:  first,  as  outlined  in  Chapter  I,  it  is  efficient  in  terms  of  hardware  invested 
compared  to  current  methods  [390]  [27 1  ]  [272]  [273]  [289][343]  [90] .  The  second  reason  is 
that  the  DC  ZEDS  model  is  well-defined  as  shown  in  Chapter  II,  which  lends  itself  well 
to  model-based  estimation  techniques[220][335].  By  combining  the  physics-based 
dynamical  models  and  sensor  infonnation,  we  compute  estimates  of  the  states  of  the  DC 
ZEDS  modules.  Then,  a  comparison  of  the  sensor  measurements  and  their  estimates  is 
made,  yielding  residuals  that  are  used  as  a  measure  of  input-output  consistency. 

Perhaps  the  most  widely  used  approach  to  generate  residuals  is  to  measure  the 
consistency  between  observations  and  their  prediction  based  on  a  mathematical  model 
[60].  Standard  approaches  are  the  observer  for  deterministic  linear  and  nonlinear  models 
and  the  Kalman  Filter  for  stochastic  models.  In  both  cases,  the  model-based  estimation 
technique  require  knowledge  of  all  inputs  u(t)  and  outputs  y(t) .  There  are  cases 
however  (especially  in  residual  generation),  where  it  is  desirable  to  estimate  the  system 
when  one  or  more  (perhaps  even  all)  inputs  are  unknown.  For  example,  it  could  be  that 
the  dynamical  model  of  a  system  can  be  separated  into  linear  and  nonlinear  parts  and  the 
nonlinear  part  could  be  modeled  as  an  unknown  input.  Also,  it  might  be  desirable  to 
remove  the  dynamical  influence  of  a  particular  input  from  an  observer  estimate  for  the 
purpose  of  detecting  and  isolating  actuator  faults.  In  either  case,  it  has  been  shown  that 
unknown  input  observers  are  well-suited  to  solve  the  problem  [1 10]  [  1 75][297] . 

One  such  method  to  obtain  an  unknown  input  observer  is  the  subspace  method 
where  the  full  state  is  estimated  based  on  dynamics  residing  in  a  subspace  of  known 
inputs  and  outputs  as  stated  in  [170]  and  [175].  The  conditions  for  UIO  existence  are 
proved  in  [60]  and  [230]  in  which  a  full-state  unknown  input  observer  is  assumed. 

To  this  goal,  consider  a  dynamic  system 
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x{t)  =  Ox + r>j  ( t ) + r2u2  ( t ) 

v(0  =  Cx(t) 


and  we  want  to  find  conditions  to  estimate  the  state  x  without  the  input  u2(t) .  This 
problem  can  be  solved  by  decomposing  the  state  in  two  components  [175] 

x  =  (/-MC)x  +  MCx 
=  q  +  My. 


The  matrix  M  is  chosen  so  that  (/-MC)  is  a  projection  matrix  and  q  =  [I -MC)x  is 
independent  on  the  input  u2(t)  .  The  component  q(t)  satisfies  the  differential  equation 

q(t)  =  (l  ~  MC)  Ox(t)  +  (/  -  MC)T1m1  (t)  +  (l-  MC)  T2u2  (t)  (0.23) 

where  u,  (t)  and  u2(t)  represents  all  the  known  and  unknown  inputs,  respectively. 

If  M  is  such  that 


(I-MC)T2  =0 


(0.24) 


q(t)  is  clearly  independent  of  u2(t)  .  For  a  more  detailed  explanation  of  the  subspace 
method,  refer  to  Appendix  C. 

Lemma  4.1  Equation  (0.24)  is  solvable  iff: 
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rank(CT2 )  =  rank  ( f2 ) 


(0.25) 


Proof:  Necessity:  When  (0.24)  has  a  solution  M  ,  it  becomes  MCY  2  =  T2  or 

(cr,)T  ht  =r2r 

which  means  1%  belongs  to  the  range  space  of  the  matrix  (CT2  )r  resulting  in: 

rank  (r/ )  <  rank |(CT2 j 
and 

rank  (T2 )  <  rank  (Cr2 ) . 

However, 

rank  ( CT2 )  <  min  | rank  ( C) ,  ran  *(r2)}<  rank (r) 

Hence,  rank  (T2 )  =  rank  (Cr2 )  and  the  necessary  condition  is  proved. 

Sufficiency:  When  rank [CT 2)  =  rank (T2)  holds  true,  CT2  is  a  full  column  rank 
matrix  (because  T2  is  assumed  to  be  full  column  rank),  and  a  left  inverse  of  Cr2  exists: 

(CTif  =((cr!)rcr2)  (cr2)r 

Clearly,  M  =  Y2  (Cr2  )  is  a  solution  to  (0.24)  [60],  The  rank  condition  also  implies  that 

the  number  of  independent  output  measurements  must  be  greater  than  the  number  of 
unknown  inputs  [175]. 

A  full-state  UIO  estimate  is  achievable  if  the  rank  condition  (0.25)  is  satisfied. 
This  condition  can  be  rather  restrictive.  A  simple  example  of  a  system  in  controllable 
canonical  fonn  would  show  that  a  necessary  condition  would  be  that  at  least  one  of  the 
transfer  functions  to  have  a  relative  degree  at  most  equal  to  one. 
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For  the  purpose  of  fault  diagnosis,  what  is  actually  needed  is  a  residual  generator, 
not  necessarily  the  full-state  estimate.  In  the  case  where  the  full-state  estimate  is  not 
achievable,  it  is  still  possible  to  compute  the  residual  based  on  the  partial  state.  The 
purpose  of  the  polynomial  method  is  to  take  advantage  of  this  property  by  generating 
residuals  based  on  the  partial  state. 

The  remainder  of  the  chapter  develops  two  unknown  input  observer  (UIO) 
methods  generalized  to  the  multivariable  case  on  a  polynomial  approach.  The  first 
process,  called  the  input  replacement  method,  is  limited  to  a  more  restrictive  class  of 
systems,  but  has  a  simpler  combined  input  and  state  Kalman  estimation  implementation 
algorithm  suitable  for  fault  diagnosis.  By  comparison,  the  second  method  (called  the 
parametrization  method),  takes  a  more  general  class  of  systems  and  is  based  on  the 
nonuniqueness  property  of  the  Diophantine  equation  to  establish  an  unknown  input 
Kalman  estimation  implementation.  In  both  cases,  a  partial-state  representation  of  the 
system  is  sought,  which  is  directly  related  to  the  controllable  canonical  form. 


A.  DYNAMIC  SYSTEMS  AND  THE  PARTIAL  STATE 
1.  Single  Input  Single  Output  (SISO)  Systems 

Consider  a  single  input  single  output  (SISO)  system  represented  as  a  proper 
rational  fraction  of  scalar  polynomials  B(p )  and  A( p)  as  shown  in  Figure  III.  1  and 
expressed  as  a  transfer  function  along  with  its  relationship  to  the  state  space  realization  in 
(0.26). 


;/(/) 

B(p) 

y(0 

A(p) 

Figure  III.l:  SISO  System  in  Transfer  Function  Form 
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H(p)  =  ^^-  =  C(pI-0)1r  +  D 

A(p )  ^  1 


(0.26) 


The  numerator  B(  p)  is  of  order  ns  or  at  least  of  order  ns  -1  (b0  =  0)  for  a  proper 
or  strictly  proper  transfer  function  H(p) ,  respectively,  and  is  of  the  form 

B(p)  =  b0pn‘  +  blp” '  +b2p"s~2  +-  +  b„s .  (0.27) 

Likewise,  let  the  denominator  characteristic  equation  be  a  monic  polynomial  of 
order  ns as 


A(p )  =  p"  +a]p'h  1  +...  +  alh  (0.28) 

where  p  is  the  differential  operator  s  in  continuous  time  and  the  time-shift  operator  z 
in  discrete  time.  If  there  are  no  common  factors  (i.e.  no  pole-zero  cancellations)  between 
B( p)  and  A(p) ,  then  the  transfer  function  is  said  to  be  coprime  and  the  system  is  a 
minimal  realization.  Now  referring  to  Figure  III. 2,  define  z(t )  as  the  partial  state  [200] 
(or  sometimes  called  the  pseudo  state  [61]). 


Figure  III.2:  Partial-state  Space  Representation  of  SISO  System 
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It  is  easy  to  see  that  it  can  be  related  to  the  input  and  output  signal,  u(t)  and  y(t) 


as 


A(p)z(t)  =  u(t)  (0.29) 

y(t)  =  B(p)z(t)  (0.30) 


The  meaning  of  the  partial  state  is  that  (0.26)  can  be  related  directly  to  the  state  of  a 
controllable  canonical  form  realization.  Again,  in  the  SISO  case  this  is  straightforward, 
since  the  state  of  the  system  is 


z(t)  =  [z(t),pz(t),...,pn°  z(t)] 


(0.31) 


where  z(t)  is  a  scalar  and  z(t)  e  M"s . 


2.  Multiple  Input  Multiple  Output  (MIMO)  Systems 

Expanding  the  discussion  of  partial  state  to  the  more  general  multivariable  case, 
consider  the  class  of  systems  represented  in  state  space  as 


px(t)  =  Ox(t)  +  ri<(t)  (0.32) 

y(t)  =  Cx(t )  +  Du(t)  (0.33) 


where  xeR"',nel"‘,vel"v  and  O, T,C,  and D  are  of  appropriate  dimensions. 
Again  p  in  the  state  dynamics  equation  (0.32)  denotes  the  differential  operator  in 
continuous  time  or  the  time-shift  operator  in  discrete  time  and  its  notation  is  consistent 
throughout.  Equations  (0.32)-(0.33)  can  also  be  represented  in  transfer  function  fonn, 


48 


only  the  “numerator”  B(p)  and  “denominator”  A(p)  are  matrices  whose  elements  are 
polynomials.  Together,  //( p)  is  called  the  transfer  function  matrix  or  simply  the 
transfer  matrix  and  can  be  written  as 


H(p)  -  B(p)A  '(p)  (0.34) 

H(p)  =  A-\p)B(p )  (0.35) 


The  right  fraction  pair  B(p )  and  A(p)  has  dimensions  ny  x  nu  and  nu  x  nu ,  respectively 
while  the  left  fraction  pair  A(p)  and  B(p)  has  dimensions  ny  x  n v  and  ny  x  nu .  The 
“denominator”  polynomial  matrices  A(  p)  and  A(p)  are  always  square.  If  the  pair 
(®,T)  is  controllable  and  (®,C)  observable  then  the  system  is  a  minimal  realization, 
and  the  polynomial  pairs  ( B(p),A(p ))  and  ( B(p),  A(p)  J  are  mutually  coprime.  In  the 

SISO  case,  this  implies  that  the  numerator  and  denominator  polynomials  have  no 
common  roots  (no  pole-zero  cancellations).  The  more  general  MIMO  case  requires  an 
appropriate  extension  of  this  concept,  as  shown  in  Appendix  D.  Under  these  conditions, 
the  three  representations 


(®,r,c,£>) 


in  state  space  form, 


(A(p),B(p)) 

in  right  polynomial  matrix  fraction  form,  and 
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( A(p),B(p )) 


in  left  polynomial  matrix  fraction  form  are  equivalent  and  one  can  be  computed  from  the 
other  [61].  Of  particular  interest  is  the  right  coprime  factorization  in  (0.34)  which  leads 
to  the  MIMO  partial-state  representation  of  the  system  shown  in  Figure  III. 3. 


Figure  III.3:  Right  Fraction  Representation  of  MIMO  System 


By  extension  to  the  SISO  systems,  the  partial  state  in  the  multivariable  case  is  related  to 
the  input  and  output  signals  u{t)  and  y(t)  of  (0.29)  and  (0.30)  and  the  system  can  be 
realized  as 

PM'Z\(t) 

\  =-A~h1Alz(t)  +  A-hlu(t)  (0.36) 

_pMn“z„S  t)_ 

y(t )  =  (B,  -  BhAh'A, )  z(t)  +  BhA~M)  (0.37) 

where  the  state  space  equations  (0.36)-(0.37)  are  related  to  the  original  system  (0.32)- 
(0.33)  by  a  nonsingular  matrix  T  such  that 

x(t)  =  Tz(t).  (0.38) 
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The  matrices  Ah ,  Bh  are  matrices  of  high  degree  coefficients  whose  dimensions  are 
nu  x  nu  and  n  x  nu ,  respectively,  and  Al ,  B )  are  the  matrices  of  all  lower  order 
coefficients  with  dimension  nu  x  ns  and  n  x  ns ,  respectively.  It  should  also  be  noted  here 

that  a  method  for  determining  T  in  (0.38)  to  transform  the  original  system  in  (0.32)- 
(0.33)  directly  to  control  canonical  form  in  (0.36)-(0.37)  is  found  in  [386].  A  more 
detailed  explanation  of  the  system  partial-state  realization  for  the  MIMO  case  is  found  in 
[61]  and  in  Appendix  D. 

Similar  to  the  vector  of  partial  states  (with  derivatives)  in  (0.31)  for  the  SISO 
case,  we  see  that  the  partial  state  is  a  vector  related  to  the  system  as 


z(t)=  Zj(0, ,  p^z^t):  z2(t), ... ,  p^~lz2{t)\  ...:zn  (t),...,pM"‘rlz  (t)  (0.39) 


where  z(f)eS"‘  M+ft +•••+«,„  ancj  ^ju2,...,jun  are  called  the  column  degrees  (i.e.  the 
maximum  degree  of  the  polynomials  in  each  column)  of  the  denominator  matrix  A(p) 
and  is  related  to  the  controllability  index  in  [200]  and  [386].  From  (0.39),  it  can  be  seen 
that  the  maximum  degree  of  the  polynomials  in  each  column  of  the  denominator 
polynomial  matrix  A ( p )  is  directly  related  to  the  full-state  estimate  of  the  observer.  It  is 
an  important  distinction,  then,  to  define  the  partial  state  alone  (without  derivatives)  as 


r  ~\T 

z(t)  =  \_zl(t),z2(t),---,z„u(t) 


(0.40) 


where  z(t)  e  M''n  has  exactly  the  same  elements  in  z_(t) ,  only  without  the  associated 
derivatives.  This  reduced  partial-state  vector  is  important  as  it  relates  to  developing 
unknown  input  observers  when  the  UIO  rank  condition  is  not  satisfied. 
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It  should  be  noted  here  that  in  Fault  Diagnosis,  the  goal  is  to  generate  residuals, 
not  necessarily  to  estimate  the  full  state.  We  will  see  that  in  many  situations,  it  is 
possible  to  estimate  the  partial  state  z(t )  but  not  the  entire  state  x(t)  =  Tz(t) ,  especially 
when  u(t)  is  unknown  or  assumed  unknown.  In  this  case,  however,  if  a  full-state 
estimate  must  be  made,  it  still  can  by  including  the  appropriate  number  of  derivatives  to 
the  partial-state  vector  through  a  variety  of  filtering  methods  to  obtain  z(7)  in  (0.39). 


B.  POLYNOMIAL-BASED  UIO  BY  INPUT  REPLACEMENT 
1.  Theory 

In  order  to  illustrate  the  UIO  concept,  let  us  consider  a  simple  case  of  a  SIMO 
system  with  equations  (0.32)-(0.33)  with  one  input  and  two  outputs  as  shown  in  Figure 
III.4.  The  goal  is  to  estimate  the  partial  state  z(t)  based  on  outputs  v,  (/)  and  y2(t)  only. 


A  (f) 

>i(0 

«(0 

A(p) 

B2(p) 

y2  (0 

A(p) 

Figure  III.4:  SIMO  System  Transfer  Function 


As  previously  discussed,  z(t )  satisfies  the  equations 

A(p)z(t )  =  u(t) 

>i(0  =  Bi(p)z(t) 

y2(t)  =  B2(p)z(t). 
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If  the  input  u(t)  is  not  available,  we  can  use  the  last  two  equations  as 


B\(p)z(t)  =  y:  (t) 
y2(t)  =  B2(p)z(t ) 


with  the  same  partial  state  z(t ) . 

From  the  definition  of  the  partial-state  vector  in  (0.31),  a  set  of  equations  can  be 
derived  from  the  input/output  relationships  as  if  v,  (/)  and  y2(t)  are  the  input  and  output 
of  the  system  as  shown  in  Figure  III.  5. 


Figure  III.5:  SIMO  Unknown  Input  System 

Figure  III. 5  shows  that  z(t )  is  the  partial  state  associated  to  the  transfer  function 
B2(p)/Bl(p)  and  therefore  it  can  be  estimated  by  a  standard  observer  or  Kalman  Filter. 

In  order  for  the  UIO  to  exist,  B{( p)  and  B2(p)  must  be  mutually  coprime  (no 

pole-zero  cancellations),  which  is  the  same  as  saying  the  system  in  Figure  III. 5  must  be 
observable.  Note  that  there  are  more  linearly  independent  outputs  (2)  than  there  are 
unknown  inputs  (1),  consistent  with  the  rank  condition  (0.25)  [230].  To  add  robustness 
to  the  state  estimation,  a  standard  Kalman  Filter  can  be  used  to  minimize  the  mean  square 
error  of  the  estimates. 
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In  general,  consider  a  MIMO  system  of  the  fonn  (0.32)-(0.33)  where  at  least  one 
of  nu  inputs  is  unknown  and  the  input  matrix  T  with  corresponding  input  vector  function 
u  are  partitioned  into  known  and  unknown  inputs  as 


Tu{t)=  Tj 


(0.41) 


where  input  is  partitioned  into  n  known  inputs  and  nu  =nu  -n  unknown  inputs.  In 
this  manner,  using  equation  (0.41)  we  can  write 


x(t)  =  Ox(t)  +  r|M|  ( t )  +  Y2u2  ( t )  (0.42) 

y(t)  =  Cx(t)  +  Du(t)  (0.43) 

where  Tj  e  M”"1  and  V 2  e  R”"2  .  For  the  unknown  input  observer  to  exist  for  full-state 
estimation,  rank (Cr2 )  =  rank (T,) ,  which  implies  that  there  must  be  at  least  as  many 

independent  outputs  as  there  are  unknown  inputs.  From  the  right  polynomial  matrix 
fraction  of  the  transfer  function  matrix  //(/;) ,  combine  the  partial-state  input/output 
equations  (0.29)  and  (0.30)  as 


A(p) 

u(t) 

z(t)  = 

Bip)  \ 

(0.44) 


where  A(p)  e  R""x"“ ,  B(p)  e  R"2’*"" ,  z(t)<=  En" ,  u(t)  e  R"“ ,  and  It  can  be 

seen  in  equation  (0.44)  that  there  are  a  total  of  nu  +ny  equations  relating  the  input  and 
output  measurements  to  the  partial  state  z(t ) .  The  main  idea  is  to  replace  all  the 


54 


unknown  input  row  equations  in  A(p)  with  the  same  number  of  output  row  equations  in 
B(p )  to  obtain  a  new  transfer  matrix 


H(p)  =  B(p)A\p) 


(0.45) 


where  "''2X,'"  and  A(p)  e  .  Right  coprimeness  of  flip)  and  column 

reducedness  of  A(p )  must  be  verified  since  there  is  no  guarantee  that  the  system  will 
have  satisfied  those  conditions. 

To  illustrate  this  argument,  suppose  a  system  described  in  (0.42)-(0.43)  where 
ns=  5,  nu=  2,  and  n  =4  and  suppose  the  second  input  u2(t)  is  unknown.  Then 

equation  (0.44)  is  represented  as 


1^10 P  1^11  P  1  ®nP  +  1^13  •  1^20 P  1^21  P  1  d 22 

«1  (0 

2®10 P  2®\lP  2^12 P  2^13  ‘  2®2oP  2®2lP  2®22 

u2(t) 

(0 

A0p3  +  AiP2  +  A2P1  +  Ai  '■  A0P2  +  A21P1  +  An 

= 

yA) 

2b\0p3  +  2bnP2  +  2bnPl  +  2bu  :  2b20p2  +  2b2lpl  +  2b22 

z2(t)_ 

y2(!) 

AioP3  +  AlP2  +  AiP*  +  3^13  :  A2oP2  +  A21P  +  3b22 

yA) 

_  AoP3  +  AlP2  +  AnP  +  4bu  :  A20P2  +  A21P1  +  A22  _ 

yA)_ 

It  is  desired  to  have  the  third  output  y3  (t)  replace  the  unknown  input  u2  (t)  in  the  input 
polynomial  matrix  A( p) .  The  rank  condition  has  been  satisfied  and  assuming  the 
leading  coefficients  on  the  polynomials  for  the  third  output,  3h10  and  ,b2(t  are  nonzero 
(which  implies  that  Ah  in  (0.36)-(0.37)  is  nonsingular),  the  replacement  can  be  made  as 
shown 
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1^10 P  \®\\P  \®YlP  +  1^13 
3^10  P  +  3^11  P  +  3^12  P  +  3^13 


A0p3  +  AnP2  +  AiA  +  An  ■ 
iKp3  +  AnP2  +  2bnPl  +  An  '■ 

AoP3  +  AnP2  +  AnP  +  An  '■ 


1^20  P  1^21  P  1^22 

«i(0 

AioP  +  Ai\P  "*■  Aii 

zA) 

t3(0 

i  Kp2  +  AnP  +  An 

= 

Ti(0 

AioP  AnP  An 

M0. 

y2(0 

AioP2  +  AnP  +  All. 

M0_ 

The  top  nu  =  2  rows  are  now  the  new  input  polynomial  matrix  A( p)  and  the  bottom 
«v-l  =  3  rows  comprise  the  new  B(p) .  The  transfer  matrix  H(p)  =  B(p)Al(p)  is 

checked  for  right  coprimeness  and  A(p)  is  verified  to  be  column  reduced.  Finally,  the 
new  system  of  equations  is  represented  in  state  space  form  for  Kalman  Filtering 
implementation.  Note  that  the  estimate  for  the  new  system  will  be  on  the  partial  state,  so 
it  must  be  transformed  through  the  linear  transformation  in  (0.38)  if  the  original  full-state 
estimation  is  desired. 


The  relationship  of  the  UIO  rank  condition  in  (0.25)  to  the  polynomial 
representation  for  full-state  estimation  implies  that  B(p)  must  have  higher  degree 
coefficients.  Consider  the  input  and  output  matrices  T,C  for  the  controllable  canonical 
form  realization  both  partitioned  as 


(0.46) 


with 


rf  =  [0  •  •  •  0  y,  ]  e  RVMt 

c.  =  [x  •••  x  c;]  e  R"yXM‘ 


(0.47) 
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for  i  =  \,...,nu.  It  turns  out  that  the  terms  Ct  are  the  coefficients  of  the  polynomial 
B(p)  .  In  particular  we  can  write 


Ap)  = 
B(p)  = 


Y\ 


Yn. 


£  Ap)+  tenns  of  lower  degree 


(0.48) 


PP  •••  c.p 


ft, -1 


+  tenns  of  lower  degree 


where 


^M(p)  =  diag  j/U1 , p**1 J  e  IT 


In  the  case  of  a  UIO  we  can  see  some  of  the  conditions  of  existence.  Let  r  be 
partitioned  as  T  =  ^Tl  T2  J  where  T2  is  associated  with  the  unknown  inputs.  In  Kudva 
et  al  [230],  the  condition  for  existence  of  the  UIO  for  full-state  estimation  is 
rank(cr2)  =  rank(T2) .  From  (0.46)-(0.48)  it  can  be  seen  that 

(0.49) 

i= 1 

with  y.  the  part  of  the  vector  y.  associated  to  the  unknown  input.  From  (0.48)  we  see 
that  a  necessary  condition  is  that  not  all  ci  can  be  zero,  which  means  that  the  polynomial 
B(p )  must  have  higher  degree  terms. 

While  the  unknown  input  observer  rank  condition  (0.25)  must  be  satisfied  for  full- 
state  estimation,  it  does  not  need  to  be  satisfied  for  estimating  the  partial  state.  This  is  the 
advantage  of  using  the  polynomial  UIO  technique  and  is  illustrated  in  the  following 
example. 
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2.  Example 

Consider  a  third-order,  single-input,  two-output  system  shown  in  Figure  III.4 
where  the  continuous  time  system  is  represented  in  right  matrix  fraction  form  as  in  (0.34) 

with  B(s)  =  [5  +  1  s  +  2]r  and  A(s)  =  s3  +6s2  +  1N  +  6  .  It  can  be  seen  that  the 
greatest  common  right  divisor  between  B(s)  and  A(s)  is  1,  which  is  unimodular.  In  [61] 
this  implies  that  the  system  transfer  matrix  H(p )  is  controllable  and  observable  and  an 
equivalent  minimal  realization  representation  of  the  system  in  state  space  form  is 


~  0 

1 

O' 

~0~ 

x(t)  = 

0 

0 

1 

x(t)  + 

0 

-6 

-11 

-6 

1 

u(t) 


y(t)  = 


1  1  0 

2  1  0 


x(t)  + 


u(t) 


(0.50) 


(0.51) 


Notice  all  the  terms  in  (0.50)-(0.51)  are  directly  related  to  the  coefficients  of  the 
polynomials  A(s),  Bt(s),  and  B2(s )  .  By  observation,  the  system  of  equations  is  indeed 

controllable  because  it  is  realized  in  controllable  canonical  form  where  the  state  vector 
x(t)  happens  to  be  the  partial-state  vector  z(/)  with  derivatives  up  to  order  2  as  shown  in 
(0.31).  Since  there  is  only  one  input,  there  is  only  one  partial  state  z(t).  It  can  be 
verified  that  the  system  of  equations  is  also  observable.  To  design  a  full-state  unknown 
input  observer,  the  rank  condition  (0.25)  must  be  satisfied.  For  this  example,  it  means 
that 


6 

"O' 

5 

6 

"O' 

X 

'1 

1  0" 

rank 

2 

1  °_ 

0 

rank 

0 

V 

1 

J 

1 

J 
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and  a  full-state  UIO  does  not  exist  for  this  system.  However,  an  unknown  input 
observer  still  exists  from  a  polynomial  approach. 

With  the  polynomial  input  replacement  UIO  method  described  in  this  chapter,  an 
estimate  of  the  partial  state  can  still  be  achieved  as  outlined  in  the  following  steps  to  the 
algorithm.  We  start  with  the  example  state  space  representation  in  (0.50)-(0.5 1). 


1)  Represent  the  system  of  equations  in  right  matrix  fraction  from  as  in  (0.34)  and 
ensure  that  H(p )  is  right  coprime,  which  implies  a  minimum  realization. 


H(s)  =  B(s)A~\s)  = 


5  +  1 
s  +  2 


1 


s: +6s  +  II5  +  6 


The  polynomial  transfer  matrix  is  right  coprime. 


2)  Represent  the  system  of  input/output  equations  in  partial-state  form  as  in  (0.44). 


s3  +6s2  +  II5  +  6 

u(t) 

5  +  1 

z(t)  = 

Ti(0 

5  +  2 

_y2(  0_ 

3)  Determine  the  right  matrix  fraction  form  of  the  UIO  where  //(/?)  =  B(p)  A  '( p) 
by  replacing  the  rows  in  the  input  polynomial  matrix  A(p )  that  correspond  to 
u2(t )  with  the  rows  in  the  output  polynomial  matrix  B(p)  that  correspond  to  the 
desired  input  replacement  equations. 
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For  this  example,  there  is  only  one  input  and  two  options  for  replacement.  Let  us 
choose  the  output  yft)  to  replace  u(t) . 


s  +  \ 

Ti(0 

z(t)  = 

s  +  2 

_t2(0_ 

This  leads  to  the  system  shown  in  Figure  III.5  with  the  new  scalar  transfer 
function  H(p )  written  as 


H(p)  = 


Bis) 

A(s) 


B2(s ) 
Bt  (s) 


s  +  2 

J+T' 


Since  H(p )  is  a  first  order  system  -  while  the  original  //(/;)  was  third  order  - 
we  can  estimate  z(t )  ,  but  not  z{t)  and  z(t)  .  This  means  that  we  can  still  create  a 
residual,  but  not  a  full-state  estimate.  The  principle  can  be  generalized  to  say  that 
for  any  column  degree  reduction  of  the  original  A(s )  to  the  new  A(s),  it  will 
result  in  an  unknown  input  observer  that  cannot  obtain  a  full-state  estimate. 

4)  Check  H{p )  =  B{p)A  '(/?)  for  right  coprimeness  and  column  reducedness. 


The  new  transfer  function  is 

H(s)  =  B(s)/A(s)  =  ^ 
s  +  l 
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since  B(s)  and  A(s )  are  scalar  polynomials.  By  inspection,  there  are  no  pole- 
zero  cancellations  of  H(s)  and  a  scalar  polynomial  is  always  column  (and  row) 
reduced  [61].  Note  that  an  equivalent  observer  (in  the  steady  state)  can  be 
determined  from  II  '(.s  ) . 


5)  Represent  H{p )  in  state  space  form  and  implement  any  standard  observer 
method.  If  the  UIO  can  estimate  the  partial  state  (i.e.  the  UIO  rank  condition  is 
satisfied),  determine  the  nonsingular  matrix  T  such  that  x(t)  =  Tzft )  to  determine 
the  original  state  estimate. 


For  this  example,  the  UIO  rank  condition  is  not  satisfied  and  a  full-state  estimate 
cannot  be  achieved.  The  state  space  representation  for  H(s)  is 


z(t)  =  -z  +  y]  (0  +  v, 

y2(t) =  z  +  y,  (t)  +  v2 


where  the  “input”  is  yl (t)  and  the  “output”  is  y2 (/) .  The  additional  terms  vl ,  v2 
account  for  possible  measurement  noise.  For  this  problem,  let  the  input 
w(0  =  sin(0  and  the  measurement  noise  variances  cr2  and  cr2  be  1%  of  the 

amplitude  of  the  input  signal.  Also,  all  measurement  noise  are  assumed  to  be 
uncorrelated.  Under  these  conditions,  the  Kalman  Gain  is  L  =  0.58 11  and  the 
observer  implementation  equation  is 


1(0  =  -1-581  lz(0  +  [0.4189 


0.5811] 


u(0 

v2(0 


and  the  output  is  the  partial-state  estimate  z(t)  itself.  From  the  form  of  the 
observer  equation,  it  can  be  easily  seen  that  the  estimate  is  independent  of  u(t) . 
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A  simulation  of  the  polynomial  UIO  input  replacement  method  is  shown 
in  Figure  III. 6  and  the  noise  terms  were  removed  in  the  simulation  to  see  the 
convergence. 


Figure  III. 6:  UIO  Input  Replacement  Method  Simulation 


In  this  example  the  replacement  of  the  unknown  input  u(t)  with  one  of  the 
outputs  yx(t)  is  straightforward,  since  the  resulting  transfer  function  B2(s)/ B] (,s)  is 

scalar.  In  the  more  general  MIMO  case,  the  system  resulting  from  input  substitution  has 
to  be  realizable,  which  might  not  always  be  the  case.  In  the  next  section  we  present  a 
general  observer  parametrization  which  has  the  unknown  input  observer  as  a  particular 
representation. 
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C.  POLYNOMIAL-BASED  UIO  BY  PARAMETRIZATION 

1.  Theory 

In  the  previous  section,  the  UIO  input  replacement  method  took  advantage  of  the 
partial-state  realization  in  right  polynomial  transfer  matrix  form  (0.34).  This  polynomial 
method,  on  the  other  hand,  is  based  on  the  left  polynomial  transfer  matrix  in  (0.35). 

In  this  section  a  polynomial  transfer  matrix  formulation  of  the  observer  or  Kalman 
Filter  yields  a  class  of  stable  solutions  easily  parametrized.  Consider  the  LTI  system  in 
(0.42)-(0.43)  rewritten  here  as 


x(t)  =  Ox(0+rlWl(0+r2M2(0  (0.52) 

y{t )  =  Cx(t)  +  Dxux  (t)  +  D2u2  (t)  (0.53) 


where  and  O,  T,  and  C  are  of  appropriate  dimensions. 

Furthermore,  T,  e  M""1  and  T2  €  M”"2  are  the  subdivided  constant  input  matrices  such  that 
ux{t)  and  u2(t)  are  vectors  containing  all  the  known  and  unknown  inputs,  respectively. 

If  the  system  (0.52)-(0.53)  is  a  minimal  realization,  it  can  be  equally  realized  in  partial- 
state  form  from  (0.44)  with  the  additional  known  and  unknown  input  partitioned  as 


u(t)  = 


«l  (0 

AC p) 

u2(t) 

_A(p)_ 

z(t) 


(0.54) 


y(0  =  B(p)z(0 


(0.55) 


Now  from  the  general  form  for  the  Kalman  estimator,  a  polynomial  representation  of  the 
observer  that  estimates  the  partial  state  z(t )  is  expressed  as 
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C(p)z(t)  =  E(p)u(t )  +  F(p)y(t) 


(0.56) 


with  C(p)  e  ,  E( p)  e  M""xn" ,  and  F( p)  e  M""x"v .  These  polynomial  matrices  can  be 
computed  either  from  the  Diophantine  Equation  (Appendix  D) 

C(p)  =  E(p)A(p)  +  F(p)B(p)  (0.57) 

or  directly  from  (0.32)-(0.33)  as  an  observer  of  the  fonn 

pz(t)  =  T~l(  ®  -  LC )  Tz(t )  +  T~lru(t)  +  TlLy(t) 
y  if )  =  Cz(t ) 


where  y* (t )  is  the  partial-state  vector  estimate  (with  derivatives)  z(/)  e  K"1  in  (0.39) 
when  the  UIO  rank  condition  is  satisfied  (C  =  In  j,  and  the  partial  state  z(t)  e  M""  in 

(0.40)  when  it  is  not.  In  the  latter  case,  C  is  an  nu  x  ns  matrix  that  contains  the  rows  of 
In  that  extracts  the  partial-state  estimate  z(t)  e  M'!"  from  the  partial-state  vector  (with 

derivatives)  z(t)  e  M"' .  Also,  T  is  the  nonsingular  matrix  transformation  to  the  partial- 
state  space  shown  as  x{t)  =  Tz(t) . 

For  convergence  of  the  observer  in  polynomial  form,  it  is  shown  that  the 
substitution  of  u(t)  and  y(t)  in  (0.54)-(0.55)  into  the  polynomial  representation  of  the 
partial-state  observer  in  (0.56),  along  with  the  Diophantine  equation  identity  in  (0.57), 
yields  the  state  error  equation 


C(p){z(t)~z(t))  =  0. 


(0.59) 
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Equation  (0.59)  implies  that  limz(t)-z(t)  =  0  provided  det(C(/?))  has  all  roots  in  the 
stable  region. 

The  useful  feature  of  the  Diophantine  equation  approach  is  that  the  solution  is  not 
unique.  In  fact,  if  E0(p)  e  and  F0(/?)el"'x"’  together  are  a  solution  of  the 

Diophantine  Equation  (0.57)  it  can  be  shown  in  [61],  [200],  and  [229]  that,  for  any 
arbitrary  polynomial  matrix  W(p )  of  appropriate  dimensions 


E(p)  =  E0(p)-W(p)B(p ) 

(0.60) 

F(P)  =  F0  (p)  +  W  ( p)A(p ) 

(0.61) 

are  solutions  to  the  differential  equations.  As  a  consequence,  from  (0.56)  and  (0.60)- 
(0.61)  we  obtain  a  family  of  observers,  parametrized  by  W(p)  e  M"“x"v  as 


C(p)z(t)  =  {E0(p)-W(p)Bip))u(t)  +  {F0(p)  +  W(p)A(p))y(t)  (0.62) 


where  A(p )  e  M”vX”v  and  B(  p)  e  are  the  left  coprime  representation  of  the  original 
system  transfer  matrix  flip)  =  A  l(p)B(p) . 

Now  partitioning  u(t)  into  known  and  unknown  inputs  we  obtain 


E(p)u(t)  = 


[ E01(P )  E(n(p)]-[WSp)  W2(p )] 


Bn  (p)  Bu(p) 
B2l(p)  B22(p) 


«i  (0 
u2(t) 


(0.63) 
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where  £01(/;)£l"'x\  £02(/))elv'\  ^(^)etvvv,  W2(p)  e  f?”"*""2 , 

2?n(/0  e  M"2'""1'2*”"1  ,  Bn(p)  e  m”v_”"2><h"2  ,  521(^)  e  M""2*""1  ,  and  522(/>)  e  M""2*""2 .  In  order 
to  make  the  observer  insensitive  to  u2(t) ,  we  choose  the  polynomial  matrix  W(p)  so 
that 


£02  (P)  ~  Wx  ( p)Bn  {p)  +  W2(p)B22(p)  =  0 . 


It  can  be  seen  that  the  row  dimension  of  Bv(p)  requires  that  ny  >  n  ,  which  means  that 

in  order  to  have  a  solution,  the  number  of  independent  outputs  must  be  greater  than  the 
number  of  unknown  inputs.  This  is  the  same  result  in  Lemma  4.1  where 
rank  (CT2 )  =  rank  (T, ) . 

If  we  choose  two  polynomial  matrices  Wx(p),  and  W2(p)  that  satisfy  the 
Diophantine  equation 


E02  (p)  =  Wx  (p)Bl2  (p)  +  W2  (p)B22  (p)  (0.64) 


then  a  viable  solution  can  be  obtained,  provided  that  Bl2(p),  and  B 22{p)  are  right 
coprime.  Finally,  combine  (0.62)-(0.64)  to  obtain 


C(p)z(t)  = 


r 

E0l(p)-W\p ) 

V 


BU(P) 

b22(p). 


«i  (0  +  ( E0(p )  +  W*  (p)A(p)'j  y(t) 


(0.65) 


with  W*(p)  =  [Wx(p)  W2(p)\  being  the  combined  solution  to  (0.64).  From  (0.65),  it  can 
be  seen  that  the  parametrized  observer  is  now  independent  of  u2{t) . 
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2.  Example 

In  summary,  the  steps  to  obtain  a  UIO  using  parametrization  is  shown  in  the 
following  example.  The  goal  is  to  demonstrate  how  the  observer  is  obtained  and 
compare  the  simulation  results  to  the  subspace  projection  method  in  [175]  and  to  a 
standard  Kalman  Filter  estimation. 

We  start  with  the  state  space  representation  in  (0.52)-(0.53)  where  the  system, 
introduced  in  [175],  is  a  fifth-order,  two-input,  four-output  model  where 


0.0000 

0.0000 

1.0000 

0.0000 

0.0000 

0.0000 

-0.1540 

-0.0042 

1.5400 

0.0000 

0.0000 

0.2490 

-1.0000 

-5.2000 

0.0000 

0.0386 

-0.9960 

-0.0003 

-0.1170 

0.0000 

0.0000 

0.5000 

0.0000 

0.0000 

-0.5000 

"  0.0000 

0.0000' 

"0 

1 

0 

0 

-f 

-0.7440 

-0.0320 

0 

0 

1 

0 

0 

0.3370 

-1.1200 

,  and  C  = 

0 

0 

0 

1 

0 

0.0200 

0.0000 

1 

0 

0 

0 

0 

0.0000 

0.0000_ 

For  this  problem,  both  inputs  are  considered  unknown  where 
u2(t)  =  [cos(t)  sin(t)]r  .  Added  white  Gaussian  noise  is  also  imposed  on  the  inputs  and 

output  measurements  where  the  variance  for  each  is  1.0%  of  the  amplitude  of  the  input. 
All  measurements  are  assumed  to  be  uncorrelated.  Given  that  the  initial  conditions  in 
[175]  are  arbitrarily  chosen  to  be 


jc(0)  =  [0.3420  0.3200  0.0718  -0.2870  -0.9497]r 
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and  the  initial  conditions  of  the  observers  are  set  to  zero,  we  set  out  to  determine  the 
polynomial-based  UIO  from  the  following  steps: 


1)  Verify  the  system  is  a  minimal  realization  and  transform  the  state  space  system  to 
controllable  canonical  form. 


The  system  is  controllable  and  observable  (and  therefore  a  minimal  realization 
[61]).  From  the  system  of  equations 


fit)  =  r'orz(t) + rf2u2  it) 


y  =  CTfyt) , 


and  the  system  matrices  ® ,  T2 ,  and  C  previously  defined,  the  nonsingular 
matrix  T  is  determined  to  be 


-0.0001 

0.0119 

-0.0003 

0.4478 

0.8941 

0.0173 

-0.0009 

0.0000 

-0.0323 

-0.4470 

-0.0012 

-0.1912 

0.0173 

0.0023 

0.2235 

-0.9950 

0.0024 

-0.0001 

0.0894 

-0.0449 

0.0034 

-0.1119 

-0.9950 

-0.0065 

0.0225 

which  transforms  the  original  system  into  controllable  canonical  form  where  the 
partial-state  vector  (with  derivatives)  is  z(t)  =  Tlxf) . 


2)  Design  a  standard  observer  or  Kalman  Filter  based  on  the  partial-state 
realization. 
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From  the  model,  the  plant  and  measurement  noise  covariance  matrices  are  set  to 
be  Qw  =0.1  In  and  Rv  =  0.01  In  ,  respectively.  The  steady-state  Kalman  Gain  is 

determined  to  be 


0.0225 

-0.0471 

0.0174 

-0.0135 

-0.0115 

-0.0377 

0.0277 

0.0098 

-0.0513 

0.0563 

-0.0015 

0.0267 

0.0937 

-0.5402 

0.0793 

-0.8628 

0.1225 

-1.3764 

0.3712 

-0.4888 

and  the  full-state  Kalman  Filter  can  be  implemented  as  in  equation  (0.58). 


3)  Check  UIO  rank  condition.  Determine  C  such  that 

a.  If  the  rank  condition  is  satisfied,  then  a  full-state  estimate  can  be  achieved 
and  C  =  I n  . 

b.  If  the  rank  condition  is  not  satisfied,  then  determine  C  such  that  nu 
partial  states  are  extracted  (no  derivatives)  from  z(t)  resulting  in  the 
partial  state  z(t )  in  (0.40). 


The  rank  condition  is  satisfied  and  therefore  for  this  example  C  =  I5. 


4)  Determine  the  polynomial  matrices  C(p),  E0(p),  and  F0(p)from  (0.57),  which 

is  the  equivalent  polynomial  representation  of  the  partial-state  observer  designed 
Step  3. 


In  the  example, 
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-s-0.0420  0.9300 


C(s)  = 


-0.0740  -s-0.1500 
-1.0000  -2.2000 
-0.0020  -0.2500 

12.0000  26.0000 


0.0720  -0.0150 

0.9600  0.0073 

-s- 1.1 000  -0.0097 

-0.7900  s+0.8700 
3.8000  -0.0560 


-0.0460 

-0.0390 

0.0560 

-0.4700 

-s-2.3000 


and  E0(s )  and  F0(s)  are  constant  polynomial  matrices 


"  0.0000 

0.0000' 

-0.0091 

0.0470 

-0.0170 

0.0130" 

0.0000 

0.0000 

0.0120 

0.0380 

-0.0280 

-0.0098 

-0.1500 

0.0130 

,  F0(s)  = 

0.0510 

-0.0560 

0.0015 

-0.0270 

0.0000 

0.0000 

0.0940 

-0.5400 

0.0790 

-0.8600 

0.2500 

-1.1000 

-0.1200 

1.4000 

-0.3700 

0.4900 

5)  In  (0.62),  determine  A(p)  and  B(p )  from  the  left  coprime  realization  of  the 
original  system,  H(p)  =  A  l{p)B{p) . 


Because  the  example  is  a  minimal  realization,  the  original  system  of  equations 
defined  by  the  matrices  ®  ,  T2 ,  and  C  can  be  realized  in  left  matrix  fraction  form 
where 


52  +  0.8455  +  0.470 


A(s)  = 


-0.200s-0.098 

0.08s+0.040 


-0.025s-0.012 


■2.400s  -2.5005-0.610 
-0.0535  +  0.035 
0.5705  +  1.300 
-0.2105  +  0.440 


0.4505- -12.0005-5.100 
0.0175  +  0.029 
0.1305  +  2.900 
-0.0495-1.100 


-0.9552 -0.415  +  0.010 
-0.0895-0.001 
-0.6805-0.005 
-0.655  +  0.002 


and 
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0.670.9-0.280  2.700.9  +  1.100 
0.130  0.066 

0.140  -0.640 

-0.054  0.240 

Notice  that  the  sum  of  the  row  degrees  (maximum  polynomial  degree  of  each 
row)  of  A(s )  is  equal  to  the  total  number  of  states  (ns  =5),  which  implies  that 

A(s)  is  row  reduced.  Since  the  degree  of  the  i'h  row  in  B(s)  is  less  than  the 
degree  of  the  i"‘  row  in  A(s)  ,  the  polynomial  transfer  matrix  is  proper  [61]. 

6)  From  the  structure  of  the  input  vector  in  (0.52),  determine  uft )  and  lift) . 
Partition  polynomial  matrices  Efp )  and  B(p)  accordingly  from  equation  (0.63) 

Since  both  inputs  are  unknown,  n  =0  and  n  =2  and  so  E0fs)  is  an  empty 
matrix  and  E0fs)  =  Efs) .  From  the  dimensions  of  the  polynomials  in  equation 
(0.63),  it  can  be  seen  that  from  the  previous  step,  B(s)  is  partitioned  as 

[-0.6705-0.280  2.700.9  +  1.100] 

B.fs)  =  ,  and  B„(s)  = 

12  [_  0.130  0.066 

7)  Verify  Bxfp )  and  B2fp)  are  right  coprime  and  solve  Diophantine  equation  to 
obtain  Wfp)  and  Wfp)  in  (0.64). 

In  the  example,  Bn(s)  and  B0fs)  are  right  coprime.  The  solution  to  the 
Diophantine  equation 


0.140  -0.640' 
-0.054  0.240 
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E01  ( P)  =  WX  (p)Bn  ( p )  +  W2(p)B22  {p) 


are  the  constant  polynomial  matrices 


Wl(s)  = 


"0.000 

0.000" 

"  0.000 

0.000" 

0.000 

0.000 

0.000 

0.000 

0.000 

-1.000 

and  W1  (s)  = 

-0.100 

0.064 

0.000 

0.000 

0.000 

0.000 

0.000 

0.099 

1.500 

-0.580 

8)  Form  W*(p)  =  \Wfp)  W2(p)]  and  from  (0.62)  go  back  to  the  state  space 
realization  of  the  partial-state  representation  for  simulation. 


The  observer  from  the  original  system  has  now  been  successfully  designed  such 
that,  in  choosing  the  polynomial  matrix  W*(p )  in  (0.62),  the  parametrized 
Kalman  Filter  in  polynomial  form  is  now  independent  of  unknown  inputs  u2{t)  . 

Changing  back  to  state  space  form  from  the  left  matrix  fraction  realization  yields 
the  following  system  of  equations: 


-0.0419 

0.9311 

0.0719 

0.0150 

-0.4610 

-0.0743 

-0.1462 

0.9586 

-0.0073 

-0.0386 

-1.0051 

-2.1930 

-1.1025 

0.0097 

0.0556 

-0.0020 

-0.2488 

-0.7881 

-0.8657 

-0.4708 

12.1667 

25.8504 

3.7543 

0.0563 

-2.3256 
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" -0.0008 

0.0000 

-0.0304 

0.0532 

0.1894 

-0.0172 

-0.0609 

0.1289 

r*  = 

-0.0568 

-0.1159 

-0.3352 

-0.1925 

-0.1118 

-0.9950 

-0.0027 

-0.6358 

0.3771 

0.6721 

3.9845 

2.4990 

c*  = 


and 


0.2405 

-0.5811 

-0.0000 

-0.9950 

0.0000 

0.1779 

2.9439 

5.1764 

-0.0449 

0.0899 

0.0000 

-0.2405 

-0.5811 

0.0000 

0.9950 

2.5932 

-5.2463 

-0.1378 

0.0906 

-0.0016 

0.1779 

2.5882 

0.0000 

-0.0449 

-0.0000 

For  implementation,  the  resulting  unknown  input  observer  is  independent  of  u(t) 
and  the  output  is  the  full  partial  state  (with  derivatives)  as 


q(t)  =  ®*q(t)  +  T*y(t) 
z(t)  =  C*q(t)  +  D*y(t ) 


where  q(t)  denotes  the  state  of  the  new  system. 

A  simulation  was  conducted  using  the  above  example  and  the  estimation 
results  on  the  original  state  x4  are  shown  in  Figure  III. 7.  The  red  line  shows  the 
actual  state  trajectory,  while  the  other  trajectories  are  from  the  three  different 
observer  designs.  It  can  be  seen  that  the  polynomial  UIO  method  (blue)  outlined 
in  this  section  performed  just  as  good  as  the  subspace  projection  UIO  method 
(green).  Both  UIO  methods,  which  only  used  output  measurements  to  estimate 
the  state,  perfonned  as  well  as  the  standard  Kalman  Filter,  which  depends  on  both 
inputs  and  outputs  to  estimate  the  state. 
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Figure  III. 7:  State  Estimate  Comparison 

When  deciding  on  a  method  for  estimation,  it  is  desirable  to  design  an  observer  so 
that  one  or  more  inputs  have  no  effect  on  the  state  estimation.  Perhaps  a  system  model 
has  an  unknown  disturbance  or  other  model  uncertainties.  In  the  case  of  residual 
generation,  perhaps  an  input  is  assumed  unknown  to  develop  a  structured  residual  set  for 
Fault  Diagnosis.  When  considering  the  subspace  UIO  methods  and  the  polynomial  ones 
described  herein,  each  is  designed  to  effectively  obtain  an  estimate  of  the  full  state  when 
the  UIO  rank  condition  is  satisfied.  The  advantage  of  the  polynomial  methods  is  that 
they  can  still  generate  a  residual  through  an  estimate  of  the  partial  state. 

As  mentioned  before,  the  input  replacement  UIO  method  is  limited  if  the 
replacement  polynomials  do  not  have  higher  degrees.  As  a  consequence,  this  method 
works  best  when  the  degree  of  the  numerator  polynomial  matrix  B(p )  is  the  same  as  the 
denominator  polynomial  matrix  A(p ) ,  which  is  the  same  as  saying  that  D  ^  0  in  the 
equivalent  state  space  realization. 
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The  parametrization  UIO  method  is  the  more  general  polynomial  method  and  is 
comparable  to  the  subspace  projection  UIO  in  complexity  and  performance.  The  main 
idea  is  to  develop  a  standard  Kalman  Filter  in  polynomial  form  and  parametrize  the 
solution  such  that  it  is  made  independent  of  the  unknown  inputs.  This  method  works  well 
in  the  presence  of  plant  and  sensor  disturbances  with  known  first  and  second  order 
statistics.  One  drawback  common  to  both  polynomial  methods  is  that  they  require  a 
linear  transformation  to  the  controllable  canonical  form  to  accommodate  estimating  the 
partial  state.  In  cases  where  system  eigenvalues  are  very  large  (i.e.  greater  than  1 06 ),  the 
inverse  of  the  transformation  matrix  T  can  be  near  singular. 

Besides  being  able  to  generate  residuals  based  on  the  partial  state,  another 
advantage  to  the  parametrization  method  is  that  since  the  solution  to  the  Diophantine 
equation  is  nonunique,  it  can  be  represented  by  a  family  of  observers.  If  one  of  the  inputs 
is  suddenly  made  unavailable  for  the  estimation,  the  observer  can  be  adjusted  adaptively. 
This  idea  is  being  explored  as  a  developing  possibility. 


D.  INPUT  OBSERVABILITY  AND  ESTIMATION 

Throughout  the  discussion  so  far  in  this  chapter  is  the  idea  that  if  the  UIO  rank 
condition  (0.25)  is  not  satisfied,  a  partial-state  estimate  is  still  possible  and  can  be 
obtained  through  polynomial  UIO  techniques.  Specifically,  it  has  also  been  shown  that 
the  column  degrees  of  the  denominator  polynomial  matrix  A(p)  (in  the  right  matrix 

fraction  realization  H(p)  =  B(p)A  ( p) )  are  directly  related  to  the  ability  to  estimate  the 
full  state  in  (0.39).  For  any  column  degree  reduction  of  the  original  A(p )  to  the  new 
observer  denominator  polynomial  A( p) ,  it  will  result  in  an  unknown  input  observer  that 
cannot  obtain  a  full-state  estimate.  This  also  implies  that  B(p )  must  have  higher  order 
terms. 

We  have  discussed  residual  generation  based  on  the  state  and  partial  state.  Now 
we  consider  it  based  on  input  and  output  measurements.  It  will  be  shown  that  in  some 
cases,  it  is  desirable  to  compute  residuals  not  only  on  the  outputs,  but  on  the  inputs  as 
well.  Depending  on  system  structure,  this  would  give  added  information  on  the  class  of 
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faults.  It  will  be  shown  that  while  output  estimation  is  straightforward  when  the  full-state 
estimate  is  determined,  the  input  estimation  problem  is  non-trivial  and  it  requires  certain 
conditions  for  existence  that  relate  to  the  system  inversion  problem. 

Compared  to  state  estimation  techniques  for  unknown  input  observers,  little 
research  has  been  conducted  on  input  estimation  [73].  Some  input  reconstruction 
methods  base  their  estimation  on  derivatives  of  the  output  measurements  [168]  [169]. 
Here  again,  it  is  shown  that  input  estimates  can  be  determined  on  a  larger  class  of 
systems  from  taking  derivatives  on  the  partial  state.  The  remainder  of  this  section 
addresses  the  necessary  condition  for  input  estimation  and  how  the  output  equations  can 
be  chosen  so  that  it  can  be  satisfied.  Similarities  between  the  unknown  input  observer 
with  input  reconstruction  and  system  inversion  will  also  be  addressed.  Moylan  [268] 
showed  a  method  for  determining  the  inverted  system  in  state  space.  It  will  be  shown 
here  how  this  result  relates  to  a  system  in  right  polynomial  transfer  matrix  form  through 
the  partial  state.  Finally,  two  DC  ZEDS  examples  will  be  shown  to  demonstrate  the 
concepts. 

1.  Input  Observability  Theory 

The  goal  of  the  input  reconstruction  problem  is  to  formulate  a  system  filter  whose 
inputs  are  the  measurements  and  known  inputs  of  the  original  system  and  whose  outputs 
converge  to  the  inputs  of  the  original  system.  We  begin  with  a  short  example.  Suppose 
we  have  a  third  order,  single  input,  2-output  system  as  shown  in  Figure  III. 4  and  is 
represented  by  the  polynomial  matrix  fraction  H(p)  =  B(p)A  l(p) .  If  a  UIO  is 
developed  (not  necessarily  polynomial-based)  such  that  a  partial-state  vector  estimate 
z_(t )  is  determined,  then  by  the  polynomial  matrices  B(p)  =  [Bx(p)  B-,(p)f  and  A(p) 
defined  as  in  (0.27)  and  (0.28),  an  estimate  of  the  input  u{t)  and  output  y(t )  can  be 
determined  from  (0.29)-(0.30)  as 


|ti(0=  A  H*)+  ,6,z(o+  ,h2z(0+  AAO 

|y2(0  =  2b0z(t)+  2^2(0+  2h2z(0+  2h3z(0 


(0.66) 
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u(t)  =  z(t)  +  axz{t)  +  a2z(t )  +  a2z(t). 


(0.67) 


Let  yx(t)  be  arbitrarily  chosen  to  be  the  “input”  to  the  estimation  u(t) ,  meaning  that  we 
solve  for  the  highest  degree  term  in  (0.66)  and  substitute  in  (0.67)  to  obtain 


u(t)  =  -j-((a1  -  A  ) z(t)  +  (a2  -  xb2 )  z(t)  +  (a3  -  A  )z)  +  -j-j^f) 


A 


A 


(0.68) 


provided  that  xb0  is  nonzero.  This  constraint  implies  that  the  polynomial  Bt(  p)  must  be 
the  same  order  as  A(p ) .  In  the  state  space  realization,  this  means  that  the  matrix  D  in 
the  state  space  model  must  be  nonzero.  If  xb0  is  zero,  then  the  input  can  be  estimated 

directly  from  (0.67)  by  taking  the  derivative  on  the  highest  order  in  z_{t )  to  obtain  z  (/) . 
Notice  that  the  equation  (0.68)  is  similar  to  the  system  inversion  process,  since  the 
“output”  is  the  system  input  estimate  u(t )  and  the  “input”  is  the  measurement  _y,  (/) 
[268]  [169]. 

Now  consider  the  more  general  MIMO  case,  again  represented  in  state  space  as 

x(t)  =  Ox(t)  +  rw(t)  (0.69) 

y(t )  =  Cx(t )  +  Du(t )  (0.70) 


where  x(t)  e  K"1 ,  u(t )  e  M"" ,  y(t  )  e  M"v  and  the  constant  matrices  O  ,  T ,  C ,  and  D  are 
of  appropriate  dimensions.  If  the  system  is  a  minimal  realization  and  a  UIO  is  designed 
to  estimate  the  partial-state  vector,  then  similar  to  solving  for  the  high  order  derivative  in 
(0.66)  for  the  SIMO  case,  we  use  the  result  in  (0.36)-(0.37)  to  relate  the  input  estimate 
u{t )  in  terms  of  the  coefficients  of  the  right  matrix  fraction  polynomials  B(  p) ,  A ( p )  and 
the  output  y(t)  as 
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(0.71) 


m  =  -(BA' )'  ( B,  -  ba'aW) +{bA'  )f  A  0 


where  Ah  e  R""x”“  and  B,  e  M”vX""  are  constant  matrices  of  high  degree  coefficients  of 
A ( p )  and  B( p) ,  respectively.  Likewise,  Aj  e  and  Bt  e  M”1’*”1  are  the  matrices  of 

all  lower  order  coefficients  of  A(p)  and  B( p) ,  respectively.  The  superscript 
represents  the  Moore-Penrose  pseudoinverse.  For  more  details  on  this  representation,  see 
Appendix  D.  Clearly,  in  order  to  have  an  inverse,  Ah  must  be  nonsingular  and  the 

product  BhAhx  (which  is  the  D  matrix  in  state  space)  cannot  be  zero  for  a  nontrivial 
solution.  Again,  a  nonzero  Bh  implies  that  B(  p)  must  have  the  same  column  degree  as 
A(p) .  If  it  does  not,  then  the  input  estimate  can  be  found  from  (0.36)  as 


u(t)  =  Ahzfi(t)  +  A,z(t) 


(0.72) 


where  z  (t)  = 


p^'zx(t) 


is  the  vector  of  p\h  n 


derivatives  on  the  partial 


state  z(t) ,  determined  by  any  standard  derivative  filter.  We  have  shown  in  Section 
IV.B.l  that  in  cases  where  the  UIO  rank  condition  is  not  satisfied,  the  partial  state  can 
still  be  estimated.  This  result  can  be  extended  to  determining  an  input  estimate  by  adding 
all  derivatives  of  all  partial  states  up  to  the  column  degree  //(,  i  =  \,...,nu. 


Input  estimation  is  very  closely  related  to  system  left  inversion  where  a  system’s 
inputs  are  detennined  from  knowledge  of  its  output.  Necessary  and  sufficient  conditions 
for  the  stable  inversion  problem  are  provided  in  [268].  Hou  and  Patton  show  the 
similarities  between  input  estimation  and  system  inversion  and  provide  a  necessary  and 
sufficient  condition  for  input  observability  [169]. 
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For  both  polynomial  methods  herein,  special  attention  is  paid  to  the  input 
reconstruction  problem  because  it  provides  more  analytical  redundancy  and  potential  for 
improving  a  Fault  Tolerant  Control  architecture. 

In  order  to  determine  conditions  for  input  observability,  consider  the  LTI  model 
in  (0.69)-(0.70)  and  express  the  system  as 


-pi + o  r 

x(t) 

"0" 

C  D 

u(t) 

I 

y(t)  • 


(0.73) 


Now  define  the  matrix 


M(/t)  = 


-XI  +  cD 

C 


r 

D 


(0.74) 


called  the  matrix  pencil  [200]  [139].  Input  observability  is  a  necessary  and  sufficient 
condition  for  the  existence  of  an  estimator  reconstructing  inputs  [169].  In  the  case  of  a 
controllable  and  observable  system  of  equations  (0.69)-(0.70)  with  unknown  initial 
conditions,  the  input  u{t )  is  observable  if  and  only  if 


rank  (M(A))  =  dim(x)  +  dim(i/),  \/X  e  C .  (0.75) 


If  the  input  observability  condition  (0.75)  is  satisfied  for  a  system,  then  one  such  input 
estimator  exists  of  the  fonn 


w(t)  =  Ftw(t)  +  F2y(t) 


(0.76) 
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(0.77) 


where  yil\t )  denotes  the  i'h  derivative  of  y(t )  [169].  The  number  k  is  called  the  index 
of  the  estimator  and  it  indicates  the  highest  derivative  needed.  The  methods  for  input 
estimation  described  here  and  elsewhere  require  high  orders  of  derivatives  on  the  output 
measurements,  which  would  make  the  estimate  sensitive  to  noise.  Also,  in  [268]  an 
iterative  algorithm  is  described  for  determining  an  input  estimation  filter. 

Another  method  for  determining  the  input  estimate  is  derived  directly  from  the 
state  space  realization  (0.69)-(0.70)  and  is  similar  to  the  system  inversion  problem  in 
[268],  only  no  iterative  process  is  necessary.  It  is  expressed  here  in  partial-state  space 
form  as 


i(0  =  T  x  ( ®-TEC)Tz(t)  +  T  xrEy(t )  (0.78) 


u(t)  =  -ECTz_(t)  +  Ey(t ) 


(0.79) 


where  all  the  system  matrices  are  defined  in  (0.69)-(0.70),  the  matrix  E  is  defined  as 


\D~l  ^ nU=ny 

|d;  if  ny>nu. 


and  the  nonsingular  matrix  T  transforms  state  x(t)  to  the  partial-state  vector  lit)  as 
before.  Notice  that  each  tenn  in  the  input  estimation  equation  (0.79)  includes  the  matrix 
E ,  which  is  directly  related  to  the  degree  of  the  numerator  polynomial  B(p )  as  shown  in 
the  equivalent  equation  (0.71)  from  the  polynomial  approach.  Once  again,  in  order  to 
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have  a  proper  system  inverse,  the  degree  of  the  numerator  polynomial  B(p)  must  be  the 
same  as  the  denominator  polynomial  A(p)  or,  in  the  state  space  realization  D  ^  0 . 

In  most  cases  for  fault  diagnosis,  the  number  of  outputs  will  be  strictly  greater 
than  the  number  of  inputs.  This  implies  that  D  is  not  square  and  only  some  of  the 
eigenvalues  of  (A-BEC)  will  depend  on  E .  If  D  is  rank  deficient,  then  the  inputs 

corresponding  to  a  zero  column  will  not  have  an  estimate.  Notice  that  the  inverted 
“output”  equation  in  (0.79)  is  identical  to  (0.71),  which  again  shows  the  equivalence  of 
minimal  realizations  of  a  system. 

If  the  system  does  not  satisfy  the  inversion  rank  condition,  it  does  not  mean  that 
an  input  estimate  cannot  be  estimated.  It  only  means  that  the  system  does  not  have  a 
proper  inverse.  An  explanation  in  [73]  and  [280]  shows  that  if  we  relax  the  rank 
condition  to 


rank(M'(T))  <  dim(x)  +  dim(w)  (0.80) 

where 


M'{X) 


-x  /+o  r 

C  0 


and  we  further  assume  that  a  full-state  estimate  is  obtainable  and  input  matrix  T  is  full 
rank  (which  it  will  be  if  the  system  is  controllable),  then  an  input  estimation  process  can 
be  established  with  derivatives  on  the  output  measurements  as  shown  in  (0.76)-(0.77). 
Another  way  to  look  at  this  result  is  to  consider  an  estimate  of  the  input  directly  from  the 
dynamic  equation  (0.32)  represented  in  the  partial  state  as 


u(t) ={t1  r)f  i(t)  -  (r_Ir)  T'®Tz(t) . 


(0.81) 
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By  comparing  equations  (0.79)  and  (0.81),  we  see  that  the  expression  (0.81)  shows 
dependence  on  derivatives  of  the  partial  state,  but  is  independent  of  D  .  This  relationship 
shows  the  similarity  of  determining  the  input  estimate  based  on  the  right  fraction  transfer 
matrix  representation  in  (0.72).  Finally,  recall  that  an  estimate  of  the  partial  state  is 
possible  even  when  a  full-state  estimate  is  not.  Then  by  extension,  an  input  estimate  in 
this  case  can  still  be  achieved  by  adding  derivatives  as  necessary  to  (0.81). 

2.  Examples 

Two  examples  are  presented  to  demonstrate  input  reconstruction  on  an  unknown 
input  observer.  The  first  shows  how  a  bank  of  unknown  input  observers  successfully 
estimates  the  inputs  while  being  systematically  independent  of  each  input.  From  the 
structured  residual  set,  it  can  be  seen  that  the  indication  of  the  presence  of  a  fault  is  more 
prominent  using  residuals  from  the  input  estimate  rather  than  based  on  the  output 
measurements.  In  addition,  for  a  certain  class  of  systems  the  residual  signal  is  an 
estimate  of  the  fault  signal  itself.  The  second  example  demonstrates  an  unknown  input 
observer  used  to  estimate  nonlinear  components  of  a  system. 

The  first  example  is  a  UlO-based  structured  residual  set  designed  to  detect  and 
isolate  actuator  faults.  Suppose  a  system  follows  the  flow  diagram  depicted  in  Figure 
III. 8.  The  system  model  shows  an  input  u(t)  being  subjected  to  an  unknown  disturbance 

rju  with  known  mean  (zero)  and  variance  ( cr  j  and  another  completely  unknown 
disturbance  fit).  The  corrupted  signal  lift)  is  the  input  to  the  known  plant  dynamics 
and  its  output  measurement  is  also  influenced  by  a  zero  mean  disturbance  with  variance 
a] .  The  goal  is  to  determine  the  characteristics  of  the  unknown  disturbance  / (/)  -  that 

is,  determine  / (t) . 
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Figure  III.8:  UIO  Example,  Signal  Flow  Diagram 

To  make  the  problem  more  general,  suppose  the  system  depicted  in  Figure  III. 8  is 
a  MIMO  system  where  u(t)  and  y(t)  are  vectors.  Let  /(/)  be  a  scalar  disturbance 
operating  on  one  of  the  inputs,  but  it  is  not  known  a  priori  which  one.  Also  let  the 
disturbances  qu  and  rj  have  influence  on  all  the  inputs  and  outputs,  respectively.  Now, 

the  goal  is  to  detennine  which  input  is  corrupted  by  f(t)  on  the  basis  of  the  input 
residual 


ru(t)  =  u(t)-u(t ) 


and,  for  comparison,  the  output  residual 


ry(t)  =  y(t)-y(t). 


The  problem  is  solved  using  unknown  input  observers  with  input  reconstruction 
(UIO-IR)  in  a  structured  residual  set  depicted  in  Figure  III. 9  where  the  bank  of  observers 
are  designed  to  reject  the  inputs  one  at  a  time.  The  system  used  in  this  example  is  the 
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sixth-order,  four-input,  twelve-output*  port  bus  model  for  the  DC  ZEDS  system  shown  in 
Figure  II. 6.,  with  parameters  listed  in  Table  II. 3.  For  the  purpose  of  this  exercise,  the  bus 
is  assumed  to  be  in  Mode  3  (where  zones  1,  2,  and  3  are  energized)  and  the  system  is 
adjusted  in  time  scale  using  arbitrary  inputs 


sin(r) 

Ux  it) 

7T 

u2(t) 

sin(r)  +  — 

w3  (t) 

sin(r)  +  n 

UA{t)_ 

...  n 

sm(0~- 

Figure  III.9:  UIO  Example,  Structured  Residual  set  for  Input  FDI 


It  turns  out  that  twelve  outputs  are  not  necessary  for  the  UIO  implementation.  A  more  manageable 
number  is  9  (bus  current,  bus  voltage,  and  capacitance  voltage  in  each  zone),  which  still  meets  all  criteria 
for  the  minimal  realization,  UIO  rank  condtion  and  input  observability  rank  condition. 
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The  bank  of  UIOs  (designed  using  the  polynomial  parametrization  method)  in 
Figure  III. 9  shows  that  each  of  the  four  observers  is  designed  to  be  insensitive  to  ux(t), 

u2(t) ,  u3 (t) ,  and  u4(t) ,  respectively,  where  the  represents  the  unknown  input.  The 
other  inputs  to  the  observer  include  the  outputs  y  (denoted  as  y f(t)  in  Figure  III. 8)  less 

the  bus  voltage  sensors  where  indicated.  For  this  problem,  the  system  is  a  minimal 
realization  and  the  UIO  rank  condition  is  satisfied  for  all  conditions.  If  only  the  bus 
current,  bus  voltage,  and  capacitance  voltage  for  each  zone  are  considered  as  outputs, 
then  the  D  matrix,  in  (0. 1 1)  would  be 


D  = 


0 

0 

0 

0 

0 

0 

-0.1 

0 

0 


0  0  0 

0  0  0 

0  0  0 

0  0  0 

0  0  0 

0  0  0 

0  0  0 

-0.1  0  0 

0  -0.1  0 


and  the  input  observability  rank  condition  in  (0.75)  is  satisfied.  From  this  structure,  it  can 
be  seen  that  rank(D)  =  3  <  nu  and  so  a  solution  to  the  input  estimate  can  be  determined 
by  (0.71)  or  (0.79)  for  all  except  u4(t ) ,  which  shows  all  zeros  in  the  fourth  column.  In 
other  words,  none  of  the  outputs  are  algebraically  dependent  on  the  input  u4(t) ,  which  is 
the  output  voltage  from  the  power  supply  in  the  DC  ZEDS  system. 

Also,  since  there  is  input  dependence  on  the  last  three  output  measurements,  they 
must  be  removed  from  the  observer  estimate  because  a  fault  on  the  input  will  also  mean  a 
fault  on  the  outputs  that  are  linearly  dependent  on  the  inputs.  These  last  three  rows  of  D 
correspond  to  the  bus  voltage  sensors. 
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In  this  example,  a  disturbance  signal  / ( t )  was  arbitrarily  chosen  to  be  a  triangle 
waveform  added  to  the  input  ux{t)  as  shown  in  Figure  III.  10. 


^(0 


Figure  III.10:  UIO  Example,  Input  Fault  Trajectory 


The  results  from  the  simulation  were  obtained  through  resuiduals  generated  from 
the  structured  residual  set  in  Figure  III. 9  and  show  a  fault  indication  on  the  far  left 
column  of  Figures  III.  11  -  III.  13.  The  columns  of  these  figures  represent  residuals 
generated  from  the  UIOs  insensitive  to  m,  (t) ,  u2(t) ,  u3(t) ,  and  u4 (t)  from  left  to  right, 

where  the  represents  the  unknown  input.  The  rows  in  each  figure  are  residuals  based 
on  estimates  of  the  inputs  or  outputs.  In  Figure  III.  11,  the  subscript  “ipz*”  are  the  load 
demands  leading  to  the  SSCMs  in  each  zone  as  shown  in  Figure  II. 6  as  “  iouU  ”,  in  Figure 


III.  12,  the  subscript  “ibus*”  are  the  bus  currents  in  each  zone  and  in  Figure  III.  13,  the 

subscript  “vb*”  is  bus  voltage  in  each  zone.  The  fourth  input  is  vout  from  the  power 
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supply  but  is  not  shown  in  Figure  III.  11  because  an  estimate  cannot  be 
determined  without  derivatives  to  the  partial  state. 

It  can  be  seen  that  the  residual  in  the  very  top  left  corner  of  Figure  III.  11  is  the 
estimate  of  the  disturbance  f(t )  and  because  it  resides  in  the  first  column,  where  no 
other  signatures  reside  going  down  the  column,  the  disturbance  is  on  the  first  input  ux{t) . 
To  make  the  figures  more  readable,  the  triangle  disturbance  signal  f(t)  is  superimposed 
in  red,  the  actual  results  are  in  blue,  and  the  mean  is  green. 


U-234 


U1-34 


U12-4 


U123- 


-1 


-1 


-1 


. 

4 

1 

0 

4 

rrWT"  Ti 

1 

0 

4 

0  10  20  '[ 

1 

10  20  '[ 

1 

i  10  20  'c 

4 

10  20 

n 

n 

0 

n 

0 

0 

•1 

4 

4 

0  10  20  'c 

1 

10  20  '  C 

4 

i  10  20  'c 

4 

10  20 

0  10  2 

n 

0 

|U|XUuuJjL||UU|U| 

10  2 

n 

n 

10  20 

-i 

0  c 

u 

-i 

o  c 

10  2( 

u 

-i 

]  c 

Time 


Time 


Time 


Time 


FigureIII.ll:  Input  Residuals  (Fault  on  ux{t)) 


If  the  estimate  ux(t)  is  determined  from  UIO. U-234  (shorthand  for  the  unknown 
input  observer  insensitive  to  ux(t)),  then  from  the  consistency  between  input  ux(t) ,  u2(l) , 
u3(t) ,  and  yf(t) ,  an  input  estimate  will  be  u(t)  =  ux(t)  +  f  (t)  and  the  residual  xru 
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(shorthand  for  input  residual  from  UIO.U-234)  would  be  xr  =  -/,'(/) ,  which  is  the 

opposite  of  the  estimate  of  the  disturbance  signal.  Any  other  combination  of  inputs  and 
outputs  will  result  in  some  signature  in  the  residuals. 

Looking  at  the  other  columns  in  Figure  III.  11,  it  is  possible  to  have  ambiguous 
results.  For  example,  column  U 1 2~4  shows  a  similar  signature  in  the  top  row  with  very 
small  residual  signatures  down  the  column.  This  is  essentially  a  problem  of  residual 
evaluation,  but  a  solution  might  be  that  if  the  disturbance  signal  (i.e.  the  fault)  is  known 
to  have  a  certain  characteristics,  then  a  matched  filter  can  be  implemented  along  with  any 
standard  noise  reduction  scheme  to  identify  signals  buried  in  noise.  Another  method 
might  be  to  look  for  fault  characteristics  in  the  frequency  domain. 
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Figure  III.12:  UIO  Example,  Output  Residuals  (Bus  Current  per  Zone) 
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Looking  at  Figure  III.  12  and  Figure  III.  13,  we  see  the  residuals  based  on  output 
measurements.  Once  again,  the  very  first  column  has  no  signature  while  all  the  other 
columns  do  (except  for  ]rv/)|  ,  which  will  be  explained  later).  This  means  that  the  fault 

indication  is  associated  with  u]  (t) .  Notice  in  these  two  figures,  the  residuals  are  not 
nearly  as  pronounced  as  they  were  in  the  residuals  based  on  input  estimation.  This  can  be 
explained  by  the  linear  combination  of  matrices  in  generating  the  residuals  from  the  input 
estimation  is  better  “tuned”,  or  simply  more  suitable,  for  detecting  input  faults.  Likewise, 
residuals  generated  from  output  estimates  are  better  suited  for  sensor  FDI. 

As  mentioned  before,  the  only  output  residual  with  any  signature  in  the  first 
column  is  the  one  in  the  top  left  corner  of  Figure  III.  13  ( ,  rvh^  j .  This  is  because  the  bus 
voltage  sensor  in  zone  1  is  linearly  dependent  on  the  input  w,  (?) .  Recall  that  this 
measurement  was  removed  from  UIO.U-234  when  determining  the  state  estimate.  This 
means  that  the  estimate  will  have  a  u]  (t)  fault  influence  on  the  output  measurement  and 

therefore  the  row  can  be  discarded. 
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Figure  III.  13:  UIO  Example,  Output  Residuals  (Bus  Voltage  per  Zone) 


In  the  second  example  the  purpose  is  to  demonstrate  how  a  UIO  can  be  used  to 
estimate  the  nonlinear  terms  in  a  system  and  then  further  estimate  unknown  input.  In  this 
case,  the  power  supply  model  with  buck  controller  is  used  where  its  topology  and 
corresponding  parameters  are  shown  in  Figure  II. 3  and  Table  II.l.  Also,  its  controller 
design  and  parameters  are  shown  in  Figure  II.4  and  Table  II.2.  The  power  supply  model 
dynamic  equation  is  shown  in  (0.4)  and  rewritten  here  as 


x  =  Apsx  +  B 


1  ,ps 


«l  (0 

+ 

u2(t) 

_M  4(0_ 

2,ps 

_m3(0_ 

(0.82) 
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where  the  known  inputs  are  [u^t)  w4(/)]7  =  42  E  iLoad  and  the  “unknown”  inputs 

are  the  nonlinear  terms  [ u2(t )  u3(t)]T  =  [j -iL  d-vdc J  .  The  UIO  design  shown  in 

Figure  III.  14  describes  the  estimation  procedure.  The  idea  is  to  design  an  observer 
independent  of  the  nonlinear  dynamics,  then  estimate  the  inputs  through  input 
reconstruction.  Finally,  divide  the  unknown  input  estimates  by  the  state  estimates  vdc 

and  iLout  to  detennine  an  estimate  on  the  control  input  d(t) .  If  the  input  can  be 
measured,  then  a  residual  rd(t)  can  be  formed.  In  this  case,  both  the  UIO  rank  condition 
and  the  input  observable  condition  are  satisfied. 


Figure  III.14:  Input  Estimation  of  a  Nonlinear  System 


The  results  are  shown  in  Figure  III.  15.  Here,  the  estimate  of  the  control  input 
d(t )  is  detennined  in  the  top  graph  while  the  bottom  graph  only  serves  to  show  that  the 
estimate  found  is  consistent  with  the  actual  control  input  trajectory. 
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Duty  Cycle  Estimate 


Time  (sec) 


Figure  I1I.15:  Results  from  Polynomial  UIO  with  Input  Reconstruction 

From  the  observer  point  of  view,  if  algebraic  equations  relating  the  unknown 
input  to  the  measured  output  are  included  in  the  output  measurements  and  certain  rank 
conditions  similar  to  the  left  inverse  problem  in  [268],  then  input  estimation  is  achievable 
without  additional  derivatives  on  the  partial  state. 
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IV.  FAULT  DIAGNOSIS  APPLIED  TO  THE  DC  ZEDS  SYSTEM 


It  has  been  shown  in  numerous  publications  that  the  physics-based  models  of  the 
DC  ZEDS  system  accurately  represent  system  dynamics  [285]  [335] [327] [328] [336]. 
Indeed,  the  DC  ZEDS  simulation  model  used  herein  was  developed  and  tested  by  the 
authors  of  the  aforementioned  publications  and  the  text  [220].  Because  the  Power 
Supply,  bus,  SSCM,  and  Oring  functions  in  Figure  II.  1  and  described  in  Chapter  II  are 
well-defined,  it  is  reasonable,  then,  to  take  a  model-based  approach  to  Fault  Diagnosis  on 
the  DC  ZEDS  system.  Chapter  II  sought  to  organize  the  DC  ZEDS  system  into  modules 
and  represent  their  dynamics  in  state  space  form.  In  cases  where  the  average  value 
models  are  nonlinear,  such  as  the  power  supply  and  SSCM  modules,  the  models  are 
separated  into  linear  and  nonlinear  parts.  Chapter  III  then  developed  several  unknown 
input  observer  approaches  that  assume  a  model  in  state  space  form  for  implementation. 

The  nature  of  the  DC  ZEDS  system  provides  high  bandwidth  system  control 
response,  but  because  it  can  be  controlled  so  tightly,  faults  can  propagate  quickly  as  well 
[327]  (on  the  order  of  tens  of  milliseconds).  The  fault  diagnosis  schemes  must 
distinguish  between  normal  and  casualty  transients  and  ultimately  be  able  to  dynamically 
reconfigure  under  stressful  battle  conditions  [390].  Because  time  is  critical,  model  based 
methods  to  detect  the  source  of  a  fault  rather  than  trace  its  effect,  is  sought. 

This  chapter  applies  a  variety  of  model-based  observer  techniques  (including  the 
polynomial  approach)  to  the  DC  ZEDS  simulation  model  for  fault  detection  and  isolation. 
It  is  desired  to  detect  faults  and  classify  them  according  to  their  characteristics.  Some  of 
the  means  for  characterization  depends  upon  the  detection  method  (i.e.  structured  residual 
set,  unknown  input  observers,  or  fault  modeling),  while  others  upon  the  signature  itself 
(i.e.  positive  or  negative,  shape,  or  relative  size).  The  goal  for  this  chapter  is  to  develop  a 
set  of  methodologies  to  identify  residuals  that  indicate  certain  classes  of  faults. 

It  is  acknowledged  here  that  standalone  sensors  perform  a  certain  level  of 
analytical  redundancy  without  computations.  However,  individual  sensors  do  not 
contribute  to  a  “global”  perspective  as  do  a  network  of  sensors  connected  together  by 
mathematical  relationships.  One  of  the  presuppositions  in  this  work  is  to  develop  ways  to 
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automate  the  fault  detection  process  quickly,  reliably,  and  accurately  -  and  a  hierarchical 
automated  fault  detection  layer  using  robust  model-based  estimation  is  one  such  way. 


A.  MODULAR  APPROACH  TO  INPUT-OUTPUT  CONSISTENCY 

The  DC  ZEDS  model  as  a  whole  contains  68  states,  and  this  does  not  consider  the 
ac  side  of  the  overall  power  system  from  prime  mover  to  constant  power  load.  From  the 
input  of  the  power  supply  to  the  input  to  the  SSIM  in  each  zone,  the  DC  ZEDS  system 
has  many  modes  and  system  configurations  that  make  it  difficult  to  estimate  as  a  whole  - 
even  under  normal  operating  conditions.  A  suitable  solution  to  handle  the  sheer  size  of 
the  problem  is  to  separate  the  system  into  smaller  logical  subsystems  and  connect  them 
together  with  common  inputs  and  outputs.  This  modular  approach  to  verifying  input- 
output  consistency  is  done  for  the  DC  ZEDS  system  and  the  topology  is  shown  in  Figure 
IV.  1. 


Figure  IV.  1:  DC  ZEDS  Module  Topology 
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Four  representative  component  module  block  types  are  configured  so  that  each 
input/output  pair  between  them  is  shared  measurements  such  that  the  output  to  one  block 
is  the  input  to  the  connecting  block.  Under  this  method,  the  entire  system  is  viewed  as  a 
dynamic  interconnection  of  subsystems  with  common  inputs  and  outputs  and  is 
consistent  with  the  DC  ZEDS  simulation  model.  Notice  also  in  Figure  IV.  1  that  the 
models  show  voltages  as  outputs  moving  toward  the  constant  power  load  in  each  zone 
and  currents  as  inputs. 

Using  various  types  of  observers,  model  consistency  is  verified  by  comparing 
module  estimates  with  measurements.  With  algebraic  equations  relating  the  states  and 
the  input  in  the  fonn  of  energy  conservation  laws,  the  observers  become  input  observable 
and  therefore  an  estimate  of  the  inputs  can  be  made  without  derivatives.  Taking  a  closer 
look  at  the  input/output  connections  shown  in  Figure  IV. 2,  we  see  that  redundancy  is 
verified  from  the  input  and  output  measurements  of  each  module  block  “A”,  “B”,  and 
“C”. 


Figure  IV.2:  System  Module  Interconnection 


Suppose  block  “B”  is  the  port  SSCM  in  Zone  1  where  the  external  inputs  and 
ufh  in  (0.19)-(0.20)  are  the  port  bus  voltage  and  load  current  and  let  yh  and  yh  be  the 

SSCM  output  voltage  and  the  Zone  1  load  requirement  to  the  port  bus,  respectively.  To 

verify  analytical  redundancy  across  the  entire  DC  ZEDS  system,  estimates  can  be  made 
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on  system  measurements  from  both  sides  of  the  block  and  compared  to  the  estimates  (or 
measurements)  from  each  connecting  block.  In  our  example,  un  and  y  from  the  port 

bus  (block  “A”)  would  represent  the  bus  voltage  and  load  current  in  (O.lO)-(O.l  1),  while 
u  and  y  from  the  Oring  function  (block  “C”)  would  represent  input  vport  and  output 

i  rt  in  (0.2 1  )-(0.22).  Residuals  can  be  generated  by  comparing  the  measurements  from 
each  common  line  as 


rAB(t)  =  ^(0-^(0 
rBA(t)  =  {la2  (0-A2(0 

nc(0  =  a,  (t)-y^(t) 

rc«(0  =  ^2(0-A2(0 


Notice  that  the  one  estimate  in  each  equation  of  (0.83)  is  independent  of  the  other. 
In  this  manner,  the  modular  interconnection  model  shown  in  Figure  IV.2  can  zoom  down 
to  any  desired  level  in  the  system,  provided  its  external  inputs  and  outputs  are  the  outputs 
and  inputs  to  connecting  modules. 

Among  the  benefits  of  this  consistency  scheme  is  increased  redundancy.  For 
example,  if  uu  is  a  control  signal  input  for  the  controller  in  Module  A,  and  the  signal 

becomes  disconnected  or  otherwise  erroneous,  then  between  the  two  estimates  u„  and 

“2 

y>h  the  fault  can  be  isolated  to  one  of  the  modules.  In  some  cases,  different  methods  of 

FDI  can  be  performed  in  separate  modules  that  draw  the  same  conclusion  and  your 
confidence  of  the  fault  indication  becomes  much  stronger. 

From  Module  A,  verifying  sensor  consistency  is  done  using  a  structured  residual 
set  of  regular  Kalman  Filters  as  shown  in  Figure  IV. 3.  Here,  an  observer  is  designed  to 
minimize  the  mean  squared  error  of  the  output  residuals  using  all  the  inputs  and  all  the 
outputs  less  one.  Each  successive  observer  in  the  set  perfonns  the  same  operations  until 
all  n  sensors  are  removed  one  at  a  time  for  a  total  of  nv  observers.  The  notation  ly{t) 
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indicates  the  output  measurement  vector  y(t )  without  sensor  1,  and  so  on.  Likewise, 
j  y(t)  and  r  (/)  are  the  corresponding  output  measurements  and  output  residuals  based 
on  observer  inputs  u(t)  and  y(t)  without  sensor  1.  The  residual  generator  compares  the 
full  output  vector  y(t)  with  all  its  estimates  for  consistency.  Suppose  for  example  that 
the  number  three  output  sensor  is  faulted. 


v(0 


Figure  IV.3:  Sensor  FDI  Structured  Residual  Set 

Then  by  the  observer  structure  in  Figure  IV.3,  the  third  observer  in  the  set  would 
accurately  estimate  the  faulted  output  signal  3,  while  the  rest  would  show  a  signature  in 
the  residuals.  This  is  because  all  the  other  observers  are  dependent  on  the  faulted 
measurement. 
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The  following  DC  ZEDS  example  demonstrates  the  concept  of  modular  FDI. 
Consider  the  nominal  operating  scenario  where  the  port  and  starboard  power  supplies  are 
energizing  the  port  and  starboard  bus,  respectively.  The  starboard  bus  is  providing  power 
to  the  load  in  Zone  1  while  the  port  bus  is  feeding  the  load  in  Zones  2  and  3  with 
operating  parameters  according  to  Table  II.  5.  The  power  supply  controller  is  regulating  a 
500  V  output  voltage  through  the  control  input  d  ,  which  is  assumed  to  be  measurable. 

Suppose  a  fault  occurs  to  the  load  current  signal  iout  in  the  port  power  supply 
controller  (shown  in  Figure  II. 4)  at  time  t fauh  =0.6  seconds.  In  the  fault  condition  the 
signal  iout  is  a  constant  zero,  thus  removing  the  feedforward  path  to  the  PI  controller. 

The  local  effect  of  the  fault  is  that  the  feedforward  path  will  communicate  to  the 
controller  that  there  is  no  required  load  current,  which  disables  the  controller’s  ability  to 
provide  fast  response  to  changes  in  the  load.  The  fault  propagates  through  to  the  port  bus 
as  a  small  transient  that  quickly  recovers  to  the  steady  state  as  shown  in  Figure  IV.4  and 
Figure  IV. 5.  The  starboard  bus  remains  unchanged. 
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Figure  IV.4:  Bus  Voltage  Showing  Bus  Current  Sensor  Fault 
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Starboard  Bus:  Bus  Current,  Zone  1 


Figure  IV.5:  Bus  Current  Showing  Bus  Current  Sensor  Fault 

The  goal  here  is  to  use  the  modular  FDI  approach  to  detect  the  fault  and 
determine  in  what  module  the  fault  originated.  Beyond  that,  it  is  desirable  to  gain  other 
fault  characteristics,  if  possible.  Looking  at  the  power  supply  module  shown  in  Figure 
IV.  1  we  see  that  it  is  connected  to  the  port  bus  through  common  signals  vout  and  iLoad . 

The  first  step  to  our  goal  is  to  conduct  sensor  FDI  on  the  sensors  within  the  bus  module. 
Developing  a  structured  residual  set  of  Kalman  Filters  as  shown  in  Figure  IV. 3  achieves 
that  task.  So  as  to  not  be  unnecessarily  exhaustive  in  showing  all  the  sensor  FDI 
residuals,  only  the  zone  current  sensors  are  plotted  here  since  the  other  residuals  are 
consistent  with  what  is  shown  here  in  Figure  IV. 6. 
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Figure  IV.6:  Bus  Sensor  FDI  Residuals  for  Bus  Current  Sensor  Fault  in  Power 

Supply 

The  residuals  in  Figure  IV.6  show  that  there  is  consistency  between  estimates  and 
measurements  since  they  are  white  and  stationary.  It  can  be  shown  that  the  other  sensors 
in  the  set  show  the  same  result.  Based  on  this  result,  we  conclude  that  the  sensor  fault  is 
not  in  the  bus  module.  Now  we  take  a  look  at  the  structured  residual  set  in  the  port  power 
supply  module.  Since  the  fault  in  question  occurs  on  the  input  to  the  power  supply,  an 
actuator  fault  diagnosis  problem  similar  to  the  one  described  in  Figure  III. 8  applies. 
Applying  the  structured  residual  set  similar  to  the  one  shown  in  Figure  III. 9  for  the  power 
supply  model  in  (0.4)-(0.5),  we  receive  the  following  residuals  on  the  input  and  output 
estimates  as  shown  in  Figures  IV. 7  and  IV. 8. 
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Figure  IV.7:  Input  Residuals  for  Bus  Current  Sensor  Fault  in  Power  Supply 

Controller 


In  contrast  to  the  sensor  residuals  in  the  bus  module  in  Figure  IV. 6,  the  residuals 
in  Figure  IV.7  indicate  the  presence  of  a  fault.  As  before,  the  residuals  are  generated 
from  the  input  estimations  indicated  by  row  and  the  unknown  input  observers  by  column. 
Looking  down  column  U123-,  it  can  be  seen  that  the  only  residual  which  is  not  white  is 
in  the  last  row,  which  corresponds  to  the  signal  in  question,  iLoad .  This  makes  sense 
because  the  estimates  in  the  column  to  the  far  right  base  their  estimate  on  all  inputs 
except  iLoad ,  which  is  where  the  fault  resides.  The  residuals  down  the  column  indicate 

input/output  consistency,  while  the  input  estimate  for  iLoad  contains  the  fault  signature. 
Similar  results  are  shown  in  the  output  residuals  in  Figure  IV. 8. 
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Figure  IV.8:  Output  Residuals  for  Bus  Current  Sensor  Fault  in  Power  Supply 

Controller 


The  result  of  the  bus  current  sensor  fault  example  shows  that  first,  the  fault  did 
not  have  any  indication  of  originating  in  the  bus  module.  It  did,  however,  have  an 
indication  in  the  power  supply  model  and,  more  specifically,  the  fault  showed  indication 
of  residing  in  the  bus  current  sensor  somewhere  between  the  bus  and  the  controller. 


B.  FDI  BASED  ON  FAULT  MODELING 

In  the  last  section  we  considered  a  class  of  faults  where  the  fault  information  was 
due  to  a  difference  between  measurements  and  estimates.  Using  the  methodical, 
structured  residual  set  for  both  input  and  output  residuals,  a  fault  was  detected.  The  key 
component  to  this  method  is  a  valid  model  that  accurately  represents  the  input-output 
relationship.  In  this  section,  we  consider  another  means  to  detect  faults  using 
mathematical  models  specifically  “tuned”  to  a  particular  class  of  fault.  In  this  case, 
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regular  Kalman  Filters  are  determined  from  models  which  are  known  to  indicate  a  fault 
condition  in  the  DC  ZEDS  system.  In  this  case,  indication  of  a  fault  in  the  residuals  is 
exactly  opposite  of  the  results  from  the  previous  section.  A  zero  residual  would  indicate 
a  match  to  the  fault  model,  providing  precise  fault  detection  and  isolation  -  provided  the 
fault  model  is  consistent  with  the  actual  fault  dynamics. 

Suppose  we  have  a  nominal  plant  configuration  similar  to  the  one  described  in  the 
last  section  and  we  wish  to  consider  what  the  effects  of  an  SSCM  switch  failure  in  the 
closed  position  would  do  to  the  DC  ZEDS  system.  Then,  the  goal  is  to  model  this  failure 
and  simulate  the  FDI  scheme. 

For  simulating  the  SSCM  switch  in  the  closed  position,  we  first  choose  the 
starboard  SSCM  in  Zone  1  depicted  in  Figure  II.l.  In  the  average  value  model,  the  fault 
is  triggered  by  setting  the  SSCM  control  signal  d  to  1.  Similarly,  in  the  detailed  model, 
the  control  switching  signal  s  is  set  to  1.  A  simulation  of  the  fault  condition  was 
conducted  and  its  effect  on  the  bus  is  shown  in  Figures  IV. 9  and  IV.  10. 
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Figure  IV.9:  SSCM  Switch  Closed  Failure:  Port/Stbd  Bus  Voltage 


From  Figure  IV.9  it  can  be  seen  that  the  SSCM  failure  has  very  little  effect  on  the 
bus  voltage.  This  is  due  to  the  sophisticated  and  robust  power  converter  control  elements 
that  provide  nearly-ideal  load  regulation  and  transient  performance  [327] [332].  The  port 
and  starboard  bus  current  is  shown  in  Figure  IV.  10  where  the  starboard  bus  shows  a  large 
transient  at  the  time  of  the  fault,  but  settles  back  to  near  the  same  level  as  before.  It  can 
also  be  seen  on  the  starboard  bus  current  figure  that  the  average  value  model  follows  very 
closely  to  the  detailed  model,  even  in  the  fault  transient.  Notice  also  that  the  port  bus 
remains  unchanged. 
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Figure  IV.10:  SSCM  Switch  Closed  Failure:  Port/Stbd  Bus  Current 


The  next  set  of  figures  shows  the  inputs  and  outputs  to  the  SSCM  under  the 
switch  closed  fault  condition.  Figure  IV.  11  shows  the  inputs  to  the  SSCM.  The  bus 
voltage  is  virtually  unchanged  as  before,  but  the  load  current  demand  oscillates  very 
rapidly  in  the  steady  state.  This  is  due  to  the  constant  power  nature  of  the  load  that  the 
SSCM  is  feeding.  If  the  voltage  to  the  load  oscillates  in  the  steady  state,  then  so  will  the 
current  according  to  the  relationship  vLoad  ■  iLoad  =  P* ,  where  P  is  the  constant  load 
requirement  shown  in  Figure  11.15.  Here  again,  the  average  value  fault  simulation  agrees 
with  the  detailed  simulation. 


106 


600 


Bus  Voltage 


400  -I 


<U 

(13 


> 


r 

_i . 

— 

— 

— 

— 

- 1 - 1 - 1 - 

Detailed  Fault 

Detailed  No  Fault 
- AVM  Fault 

- AVM  No  Fault 

_ 

_ 

_ 

_ 

0  0.1  0.2  0.3  0.4  0.5  0.6  0.7  0.8  0.9  1 

Load  Current 


0  0.1  0.2  0.3  0.4  0.5  0.6  0.7 

Time  (sec) 


0.8  0.9 


FigureIV.ll:  SSCM  Switch  Closed  Failure:  SSCM  Inputs 


Looking  at  the  SSCM  circuit  diagram  in  Figure  II.  10,  it  can  be  seen  that  if  the 
switch  is  closed  permanently,  then  the  steady-state  output  should  be  very  close  to  the 
input.  This  is  indeed  the  case  as  shown  in  Figure  IV.  12  where  the  dc  value  of  the  SSCM 
output  is  close  to  the  bus  voltage  (500  V). 
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Figure  IV.12:  SSCM  Switch  Closed  Failure:  SSCM  Outputs 

Next  we  examine  what  is  happening  inside  the  controller  depicted  in  Figure  11.12. 
As  it  turns  out,  the  position  of  the  control  signal  d  becomes  important  for  detecting  the 
fault.  First,  it  makes  sense  that  the  duty  cycle  d  is  around  80%  under  no-fault  conditions 
because  the  regulated  output  voltage  is  80%  of  the  input  voltage. 

If  the  SSCM  switch  fails  closed,  then  we  see  that  the  output  voltage  becomes 
close  to  the  input  voltage.  The  controller  in  Figure  11.12  shows  a  large  error  signal 
between  the  output  voltage  (-500  V)  and  the  reference  voltage  (400  V),  resulting  in  the 
controller  commanding  the  switch  to  open.  This  fault  condition  means  that  the  command 
signal  d  is  a  constant  zero,  while  the  actual  operating  condition  of  the  switch  is  stuck  in 
the  opposite  state.  These  command  signal  dynamics  are  shown  in  Figure  IV.  13. 
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Buck  Control  Duty  Cycle 


Figure  IV.13:  SSCM  Switch  Closed  Failure:  SSCM  Control  Signal  d 


In  what  follows  we  derive  the  equations  that  model  the  fault.  Consider  the  SSCM 
model  equations  in  (0.19)-(0.20)  where  the  dynamic  equation  is  rewritten  here  as 


X  =  ASScmX  +  BU 


u. 

m9 

1 

2  ,sscm 

z 

U4_ 

u3 

(0.84) 


Recall  that  the  inputs  u2(t)  and  u2{t)  are  the  nonlinear  terms  relating  the  control  signal 
d  and  the  states  of  the  system  as 


u2  (t)  =  d  ■  x5  ( t ) 
u3(t)  =  d  -x3(0 
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(0.85) 


Under  the  SSCM  switch  closed  fault  condition,  the  control  signal  d  =  1  and  the 
nonlinear  tenns  in  (0.84)  can  now  be  expressed  as  a  part  of  the  linear  dynamics  in  the 
new  fault  dynamics  state  space  equation 


x (t)  =  A*  x(t)  +  B, 

V J  sscm  V  /  1  ,sscm 


«i  (0 
U4(t) 


(0.86) 


where 


0  0  0  0  0 
0  0  0  0  0 


A  —  A  + 
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The  fault  model  in  (0.86)  along  with  the  output  equations  in  (0.20)  is  controllable, 
observable,  and  input  observable.  Therefore,  a  standard  Kalman  Filter  is  applied  to  the 
SSCM  switch  closed  fault  model  and  simulated  using  the  detailed  model  data.  The 
results  in  Figure  IV.  14  show  the  opposite  trend  in  a  positive  identification  of  the  fault 
compared  to  previous  methods,  which  is  expected.  What  we  have  is  essentially  a  filter 
that  is  “tuned”  to  the  particular  fault  where  the  residuals  will  be  zero  when  the  fault 
occurs,  which  indicates  input/output  consistency  in  the  fault  condition.  All  five  states  of 
the  SSCM  model  reveal  a  constant  near-zero  mean  residual,  but  only  4  of  the  5  are  shown 
in  Figure  IV.  14. 
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Figure  IV.14:  SSCM  Switch  Closed  Failure:  Fault  Model  Residuals 


Finally,  we  go  back  to  the  nonlinear  input  estimation  problem  described  in  Figure 
III.  14.  The  SSCM  model  (0.19)  has  a  similar  model  equation  structure  as  the  power 
supply  model  (0.4)  and  so  its  control  signal  d  can  be  estimated  in  like  manner.  Recall 
that  in  the  failed  condition  where  the  output  voltage  is  far  greater  than  the  reference 
voltage,  the  SSCM  controller  is  continuously  sending  the  command  to  have  the  switch 
opened  as  seen  in  Figure  IV.  13.  Without  the  use  of  the  fault  model,  a  fault  indication  can 
still  be  achieved  if  a  UIO  scheme  is  set  up  such  that  it  estimates  the  command  signal  d  . 
If  the  control  signal  is  measurable,  then  a  residual  on  the  command  signal  can  indicate  the 
presence  of  an  SSCM  switching  fault  as  shown  in  Figure  IV.  15. 
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SSCM  Switch  Fault  (Closed):  Residual  on  SSCM  Control  Signal  D 


Figure  IV.15:  SSCM  Switch  Closed  Failure:  Residual  on  Control  Signal  d 

One  more  example  for  FDI  based  on  fault  modeling  is  presented.  Suppose  only 
the  port  power  supply  is  on  line  providing  power  to  all  three  zones  and  an  open  circuit  on 
the  bus  occurs  so  that  Zone  3  is  dead.  The  scenario  was  simulated  and  the  results  are 
shown  in  Figure  IV.  16.  The  bus  current,  is  shown  to  drop  from  the  no  fault  Mode  3 
condition  down  to  the  Mode  2  condition.  This  is  consistent  with  what  is  expected  since 
the  load  in  Zone  3  is  completely  dropped. 
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Figure  IV.16:  Bus  Open  Circuit  in  Zone  3:  Bus  Current 

We  use  the  three  bus  mode  models  shown  in  Figure  II. 6  and  Table  II.4  as  a  means 
to  detect  the  bus  open  circuit.  Using  the  inputs  and  output  measurements,  we  process  the 
data  through  three  Kalman  Filters  tuned  to  each  bus  mode  configuration.  The  model  that 
is  most  consistent  with  the  input  and  output  data  will  show  a  zero  residual.  The  observers 
are  based  on  bus  model  equations  (O.lO)-(O.l  1),  and  the  results  are  shown  in  Figure 
IV.  17.  From  the  residuals,  it  can  be  seen  that  the  system  began  in  Mode  3  operation  (i.e., 
all  zones  on  line)  and  then  suddenly  switched  to  zone  2,  indicating  an  open  circuit  in  zone 
3. 
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Figure  IV.17:  Bus  Open  Circuit:  3  Mode  Residuals 


To  conclude  this  section,  it  is  shown  that  FDI  using  fault  models  shows  good 
results.  But  like  all  model-based  methods,  they  are  dependent  on  the  accuracy  of  the 
model.  Fortunately,  these  models  have  been  proven  to  be  accurate. 
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V.  CONCLUSION 


A.  SUMMARY  OF  FAULT  DIAGNOSIS  RESULTS  TO  THE  DC  ZEDS 

MODEL 

Model-based  Fault  Diagnosis  is  an  excellent  method  for  validating  system 
input/output  consistency  and  to  detect  fault  indicators  in  the  residuals,  which  is  the 
fundamental  principle  of  analytical  redundancy.  However,  these  methods  are  only  as 
good  as  the  models  upon  which  they  are  based.  Fortunately,  it  has  been  shown  that  the 
DC  ZEDS  model  is  an  accurate  power  electronics-based  distribution  system  and  is  very 
similar  in  design  to  the  reduced  scale  hardware  model  implemented  in  the  Naval  Combat 
Survivability  (NCS)  Testbed  located  at  Purdue  University,  which  is  specifically  designed 
for  concept  studies  such  as  the  ones  in  this  work  [335].  The  DC  ZEDS  system  is  a 
reduced-scale  model,  yet  it  contains  enough  of  the  complex,  nonlinear,  dynamically 
interdependent  nature  to  make  it  quite  suitable  for  scholarly  research 
[56]  [327]  [328]  [245][246]  [335]  [285]. 

Two  models  were  used  side-by-side  for  the  results  of  this  research:  the  average 
value  model  and  the  more  accurate  detailed  model.  The  difference  between  the  two 
models  is  the  20  kHz  switching  dynamics  of  the  power  supply  and  SSCM.  The  average 
value  models  are  nonlinear,  but  the  switching  dynamics  are  removed  and,  through  the  use 
of  unknown  input  observers,  accurate  state  and  input  estimates  are  determined  for  robust 
residual  generation.  In  these  fault  scenarios,  most  of  the  FDI  results  aim  at  detecting  the 
source  of  a  fault  rather  than  studying  its  propagation  effects. 

Based  on  the  results,  the  average  value  model  accurately  matched  the  detailed 
model  for  the  fault  scenarios  studied  herein  -  even  in  the  transients.  It  is  acknowledged, 
however,  that  a  class  of  faults  caused  by  switching  dynamics  may  not  be  accurately 
modeled  (and  therefore  not  accurately  estimated)  using  the  average  value  model  and  so 
they  were  not  addressed. 

The  closest  one  of  the  fault  scenarios  came  to  studying  switching  dynamic  faults 
is  the  one  with  the  SSCM  switch  failure  in  the  closed  position.  Here,  the  nonlinear 
dynamics  are  essentially  removed  by  a  constant  connection,  and  a  model  of  the  fault  is 
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formed.  The  FDI  scheme  in  this  scenario  used  an  observer  based  on  the  fault  model  to 
generate  a  residual.  Similarly,  another  fault  model  scheme  presented  in  the  results  uses  a 
bank  of  Kalman  Filters  specifically  “tuned”  to  detect  an  open  circuit  in  various  places 
along  the  bus.  Again,  the  indication  of  this  fault  will  be  a  zero  in  the  residuals,  meaning 
that  there  is  input/output  consistency  in  the  fault  condition.  As  long  as  the  models 
accurately  represent  the  system  and  the  fault,  the  fault  model  FDI  scheme  works  quite 
well.  One  drawback  to  this  method  might  be  that  it  is  specifically  designed  to  detect  a 
particular  fault  while  other  FDI  methods  (using  a  structured  residual  set,  for  example)  are 
able  to  detect  a  broader  class  of  faults.  This  might  not  be  undesirable,  though,  if  several 
FDI  methods  are  applied  concurrently;  for  example,  one  FDI  method  might  be  used  for 
detecting  a  broader  class  of  faults  with  less  accuracy  and  the  other  method  might  be  used 
for  detecting  (with  greater  accuracy)  a  very  specific  fault  that  is  known  to  produce  grave 
consequences. 

For  a  more  general  FDI  scheme,  it  is  shown  that  model  consistency  for  the  overall 
large-scale  system  can  be  successfully  subdivided  into  modules  with  common  inputs  to 
and  outputs  from  connecting  modules.  In  this  manner,  classes  of  faults  (such  as  sensor 
faults  or  communication  faults  between  modules)  can  be  isolated  to  a  particular  module 
by  verifying  input/output  consistency  between  modules.  This  modular  FDI  approach  has 
the  effect  of  being  able  to  subdivide  into  arbitrarily  small  submodules,  provided  that  the 
inputs  and  outputs  are  common  to  the  connecting  modules. 


B.  SUMMARY  OF  POLYNOMIAL  UIO  RESULTS 

A  special  type  of  observer  that  can  be  made  insensitive  to  any  combination  of 
inputs  (provided  certain  conditions  are  satisfied)  is  used  extensively  for  the  FDI 
application.  It  is  shown  in  several  examples  that  a  structured  residual  set  of  unknown 
input  observers  can  accurately  detect  and  isolate  actuator  faults  and  in  some  cases  can 
estimate  the  fault  signal  itself.  Here,  two  UIO  methods  based  on  a  polynomial 
representation  of  a  system  are  developed  and  compared  to  the  subspace  projection  UIO 
method.  Since  these  methods  take  advantage  of  the  partial-state  representation  of  a 
system,  it  is  shown  that  in  some  cases  an  estimate  of  the  partial  state  can  still  be 

determined  even  if  the  UIO  rank  condition  is  not  satisfied,  which  is  a  necessary  condition 
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for  the  subspace  projection  UIO  method.  For  the  purpose  of  residual  generation,  a 
partial-state  representation  could  possibly  be  all  that  is  needed.  If  the  full-state  estimate 
is  required,  then  it  can  still  be  done  by  taking  derivatives  of  the  partial  state  -  a  process 
that  is  sensitive  to  noise. 

The  input  replacement  UIO  method  is  straightforward  to  implement,  but  is 
restricted  to  a  smaller  class  of  systems  where  higher  degrees  of  the  numerator  polynomial 
is  necessary.  One  way  to  ensure  the  higher  numerator  polynomial  degree  is  to  include 
output  measurements  that  are  dependent  on  the  states  and  inputs.  In  the  state  space 
realization,  this  implies  that  the  D  matrix  is  nonzero.  The  input  replacement  UI  method 
is  also  similar  to  the  input  inversion  problem. 

The  parametrization  UIO  method  is  comparable  in  performance  to  the  subspace 
projection  method  and  applies  to  a  larger  class  of  systems.  It  is  based  on  the  nonunique 
solution  of  the  Diophantine  equation,  where  the  particular  solution  we  choose  separates 
the  known  system  dynamics  from  the  unknown.  It  is  believed  that  this  concept  can  be 
extended  to  an  adaptive  polynomial  UIO  method  where  the  parameter  vector  is 
determined  adaptively  by  minimizing  the  output  residual.  However,  precise  conditions 
for  existence  of  a  solution  still  have  to  be  detennined.  In  some  applications  to  the  DC 
ZEDS  models,  the  parametrized  solution  of  the  UIO  provided  a  nonrealizable  estimator. 

Having  a  feed-through  D  matrix  in  the  system  of  equations  of  the  state  space 
realization  (which  implies  higher  order  numerator  polynomials  in  the  right  matrix 
polynomial  realization)  also  allows  for  input  estimation  without  having  to  take 
derivatives  of  the  partial  state.  Input  estimation  is  quite  useful  for  fault  diagnosis  as 
shown  in  several  examples  where  the  input  to  a  nonlinear  model  was  estimated.  It  is  also 
a  key  element  to  the  modular  approach  to  FDI  as  it  enables  each  module  to  estimate  both 
outputs  and  inputs  for  consistency  validation.  It  has  also  been  shown  that  in  some  cases 
an  estimate  of  the  fault  signal  itself  can  be  detennined. 
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C.  FUTURE  WORK 

For  follow-on  work,  other  non-model  based  fault  diagnosis  strategies  could  be 
developed  and  compared  to  these  results.  It  would  also  be  valuable  to  validate  the  fault 
diagnosis  schemes  using  the  Naval  Combat  Survivability  (NCS)  Testbed;  it  is  very 
similar  in  design  to  the  DC  ZEDS  model  described  in  Chapter  II  and  has  been  developed 
by  the  same  team  of  engineers  at  Purdue  University  [335].  Here,  real-time  FDI  on  an 
actual  hardware  DC  ZEDS  system  could  be  compared  with  simulation  results.  The  work 
would  require  modeling  and  simulation  on  a  validated  truth  model  using  Advanced 
Continuous  Simulation  Language  (ACSL). 

At  the  reconfiguration  level  shown  in  Figure  1.1,  raw  residual  data  is  analyzed  and 
decisions  are  made  about  positive  fault  identification.  From  the  innovations  determined 
in  Chapter  IV  on  the  stochastic  models  described  in  Chapter  II,  residual  evaluation 
techniques  can  be  applied  to  detennine  fault  probabilities  or  to  develop  an  intelligent 
fault  decision-making  process.  Methods  for  residual  evaluation  may  include: 


•  Generalized  likelihood  ratio  (GLR)  testing 

•  Fuzzy  logic 

•  Adaptive  thresholds  (time  and/or  freq  domain) 

•  Hypothesis  testing 

•  Multiple  Model  Adaptive  Estimation  (MMAE) 


In  this  systematic  manner  of  evaluating  residuals,  a  positive  fault  indication  can  then  be 
sent  to  a  control  reconfiguration  scheme  to  achieve  the  ultimate  goal:  continuity  of 
service  despite  combat  damage  involving  cascading  failures. 

Other  areas  for  FDI  research  could  be  to  introduce  multiple  simultaneous  faults  or 
conduct  similar  fault  simulations  on  the  DC  ZEDS  system  with  a  droop  a  mode  capability 
where  power  in  each  zone  is  brokered  between  the  port  and  starboard  busses. 
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In  [60],  authors  Chen  and  Patton  note  that  robust  FDI  is  still  an  “open  problem” 
where  the  most  significant  challenge  is  due  to  conflicting  criteria  in  model-based  robust 
residual  generation.  Specifically,  residual  generation  methods  must  be  sensitive  to  faults, 
yet  insensitive  to  modeling  uncertainties.  In  order  to  account  for  this  in  the  residuals,  a 
multicriteria  optimization  approach  can  be  investigated  where  fault  effects  are  maximized 
and  modeling  uncertainties  are  minimized.  Some  general  optimal  FDI  work  has  been 
studied  using  a  single  objective  function  codified  as  a  weighted  sum  of  objectives 
[59] [60].  In  many  of  these  cases,  however,  the  solutions  do  not  adequately  address  the 
original  problem  [322].  What’s  more,  FDI  is  most  often  considered  separately  from  the 
rest  of  the  system  design.  A  better  approach  might  be  to  develop  a  multicriteria  optimal 
solution  by  investigating  trade  offs  between  design  variables,  system  constraints,  control 
laws,  sensor  allocation,  and  fault  sensitivity. 
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APPENDIX  A:  IFAC  SAFEPROCESS  FAULT  DIAGNOSIS 

DEFINITIONS 


The  following  keywords  are  common  tenns  applied  to  the  field  of  Fault  Diagnosis  and 
Fault  Tolerant  Control.  They  were  developed  by  the  Fault  Detection,  Supervision,  and 
Safety  for  Technical  Processes  (SAFEPROCESS)  Technical  Committee  and  published  in 
Isermann  and  Balle,  1997;  [178].  The  definitions  are  widely  accepted  in  published  books 
on  the  subject  (Blanke  et.  al.,  2003;  [44],  Chen  and  Patton,  1997;  [60]). 


Active  fault-tolerant  A  fault-tolerant  system  where  faults  are  explicitly  detected  and 
control  system  accommodated.  Opposite  of  a  passive  fault-tolerant  system. 


Analytical 

redundancy 


Use  of  two  or  more,  but  not  necessarily  identical  ways  to  determine  a 
variable  where  one  way  uses  a  mathematical  process  model  in  analytical 
form. 


Availability 


Constraint 


Dependability 


Diagnostic  model 


Probability  of  that  a  system  or  equipment  will  operate  satisfactorily  and 
effectively  at  any  point  in  time. 

The  limitation  imposed  by  nature  (physical  laws)  or  man.  It  permits  the 
variables  to  take  certain  values  in  the  variable  space. 

A  form  of  availability  that  has  the  property  of  always  being  available  when 
required.  It  is  the  degree  to  which  a  system  is  operable  and  capable  of 
performing  its  required  function  at  any  randomly  chosen  time  during  its 
specified  operating  time,  provided  that  the  item  is  available  at  the  start  of 
that  period. 

A  set  of  static  or  dynamic  relations  which  link  specific  input  variables  -  the 
symptoms  -  to  specific  output  variables  -  the  faults. 


Discrepancy 


An  abnormal  behavior  of  a  physical  value  or  inconsistency  between  more 
physical  values  and  the  relationship  between  them. 


Disturbance 


An  unknown  (and  uncontrolled)  input  acting  on  a  system 


Error 


Deviation  between  a  measured  or  computed  value  (of  an  output  variable) 
and  the  true,  specified,  or  theoretically  correct  value. 


Fail  safe  The  ability  to  sustain  a  failure  and  retain  the  capability  to  make  a  safe  close¬ 

down.  A  system  where  the  occurrence  of  a  single  fault  can  be  determined 
but  not  isolated  and  where  the  fault  cannot  be  accommodated  to  continue 
operation. 
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Fail-operational 

Failure 

Failure  effect 

Failure  mode 
Fault 

Fault 

accommodation 
Fault  detection 
Fault  detector 
Fault  diagnosis 

Fault  estimation 
Fault 

identification 
Fault  isolation 

Fault  modeling 

Fault  propagation 
analysis 

Fault  tolerant 
control 


A  fault-tolerant  system  has  the  property  that  a  single  fault  does  not  develop 
into  a  failure  of  the  closed-loop  system.  A  Fail-operational  system  is  one 
that  the  performance  of  the  faulted  system  remains  the  same  (original 
objective  of  the  nominal  system  is  met). 

Permanent  interruption  of  a  systems  ability  to  perform  a  required  function 
under  specified  operating  conditions. 

The  consequence  of  a  failure  mode  on  the  operation,  function,  or  status  of  an 
item. 

Particular  way  in  which  a  failure  can  occur. 

Unpennitted  deviation  of  at  least  one  characteristic  property  or  parameter  of 
a  system  from  its  acceptable/usual/standard  condition.  [Blanke  et.  al.  adds: 

A  fault  is  the  occurrence  of  a  failure  mode  (Blanke,  2003)] 

A  correcting  action  (reconfiguration  or  a  change  in  the  operation  of  a 
system)  that  prevents  a  certain  fault  to  propagate  into  an  undesired  end- 
effect. 

Determination  of  faults  present  in  a  system  and  time  of  detection. 

An  algorithm  that  perfonns  fault  detection  and  isolation. 

Detennination  of  kind,  size,  location,  and  time  of  occurrence  of  a  fault. 

Fault  diagnosis  includes  fault  detection,  isolation,  and  estimation. 


Detennination  of  the  size  and  behavior  of  a  fault  over  time. 

Detennination  of  the  size  and  time-variant  behavior  of  a  fault.  Follows  fault 
isolation 

Detennination  of  kind,  location,  and  time  of  detection  of  a  fault.  Follows 
fault  detection. 

Detennination  of  a  mathematical  model  to  describe  specific  fault  effect. 


Analysis  to  detennine  how  certain  fault  effects  propagate  through  the 
considered  system. 

Aims  at  changing  the  control  law  so  as  to  cancel  the  effects  of  the  faults  or 
to  attenuate  them  to  an  acceptable  level. 
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Fault  tolerant 

system 

Hardware 

redundancy 

Incipient  fault 

Malfunction 

Monitoring 


Passive  fault- 
tolerant  system 

Perturbation 


Protection 


Qualitative 

model 

Quantitative 

model 

Recoverability 

Reliability 

Remedial  action 


A  system  where  a  fault  is  accommodated  with  or  without  performance 
degradation,  but  a  single  fault  does  not  develop  into  a  failure  on  subsystem 
or  system  level.  Fault  tolerant  system  has  the  property  that  a  single  fault 
does  not  develop  into  a  failure  of  the  closed-loop  system  (p.  9). 

Use  of  more  than  one  independent  instrument  to  accomplish  a  given 
function 

A  fault  where  the  effect  develops  slowly  e.g.  clogging  of  a  valve.  In 
opposite  to  an  abrupt  fault. 

An  intermittent  irregularity  in  the  fulfillment  of  a  system's  desired  function. 


A  continuous  real-time  task  of  determining  the  conditions  of  a  physical 
system,  by  recording  infonnation,  recognizing  and  indicating  anomalies  in 
the  behavior 

A  fault-tolerant  system  where  faults  are  not  explicitly  detected  and 
accommodated,  but  the  controller  is  designed  to  be  insensitive  to  a  certain 
restricted  set  of  faults.  Contrary  to  an  active  fault-tolerant  system. 

An  input  acting  on  a  system,  which  results  in  a  temporary  departure  from  the 
current  state. 

Means  by  which  a  potentially  dangerous  behavior  of  the  system  is 
suppressed  if  possible  or  means  by  which  the  consequences  of  a  dangerous 
behavior  are  avoided. 

Use  of  static  and  dynamic  relations  among  system  variables  and  parameters 
in  order  to  describe  a  system’s  behavior  in  qualitative  tenns  such  as 
causalities  or  if-then  rules. 

Use  of  static  and  dynamic  relations  among  system  variables  and  parameters 
in  order  to  describe  a  system’s  behavior  in  quantitative  mathematical  terms. 

Possibility  to  accommodate  the  fault  or  to  reconfigure  the  system  if  fault 
occurs. 

Ability  of  a  system  to  perfonn  a  required  function  under  stated  conditions, 
within  a  given  scope,  during  a  given  period  of  time. 

A  correcting  action  (reconfiguration  or  a  change  in  the  operation  of  a 
system)  that  prevents  a  certain  fault  to  propagate  into  an  undesired  end- 
effect.  Synonymous  to  fault  accommodation. 
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Residual 


A  fault  indicator,  based  on  a  deviation  between  measurements  and  model- 
equation-based  computations 


Safety 


Ability  of  a  system  not  to  cause  danger  to  persons  or  equipment  or  the 
environment 


Safety  system 


Electronic  system  that  protects  local  subsystems  from  permanent  damage  or 
damage  to  environment  when  potential  dangerous  events  occur. 


Sensor  fusion 
Supervision 

Symptom 

Threshold 


Integration  of  correlated  signals  from  different  [sources]. 

Monitoring  of  a  physical  system  and  taking  appropriate  actions  to  maintain 
the  operation  in  the  case  of  faults. 

A  change  of  an  observable  quantity  from  normal  behavior. 

Limit  value  of  a  residual's  deviation  from  zero,  so  if  exceeded,  a  fault  is 
declared  as  detected. 
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APPENDIX  B:  POWER  SUPPLY  AVM  TRANSFORMATION 


The  average  value  model  begins  with  several  reference  frame  transfonnations  of 
the  input  voltage  to  obtain  the  time-dependent  nns  voltage  of  the  input.  The  source  for 
the  transformation  is  shown  in  Krause  et  al  [220]. 

The  3-phase  ac  input  voltage  is  expressed  as 


vag=42Ecos(9g) 

(.87) 

vbg  =42Ecos[og-^- 

V  5  9 

(.88) 

vcg  =  V2£cos  f#g  +  ^- 

V  5  9 

(.89) 

and  the  arbitrarily  rotating  real  and  reactive  (q-d)  axis  representation  is  expressed  as 


Vqd0s 


=  Kv 


s  v abc  ? 


(.90) 


where 


Ks=- 

3 

cos  6g 

sin^ 

cos(<9,,  -2;r/3) 
sin(6(,  -2;r/3) 

cos(<9?  +2;r/3) 
sin(#e  +2^/3) 

1/2 

1/2 

1/2 

If  the  input  deviates  from  being  balanced,  then  the  third  of  3  outputs  of  the 
transformation  will  be  nonzero. 
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A  second  transformation  is  applied  to  the  arbitrarily  rotating  orthogonal  q-d  axis 
to  a  reference  frame  where  the  reactive  voltage  is  identically  equal  to  zero.  In  doing  so, 
the  transformation  calculates  the  rms  voltage  of  the  input,  y/2E . 


<1 

.Vfj 

0 

i  Kg^qdOs 


(.91) 


where 


a 


K 


COS  (<pg)  sin(^g)  0 
-sin(^g)  cos(^g)  0 
0  0  1 


and 


(Pg  =  tan 


0*. 


(.92) 


The  average  value  model  of  the  3 -phase  bridge  rectifier  is  developed  in  [Krause, 
et  al]  and  the  schematic  of  the  Port  Power  Supply  in  the  DC  ZEDS  model  is  shown  in 
Figure  1  along  with  its  Buck  Converter  controller.  The  input  voltage  expression  takes 
into  account  the  line-to-line  voltage  of  the  balanced  3 -phase  ac  input  and  represents  it 
with  the  rms  voltage.  The  3  phase  full  wave  rectifier  equipped  with  thyristors  allows  for 
full  control  of  the  output  dc  voltage  through  the  control  input  firing  angle  a  . 
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APPENDIX  C:  SUBSPACE  PROJECTION  UNKNOWN  INPUT 

OBSERVER 


The  purpose  of  Appendix  C  is  to  develop  the  subspace  method  Hui  and  Zak’s 
unknown  input  observer  from  a  slightly  different  approach  presented  in  [175].  The 
resulting  observer  design  allows  for  any  combination  of  inputs  to  be  made  insensitive  to 
the  observer  -  or  all  of  them  -  provided  that  certain  rank  conditions  are  met.  This  design 
freedom  makes  this  observer  type  useful  for  Fault  Diagnosis.  The  paper  does  not 
explicitly  address  output  measurements  that  are  dependent  on  the  inputs  (i.e.  D  ^  0 ), 
only  to  the  extent  that  the  input  can  be  subtracted  from  the  output  measurements  before 
the  UIO  is  implemented  and  then  added  back  to  obtain  an  output  estimate. 

The  class  of  dynamical  models  considered  herein  is  of  the  form 


x  =  Ax  +  Bu 

(.93) 

8 

ii 

(.94) 

where  state  vector  x  e  M"’ ,  input  u  e  R''” ,  and  output  y  e  R"'  and  matrices  A  ,  B  and  C 
are  completely  known  with  appropriate  dimensions.  For  a  system  where  at  least  one  of 
nu  inputs  is  unknown,  the  input  matrix  B  and  corresponding  input  vector  function  u  are 
partitioned  into  known  and  unknown  inputs  as 


Bu  =  [f?, 


where  the  first  n  known  inputs  are  separated  from  the  remaining  nu  =  nu  -  nlh  unknown 

inputs.  In  this  manner,  the  inputs  of  (  .93)  can  be  represented  with  separated  known  and 
unknown  inputs  with  the  same  output  equations  as 
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(.95) 


x  =  Ax  +  Bxux  +  B2u2 
y  =  Cx 


where  Bx  e  M”"1  and  B2  e  M”"2  .  It  should  be  noted  up  front  that  the  well-known  rank 
condition  must  be  satisfied  in  order  to  check  for  the  existence  of  the  observer  [230].  For 
the  observer  to  exist,  rank  ( CB2 )  =  rank  ( B2 ) ,  which  implies  that  there  must  be  at  least  as 
many  independent  outputs  as  there  are  unknown  inputs. 

The  Hui  and  Zak  UIO  algorithm  is  based  on  the  idea  that  we  want  to  project  the 
unknown  component  in  M"’  onto  the  subspace  of  known  state  components.  Then,  an 
observer  is  constructed  within  this  subspace  to  estimate  the  entire  state.  Let  P  be  a 
projection  matrix  such  that 


P  =  MC 


(.96) 


where  P  e  W  and  M  e  M”'  .  The  projection  matrix  has  two  important  properties.  First, 
P  is  symmetric  ( PT  =P).  Second,  the  square  of  P  is  itself  ( P 2  =P).  The  second 
property  also  implies  that 


P"  =P  Vw  >  1 


since  any  positive  integer  exponent  n  can  either  be  divided  by  two  with  a  remainder  of  0 
or  1.  When  expanding  the  exponents  of  P  the  result  will  always  be  a  series  of 
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multiplications  of  the  projection  matrix  raised  to  the  first  or  second  power,  which  will 
always  reduce  to  the  second  property  of  the  projection  matrix.  Using  the  projection 
matrix,  the  state  vector  can  be  separated  as 


x  =  (l  -P)x  +  Px 


(.97) 


Notice  that  the  expression  ( I  —  P )  is  a  projection  matrix  itself  because  it  is  symmetric 
and 


(I-P)2  =  I2-2P  +  P2 
=I-2P+P 


=(i-p). 


We  can  also  conclude  that 


P(I -P)  =  P-P2  =0. 


Substituting  the  output  equation  ( .94)  and  ( .96)  into  equation  ( .97)  we  get 


x  =  q  +  My 


(.98) 


where 


q  =  ( /  -  MC  )  x . 


(.99) 
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The  projection  matrix  I -MC  is  constant  so  the  derivative  of  ( .99)  is  simply 


q  =  (/ -MC)x 


(.100) 


Notice  that  q  can  be  viewed  as  the  projection  (not  necessarily  orthogonal)  of  x  onto  the 
subspace  of  known  components  in  M"’  .  If  the  subspace  state  q  can  be  estimated,  then  an 
estimate  of  the  entire  state  x  can  be  determined  by  equation  (  .98)  since  the  output 
measurement  vector  y  is  known.  To  detennine  M  ,  the  constraint 


(I-MC)B2=  0 


(.101) 


must  be  met,  which  essentially  removes  the  influence  of  the  unknown  input  in  (  .95).  An 
observer  is  then  designed  within  the  subspace.  Provided  that  the  rank  condition 
rank(CB2 )  =  rank (  B2 )  is  met  (Kudva  el  al  [230]),  the  known  state  subspace  exists  and 
the  matrix  M  is  detennined  to  be 


M  =  &((CR)'  +//„(/,  -(CA)(CA)'))  ( .102) 


where  the  superscript  f  denotes  the  Moore-Penrose  pseudo-inverse  operation  and 
H0  is  [nu  xn)  is  a  design  parameter  matrix  introduced  in  [175]  (default  is  zero).  To 

determine  the  dynamic  state  space  model  for  the  subspace  state  q  ,  we  substitute  equation 
(  .95)  into  ( .100)  to  get 

q=(l-MC)A(q  +  My)  +  (l-MC)Blu1+(l-MC)B2u2  (  .103) 
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Recall  that  M  is  chosen  such  that  the  constraint  (  .101)  is  met  and  so  the  last  tenn 
in  equation  -  the  term  containing  all  unknown  inputs  -  is  zero.  The  resulting  state  space 
equation  is  of  the  form 


q  =  Aq  +  Bu  (-104) 

v  =  Cq  (.105) 


where 


A  =  (I-MC)A 
B  =  (I-MC)[AM  Bx] 

u  =  '  and 

Wj_ 

y  =  y-CMy 


Notice  the  new  input  u  is  the  vector  function  of  output  measurements  and  known 
inputs,  which  is  fully  known.  From  here,  it  can  be  shown  that  a  standard  proportional 

gain  observer  can  be  applied  to  equations  (  .  104)-(  .105)  provided  the  pair  (z(,c)  is  fully 
observable.  Note  that  it  is  possible  that  the  original  system  might  be  fully  observable 
while  the  pair  is  not.  Choosing  the  closed  loop  gain  matrix  L  to  be 


L  =  (I-MC)L, 


where  L  is  the  observer  gain  matrix  chosen  for  the  original  system  in  (  .93)  and  (  .94). 
Finally,  the  estimate  of  the  subspace  state  q  is  determined  to  be 
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q  =  Aq  +  Bu  +  L[y - Cq^ . 


(.106) 


Now  to  show  convergence  of  the  subspace  method,  choose  the  state  error  to  be  e  =  q-q 
so  that  the  error  dynamics  are  determined  by  the  equation 


e  =  (A-LC)e. 

/ 


Knowing  that  A  =  (/ - MC) A  and  L-[1-MC)L  defined  in  equations  (  .104)  and 

(  .105),  the  error  dynamics  can  be  expressed  as  the  error  of  the  original  observer 
projected  onto  the  known  system  subspace  as 


e  =  (l-MC){A-LC)e. 


(.107) 


Looking  at  equation  (  .107),  it  is  not  readily  apparent  that  the  error  converges  to  zero. 
With  a  few  observations,  though,  it  can  be  shown  that  it  does.  First,  pre-multiply 
equation  ( .107)  by  the  projection  matrix  MC  ,  which  yields 


MCe(t )  =  0  Vf 


(.108) 


since 


MC(I-MC)  =  MC-(MC)2  =  MC  -  MC  =  0 , 
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which  implies  that 


MCe(t )  =  constant  Vt . 


Next,  pre-multiply  the  projection  matrix  MC  on  q  in  equation  ( .99)  ,  yielding 


MCq(t)  =  0  Vf 


for  the  same  reason  as  in  (  .108).  Choose  the  initial  condition  on  the  estimate  q  in  the 
expression  as 


MCq(  0)  =  0, 


then 


MCe(t )  =  0  VC 


For  simplicity,  define  a  projection  matrix  P'  and  system  matrix  A'  as 


P'  =  I  -  MC  and 
A'  =  A-LC 


then  the  expression  for  the  error  equation  in  ( .107)  can  be  expanded  as 
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P'e(t)  =  P' 


=  P' 


I  +  P'A't  + 


P'2A'2t 2  A 

-  +  ... 


2! 


.(0) 


,  A'2‘ 2  A 

I  +  At  + - +  ... 


2! 


.(0) 


J 


=  P'eAte(  0) 


Now  choose  L  so  that  the  desired  plant  characteristics  A'  =  (A-LC)  converge  to 
zero  so  that  (/  -MC)e(t)  — >  0  .  Putting  it  all  together,  it  can  be  seen  that  the  error  e(t) 
approaches  zero  because 


(/  -  MC )  e{t )  =  e(t )  -  MCe{t )  — »  0 

e(t)  — >  0. 

Finally,  since  it  has  been  shown  that  the  estimate  q  of  the  known  component 
subspace  state  converges  to  the  actual  state  q ,  an  estimate  x  for  the  whole  state  can  be 
detennined  based  on  the  relationship  between  q  and  x  in  equation  ( .98)  as 


x  =  q  +  My  . 
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APPENDIX  D:  POLYNOMIAL  MATRIX  DEFINITIONS  AND 

THEOREMS 


The  purpose  of  Appendix  D  is  to  provide  definitions  to  polynomial  matrices  for 
which  the  multivariable  unknown  input  observer  method  is  based.  Most  of  the 
definitions  herein  are  documented  in  [61], [200],  and  [229]. 


Definitions 

Definition  D.l  (field).  A  field  consists  of  a  set,  denoted  by  F ,  of  elements  called  scalars 
and  two  operations  called  addition  "  +  "  and  multiplication  the  two  operations  are 
defined  over  IF  such  that  they  satisfy  the  following  conditions: 

1.  To  every  pair  of  elements  a  and  ft  in  F ,  there  correspond  an  element  a  +  f  in 
F  called  the  sum  of  a  and  ft ,  and  an  element  a  ■  f)  or  af  in  F ,  called  the 
product  of  a  and  (3 . 

2.  Addition  and  multiplication  are  respectively  commutative:  For  any  a,  f  in  F , 

a  +  f  =  f  +  a  a-p  =  p-a 

3.  Addition  and  multiplication  are  respectively  associative:  For  any  a,  P,  y  in  F , 

(a  +  p)  +  Y  =  (a-p)-y  (a  ■  p)-y  =  a  \p -y) 

4.  Multiplication  is  distributive  with  respect  to  addition:  For  any  a,P,y  in  F , 

a\p  +  y)  =  (a-p)  +  (a-y) 

5.  F  contains  an  element,  denoted  by  0,  and  an  element,  denoted  by  1,  such  that 
a  +  0  =  a,  1-  a  -  a  for  every  a  in  F . 

6.  To  every  a  in  F ,  there  is  an  element  [3  in  F  such  that  a  +  P  -  0 .  The  element 
P  is  called  the  additive  inverse. 
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7.  To  every  a  in  T  which  is  not  the  element  0,  there  is  an  element  y  in  T  such 
that  a  •  y  =  1 .  The  element  y  is  called  the  multiplicative  inverse. 

Note  1 :  The  symbol  M  and  C  is  used  to  denote  the  field  of  real  numbers  and  the  field  of 
complex  numbers,  respectively. 

Note  2:  The  set  of  integers  and  the  set  of  polynomials  do  not  form  a  field  because  they 
have  no  multiplicative  inverse  defined  within  the  field. 


Definition  D.2  (ring).  More  precisely:  a  commutative  ring  with  multiplicative  identity. 
A  set  with  all  properties  of  a  field  except  for  property  7  in  Definition  D.  1 . 


Note:  The  set  of  integers  and  the  set  of  polynomials  with  real  coefficients  form  a  ring. 


Definition  D.3  (linear  space  over  a  field).  A  linear  space  over  a  field  T ,  denoted  by 
(hi, IF ) ,  consists  of  a  set,  denoted  by  Tt,  of  elements  called  vectors,  a  field  T ,  and  two 

operations  called  vector  addition  and  scalar  multiplication.  The  two  operations  are 
defined  over  hi  and  F  such  that  they  satisfy  all  the  following  conditions 

1.  To  every  pair  of  vectors  xl  and  x2  in  hi ,  there  corresponds  a  vector  xl  +x2  in 
hi ,  called  the  sum  of  x3  and  x2 . 

2.  Addition  is  commutative:  For  any  x, ,  x2  in  hi ,  x,  +  x2  =  x2  +  x, . 

3.  Addition  is  associative:  For  any  xt,  x2,  x3  in  hi,  (x,  +  x2 )  +  x3  =  x,  +  (x2  +  x3 ) 

4.  hi  contains  a  vector,  denoted  by  0,  such  that  0  +  x  =  x  for  every  x  in  hi . 

The  vector  0  is  called  the  zero  vector  or  the  origin. 

5.  To  every  x  in  hi ,  there  is  a  vector  x  in  hi ,  such  that  x  +  x  =  0  . 

6.  To  every  a  in  T ,  and  every  x  in  hi,  there  corresponds  a  vector  ax  in  hi 
called  the  scalar  product  of  a  and  x . 
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7.  Scalar  multiplication  is  associative:  For  any  a ,  P  in  T  and  any  x  in  Ft, 
<x(/?x)  =  (<x/?)x . 

8.  Scalar  multiplication  is  distributive  with  respect  to  vector  addition:  For  any  a  in 
T  and  any  x, ,  x2  in  Ft,  a  (xj  +  x2 )  =  axx  +  ax2 . 

9.  Scalar  multiplication  is  distributive  with  respect  to  scalar  addition:  For  any  a , 
P  in  F  and  any  x  in  FI,  (a  +  P)  x  =  ax  +  /?x . 

10.  For  any  x  in  Ft ,  1  x  =  x ,  where  1  is  the  element  1  in  T . 


Note  1:  If  J~=  R ,  then  is  called  the  //-dimensional  real  vector  space.  For 

brevity,  the  field  of  real  numbers  is  implied  and  therefore  is  denoted  as  M" . 

Note  2:  [  p]  is  the  set  of  all  polynomials  of  degree  less  than  n  with  real  coefficients. 

Note  3:  If  a  field  is  replaced  by  a  ring  E(. ,  then  is  called  a  module  over  the  ring. 

A  n  x  1  or  lx//  polynomial  vector  can  be  considered  as  an  element  of  the  rational  vector 
space  (M"(/?),IR(//)) ,  or  an  element  of  the  module  (M”[//],E[//]) .  A  set  of  polynomial 
vectors  is  linearly  independent  over  the  field  R(/?)  if  and  only  if  the  set  is  linearly 
independent  over  the  ring  R[//] . 


Definition  D.4  (normal  rank).  The  normal  rank  of  a  polynomial  matrix 
A(p)  e  M"!X”  [//]  is  defines  as  e.g. 

1.  The  number 


rank  A(p)  =  max  rank  A(p) 

ssC  v  7 

2.  The  number  of  linearly  independent  columns  (or  rows)  of  A(p) 

3.  The  number  of  invariant  polynomials  of  the  Smith  form  of  A(p) 
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4.  The  dimension  of  the  space  spanned  by  the  rows/columns 


Definition  D.5  (rank  of  a  rational  matrix).  The  rank  of  a  matrix  M  ( /; )  e  M"'x"  ( p)  is 
defined  as  e.g. 

1.  The  number  of  linearly  independent  rows/columns  with  rational  coefficients 

2.  The  (normal)  rank  of  N  ( p )  where  N  ( p )  is  the  numerator  in  a  left  or  right  MFD 
of  M  (p) 

3.  The  dimension  of  the  space  spanned  by  the  rows/columns 

Definition  D.6  (matrix  pencil).  A  matrix  pencil  is  a  polynomial  matrix  M  ( p )  of  degree 
one.  It  is  often  written  as 

M[p)  =  pE  +  F 

where  E  and  F  are  constant  matrices  and  p  is  a  differential/difference  operator  in  the 
Laplace/Z  domain  in  continuous/discrete  time. 

Definition  D.7  (regular  matrix  pencil).  A  matrix  pencil  M  ( p )  is  said  to  be  regular  if 
M  ( p)  is  square  and  full  (normal)  rank. 


Definition  D.8  (controllability  index).  The  fact  that  the  controllability  matrix 

C  =  [b  AB  —  A”~1B~\ 

has  rank  n  means  that  there  are  n  linearly  independent  columns  in  C  .  It  could  be  that 
n  such  columns  in  the  partial  controllability  matrix 

Cq=\_B  AB  —  Aq-'B],  1  <  q  <  n 
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can  be  found  and  that  the  smallest  q  is  denoted  as  p .  The  controllability  index  p  of 
the  pair  (A,  B)  is  said  to  be  the  smallest  q  such  that  Cq  has  rank  n. 


Definition  D.9  (unimodular  matrix).  A  square  polynomial  matrix  M{p)  is  called 
unimodular  if  its  determinant  is  nonzero  and  independent  of  p  . 


Definition  D.10  (right/left  divisor,  right/left  multiple).  Consider  the  polynomial 
matrix  equation 

A(p)  =  B(p)C(p) 

where  A(p),  B(p),  and  C(p)  are  polynomial  matrices  of  appropriate  orders. 

1.  C(p)  is  a  right  divisor  of  A(p ) 

2.  B(p)  is  a  left  divisor  of  A(p) 

3.  A(p)  is  a  left  multiple  of  C{p ) 

4.  A(p)  is  a  right  multiple  of  B(p) 


Definition  D.ll  (greatest  common  right  divisor  (gcrd)).  A  square  polynomial  matrix 
R(p)  is  a  greatest  common  right  divisor  (gcrd)  of  N( p)  and  D(p)  if  R(p)  is  a  common 
right  divisor  of  N(p)  and  D(p )  and  is  a  left  multiple  of  every  common  right  divisor  of 
N(p)  and  D(p). 


Definition  D.12  (greatest  common  left  divisor  (geld)).  A  square  polynomial  matrix 
Q(p)  is  called  a  greatest  common  left  divisor  (geld)  of  A(p)  and  B(p)if  Q(p)  is  a 
common  left  divisor  of  A(p)  and  B(p)  and  Q(p)  is  a  right  multiple  of  every  common 
left  divisor  Qfp )  of  A( p)  and  B(p). 
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Definition  D.13  (right  coprime).  If  a  square  polynomial  matrix  R{p )  is  a  gcrd  of 
N(p )  and  D(  p)  and  R(p)  is  a  unimodular  matrix,  then  N(p)  and  D(p )  are  said  to  be 
right  coprime. 


Definition  D.14  (left  coprime).  If  a  square  polynomial  matrix  Q{p )  is  a  geld  of  A(p) 
and  B(p)and  Q(p )  is  a  unimodular  matrix,  then  A(p)  and  B(p)  are  said  to  be  left 
coprime. 


Definition  D.  15  (Sylvester  matrix).  Let  W(p)  be  a  polynomial  matrix  of  degree  d  such 
that 

W{p)  =  ±Wtp' 

i= 0 

with  Wd  0  .  Then  the  qth  order  Sylvester  matrix  for  W(p)  is  defined  as 

W0  W]  ...  Wd  0  0  0  " 

0  W0  Wx  •••  Wd  0  0 

0  0  0 

0  0  0  w0  wx  •••  wd 


S  (W  (p),q)  = 


Definition  D.16  (characteristic  polynomial).  The  characteristic  polynomial  of  a  proper 
rational  matrix  G(p)  is  defined  as 

1.  The  least  common  denominator  of  all  minors  of  G(p) .  The  degree  of  the 
characteristic  polynomial  is  defined  as  the  McMillan  or,  simply,  the  degree  of 
G(p)  and  is  denoted  by  SG(  p) . 

2.  Consider  the  proper  rational  matrix  G(p )  factored  as 

G(p)  =  N  (p)D\p)  =  Dfp)N(p) 
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where  N(p )  and  D(p )  are  right  coprime,  and  N(p)  and  D(p )  are  left  coprime.  Then 
the  characteristic  polynomial  of  G(p)  is  defined  as 

det  D(p)  or  det  D(p) 

and  the  degree  of  G(p)  is  defined  as 

5G{p )  =  S(  det  D(p))  =  <7(  det  D(p)) 


Definition  D.  17  (Right  polynomial  fraction).  Every  qxp  proper  rational  matrix  G(p ) 
can  be  expressed  as 

G(p)  =  N(p)D\p) 

where  N{ p)  and  D{p)  are  qxp  and  pxp  polynomial  matrices,  respectively.  G(p) 
in  this  form  is  called  a  right  polynomial  fraction  or,  simply,  right  fraction. 


Definition  D.18  (Left  polynomial  fraction).  The  dual  expression  of  G(p )  in  Definition 
D.16  where  the  same  proper  rational  matrix  can  be  equally  expressed  as 

G(p)  =  D-1(p)N(p ) 

where  D(p)  and  N(p)  are  qxqand  qxp  polynomial  matrices,  respectively.  G( p)  in 
this  form  is  called  a  left  polynomial  fraction  or,  simply,  left  fraction. 


Definition  D.19  (column/row  reduced).  A  nonsingular  polynomial  matrix  M{p)  is 
column  reduced  if 

£(det  M(p)}  =  sum  of  all  column  degrees 


It  is  row  reduced  if 


8 (det  M(p ))  =  sum  of  all  row  degrees 
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Definition  D. 20  (Bezout  Equation).  Consider  two  polynomial  matrices  A(p),B(p)  and 
four  polynomial  matrices  P(p),Q(p),R(p),S(p )  such  that 


P(p ) 

Q(p) 

A(p) 

G(p) 

R(p) 

S(p)_ 

B(p)J 

0 

then  the  following  can  be  said: 

1.  A(p),B(p)  are  left  coprime  if  and  only  if  G(p)  is  unimodular 

2.  R(p)  and  S(p)  are  right  coprime 

3.  B(p)Afp)  =  -Sl(p)R(p) 

This  implies  that  if  G( p)  is  unimodular,  then  right  and  left  polynomial  fractions  of  the 
same  proper  rational  matrix  are  equivalent  as 

H(p )  =  B(p)Afp)  =  Afp)B(p) 

where 

A(p)  =  -S(p) 

B(p )  =  R(p). 


Definition  D.21  (Diophantine  Equation).  Given  three  polynomial  matrices 
A(p),B(p),C(p)  that  are  mutually  right  coprime  and  deriving  equations  in  a  Bezout 
Equation  form,  define  the  Diophantine  Equation  as 

C(p)  =  E(p)A(p)  +  F(p)B(p) 

from  the  top  row  of  the  Bezout  equation  and  a  left  coprime  representation  such  that 

R{p)A{p)  +  S(p)B(p)  =  0 

from  the  bottom  row  of  the  Bezout  equation.  A  solution  E0(p),F0(p )  exists  such  that  all 
solutions  (nonunique)  can  be  expressed  as 

E(p)  =  Efp)  +  W  (p)R(p) 

E(p)  =  F0(p )  +  W(p)S(p) 
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with  W(p )  an  arbitrary  matrix  of  appropriate  dimensions. 

Definition  D.22  (polynomial  degree  operator).  Given  a  polynomial  matrix  M(p) , 
define 

8cjM{p )  =  degree  of  i‘h  column  of  M(p) 

8rjM(p )  =  degree  of  ith  row  of  M(p). 

Denote  8ci  the  column  degree  and  8n  the  row  degree. 

Definition  D.23  (right  polynomial  transfer  matrix  and  control  canonical  form) 

Let  a  minimal  realization  of  the  system 

px(t)  =  Ox(t)  +  Tu(t)  (-109) 

y(t)  =  Cx(t)  +  Du(t )  (.110) 

where  x(t)eRn%  u(t)  e  R"“ ,  y(t)  e  E'iy  and  (O,  T,C,  D)  are  of  appropriate  dimensions. 
Because  the  system  is  a  minimal  realization,  an  equivalent  representation  of  (  .109)- 
(  .110)  is  the  right  polynomial  transfer  matrix  fonn  as  G(p)  =  B(p)A  '(p) .  Let 
jul,ju2,...,pn  be  the  column  degrees  (i.e.  the  maximum  degree  of  the  polynomials  in  each 
column)  of  the  matrix  A{p)  and  p  is  the  differential  operator  in  continuous  time  or  the 
time  shift  operator  in  discrete  time.  If  z(t )  is  defined  as  the  partial  state  [200]  (or 
sometimes  called  the  pseudo  state  [61]),  then  it  is  easy  to  see  that  it  can  be  related  to  the 
input  u(t)  and  output  y(t)  as 

u(t)  =  A(p)z(t)  (-111) 


y(t)  =  B(p)z(t) . 
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(.112) 


The  partial-state  equations  in  (  .1 1 1)-(  .112)  can  be  expressed  in  state  space  form  by  first 
defining  the  matrices 


^Ap)  =  diag  pM' , p"2,..., €  Rn"™“ 


(.113) 


Lip)  =  diag 


1 

P 


1 

P 

ito- 


1 

P 


Mn»- 1 

p  % 


if? 


(A  +j“2  +~+tlnu  )xnu 


(.114) 


and  expressing  the  polynomial  matrices  of  ( .  1 1 1)-(  .  1 12)  as 

A(p)  =  Al£/l(p)  +  AlL/l(p)  (.115) 

B(p)  =  BhZ/l(p)  +  BlL/l(p)  (.116) 


where  ( Ah ,  Bh )  are  matrices  of  high  degree  coefficients  whose  dimensions  are  nu  x  nu 
and  nv  x  nu ,  respectively,  and  ( 4, ,  B, )  are  the  matrices  of  all  lower  order  coefficients 
with  dimension  nu  x  ns  and  n  x  /?a  ,  respectively.  The  partial-state  vector  with  dimension 
ns  =  //j  + ...  +  pn  ,  can  now  be  expressed  as 


z(t)  = 


z\(t),  _1Zi(t)  :  z2(t), ...  ,pM2^zx(t) : ... \zn  (t),...,pM"‘rAn  (t)  (.117) 


can  now  be  expressed  as 


z(t)  =  Lfl{p)z(t). 


(.118) 
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Substitute  (  .1 15)-(  .116)  into  (  .1 1 1)-(  .112)  and  assume  Ah  is  nonsingular  to  obtain  the 
state  space  realization  of  the  original  system  in  (0.32)  and  (0.33)  as 


pM'zx(t) 


a"x..  (0 


=  -Ah'A,z(t)  +  Ah'u(t) 


y(t)  =  (Bi~  BhAhAi )  z(0  +  BhA;xu(t) 


(.119) 


(.120) 


It  can  be  shown  [61]  that  the  ensuing  state  space  structure  of  (  .1 19)-(  .120)  is  in 
controllable  canonical  form  for  which  the  nonsingular  matrix  T  can  be  found  for  MIMO 
systems  using  a  method  such  as  the  one  published  in  (Zak,  2003;  [386]). 


Theorems 

Note:  Theorems  herein  are  stated  without  proof.  Source  information  is  provided. 


Theorem  D.l  ([61];Theorem  G-6).  A  square  polynomial  matrix  is  unimodular  if  and 
only  if  its  inverse  is  a  polynomial  matrix 


Theorem  D.2  ([61];  Theorem  7.1).  The  controllable  canonical  form  is  observable  if  and 
only  if  D(p)  and  N(p)  are  coprime. 
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Theorem  D.3  ([61];  Theorem  7.M2).  A  state  equation  ( A,B,C,D )  is  a  minimal 
realization  of  a  proper  rational  matrix  G(p )  if  and  only  if  the  ( A,B )  is  controllable  and 
the  ( A,C )  is  observable  or  if  and  only  if 

dim [A)  =  5G{p) . 


Theorem  D. 4  ([61];  Theorem  7.8).  Let  N{p)  and  D{p )  be  qxp  and  px p  polynomial 
matrices,  and  let  D(p)  be  column  reduced.  Then  the  rational  matrix  N(p)D  fp)  is 
proper  (strictly  proper)  if  and  only  if 

8ciN{p)  <  SciD(p )  ( 8ciN{p )  <  SdD(p )) 

for  i  =  1,2,..., p  . 


Theorem  D.5  ([61];  Theorem  7.M3).  All  minimum  realizations  of  G(p)  are  equivalent. 
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